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Résumé. — La fréquence d’une oscillation fondamentale varie avec l'amplitude de l’agitation 


thermique. Les oscillations secondaires sont de fréquences [v, + vil; [v, + ve + vil, .. 


+ 5 Vs Va, 


Vs, --., étant les fréquences des oscillations fondamentales. Elles prennent des amplitudes propor- 
tionnelles aux produits des amplitudes des oscillations fondamentales. 


Abstract. — The frequency of a fundamental oscillation varies with the amplitude of the thermal 


agitation. 


Secondary oscillations have the frequencies [v, + vel ; [vi + vo + vil, …. 
V3, --., being the frequencies of fundamental oscillations. 


OV Ve, 


They take amplitudes proportional 


to the products of the amplitudes of the fundamental oscillations. 


V. Seconde analyse. 


En seconde analyse (équations 37) les oscillations 
fondamentales sont déterminées par les équations 


. M MD 6 m 
Us + 2 6 Ë Uk Sr — E sr: (39) 
(a 


À. Seconde opération (1). — Les forces élémen- 
taires, constantes pendant des temps limités, com- 


prises dans la force Pa 5 (32), sont de deux sortes 
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ja = HP EEE D af 


X cos 27[S(p — q) — qug.sr + Oiy.sr] (40) 
et 
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St € re 
x cos 27[{S(p + q) — oxp,sr— piy,sr + 0— sr — Os]. (41) 


(2) Rappelons que la première opération qui constitue la 
première analyse consiste en la résolution des équations 
séculaires (27), obtenues en faisant état de forces de rappel 
conformes à la loi de Hooke. 
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Développée par une seule oscillation fonda- 
m FRE, 

mentale Sr, toute force 55 . (40) se répète 
sur tous les atomes en même position dans le 
. . . m : . or T: 
motif cristallin : Jo ç. — Josa.ss. Son intensité 
augmente ou diminue chaque fois que l’oscil- 
lation Sr gagne ou perd un quantum d’énergie Avs.. 
Mais ces variations d’intensité diminuent, en gran- 
deur relative, quand la température s’élève. Aux 
moyennes et aux hautes températures les 


forces fo? restent sensiblement constantes au 
cours du temps. 
m . 
Chaque force fo5 8 (41) est produite par 


deux oscillations fondamentales S+ et — Sr, de 
même fréquence : v_s — vg, mais pilotées par 
deux vecteurs d'onde, $ et — $, égaux et oppo- 
sés (29-30). En conséquence, les forces [TER 
développées par le même couple d’oscillations fon- 
damentales, Sr et — S+, prennent des intensités 
différentes sur les différents atomes : elles varient 
avec m et j. 

Toute oscillation fondamentale S+ change de 
phase temporelle, 2x04., à chaque fois qu’elle prend 
ou cède un quantum d’énergie kvs.. De la sorte, les 
deux phases 04. et 0. varient. Et leurs variations 
étant indépendantes, la différence 2r(0_ 57 — 05) 
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prend uniformément au cours du temps toutes les 


grandeurs entre 0 et 2x. Donc les forces fes Fu 
CM Re ENER DO SE) sontroullensen 
moyenne dans le temps. Et comme leurs phases 
temporelles 2r(8_. — 04.) oscillent en outre, indé- 
pendamment les unes des autres, leur résultante, 

m / m 
fre GS # : & fe ST, —S7) à 
reste sensiblement nulle à tout instant. Ainsi, les 
forces Re n’apportent pas de composantes 
permanentes aux tensions thermiques. Je n’en ferai 
donc plus état, ni de toutes les forces semblables, 
développées par deux oscillations harmoniques de 
même fréquence et de vecteurs d’onde égaux el 
opposés, constantes pendant des durées variables 
mais nulles en moyenne dans le temps. 


Soit ge (t) la force obtenue en retranchant de la 


m # 2 LS 
force 5 (32) toutes les forces élémentaires 
m m 
Îo2 Le et fre Fes COnStantes pendant des temps 
limités. 
Les équations 
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admettent la solution 
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PU 


jZ 
p’.sres est le module et 2r(@,,5.15.") l'argument 
œ 
de C'srtsr. 


Les coordonnées vectorielles Cuve U = 1 
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2, ..., g;x — 1, 2, 3) sont déterminées par les 
équations séculaires : 


: o :% ik, ç 7 7x0 
(ose + og)? D ,sT+ 5/7 -È LABS + S) CE sr+ sr 
e { 
; mp 
| DA 
22 [E PT exp. 27 Slim —p) 
kit 6y La Vu; x LU 


exp.127 S’(in — a| tb se CV grer: 


pi ik, ç@ + 7 
(@,$7—0 5/7)" GPAPE à Legs Té S") CP ge 8x 
CE 


PaV EE; Uy LU 


3) 
= [E nn exp.127 Sim — p) 
kel RY 


s ÿ/ ” B Ty 
exp.—127 S'{in — q)| CE sr se 


ROUE PPT ER PRE een, (EE) 

Les fréquences vs; prises en compte dans les 
équations (45) sont celles qui sont déterminées par 
les équations de d’Alembert (39) de la seconde 
analyse. À vecteur d’onde $ égal, elles diffèrent des 
fréquences v$- données par les équations sécu- 
laires (27) de la première analyse. Je suppose que les 
vecteurs d’ondes fondamentaux, S, sont déter- 
minés. De la sorte, les fréquences vs. données par 
les équations (27) et par les équations (39) forment 
deux suites discrètes qui sont distinctes. La proba- 
bilité d’une égalité 


los ‘a &,8/z"| —= O,$4$/,7 


€) 


= PURES ) D] 
PAT T SD es + OU 


(45) 
est nulle ou infime. La matrice 
(o,sr + ©,87)? 1 — F(s+s’) 


est toujours, ou presque toujours régulière (2). 
Posons 


Gsrksr — [(o,sr + ©,82)? 1 —Tists)] 1. 
Nous avons (°) : 


nv pa 
wi l 


; ; k /! 
CP pepe = > PAST DOS b EST 
on 


Ke BY pa vu ty tu 


< eXp.127 S(m — p) exp.127 S’(m al cxê sr sp. 


(?) J'est la matrice unité d’ordre 3g X 3g. 

(*) Les oscillations secondaires que je définis sont celles 
qui sont produites par les oscillations fondamentales en 
états stationnaires. 


Je suppose déterminées les fréquences vs- et, par consé- 
quent, les énergies, 


(n + 1/2) hvsz, DU TP 0 
des oscillations fondamentales en états stationnaires. Mais 
cela n'implique pas que les quantums d'énergie kvs- perdus 
ou gagnés par les oscillations fondamentales soient rigou- 
reusement définis. 

Les états stationnaires des oscillations fondamentales ne 
sont pas durables. La probabilité de l’existence d’une oscil- 
lation fondamentale S+ dans un état stationnaire donné — 
associé à l’énergie (n +- 1/2) hkvs (n constant) — décroît au 


ni he 0 “fn AS à. 


nt pl Cd HOPPER 
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… Pour obtenir les coordonnées vectorielles 4. 
il suffit de remplacer, dans l’expression précé- 


dente (47), GP gr+sr par GC. 6, S'par—S’ et 
Cor" par CP se (°). 


Notons que l’oscillation globale u,* (43) n’a pas 

! (2 (02 
Le composantes 03% \sr—se OÙ WU ,sr—487 €b que 
l’amplitude des oscillations u,%* 4 est deux fois 


plus petite que l’amplitude définie par lPexpres- 


sion (44). 
Suivant la règle adoptée pour désigner les oscil- 
lations fondamentales, j’appellerai oscillation 


secondaire Sr + S’+’ l’ensemble des oscillations de 
fréquence |vs + vs] faites par tous les atomes du 
cristal. Cette définition admise, on trouve, pour 
SOU oscillations fondamentales, 


SIU(3IN — 1/2) = 3912 
oscillations secondaires. 
Les oscillations u,"s4s. (44) forment les pre- 


mières composantes — celles qui sont de degré 2 — 
des oscillations secondaires S+ + S’x’. 


cours du temps, suivant la loi exponentielle, comme la pro- 
babilité qu'a un atome de rester dans un état excité. Et, 
selon les relations d’incertitude (Heisenberg), tout quantum 
d'énergie hvs- cédé ou reçu par une oscillations fonda- 
mentale comporte, de ce fait, une indétermination A vs.. 
Cette indétermination affecte tous les échanges d’énergies 
par quantum hvs., en particulier, les échanges d’énergie 
entre les oscillations fondamentales. Aïnsi, supposer les 
fréquences vs- déterminées par le champ de forces existant 
dans le milieu cristallin est une hypothèse qui n’exclut pas 
l’indétermination hôvs. envisagée par les théoriciens qui 
traitent de la propagation de la chaleur dans les cris- 
taux [7-8]. 

En outre, attribuer aux oscillations fondamentales les 
énergies (7 ++ 1/2) kvs+, c’est supposer ces oscillations indé- 
pendantes les unes des autres. Les oscillations fonda- 
mentales, véritables (définies approximativement par les 
équations (39-58-...) sont liées entre elles par les 


m . . 
forces fs 5 ..… Etsi l’on recherche une expression 


ST 15 ST) 
plus précise de l'énergie thermique, par la méthode quan- 
tique des perturbations, on trouve des énergies de couplage 
entre les oscillations fondamentales qui s’ajoutent aux 
énergies propres (x + 1/2) hkvs.. En conséquence, même 
abstraction faite des relations d’incertitude, les échanges 
d’énergies entre les oscillations fondamentales ne se font pas 
exactement par quantum Avs+. 
(4) Les champs de force existant dans le milieu cristallin 
recèlent encore beaucoup d’inconnues ei, pour quelque 
champ particulier, la probabilité de Pégalité (46) peut ne 
pas être nulle. Supposons cette égalité satisfaite et les 
vecteurs d’onde $ conformes aux conditions cycliques donc 
égaux aux translations d’un réseau. En ce cas, la force élé- 
m 
mentaire, périodique, fe 7 &p se 
oscillations fondamentales $7 et $’+’ a la fréquence, 


Vers = |Vsr + vs, 
d’une autre oscillation fondamentale S$”+”, pilotée par le 
vecteur d’onde + (S.+ S’) ; elle est une composante de la 


force de rappel qui entretient l’oscillation ,$”+”. Il convient 


m 
donc de l’ajouter à la force fai de dans le second membre 
de l’équation (39). 


développée par les deux 
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B. Troisième opération. — [. Développée en 
série, l’expression de la force si (34) contient des 
termes qui comportent la fréquence vs. La 


force f35 . est définie par la somme de tous ces 
termes : 


m m ve 
j af j / RE 
Ja &, ST 30 &,ST F/3n &,ST ? 


el, compte tenu des expressions (23), (44), et des 
relations (31) nous avons : 


(is 
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FT PkB,sT—s/r" + Œiy,sr + S’(p — q) — spi]. (49) 


20 
m A$- 799 Vo 
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X ch 57 cos 2n{vsr t + Osr + onp,5r — spi 
et 
aie 5 
MTS Sr = (A5,) a se cs ,S’T’ 
X_C0S 2T[S"(4 — 7) — qiyss + pss,sr]. (50) 
Reportons ces expressions (48-49-50) des 


m 4 . [s rl: . 
forces f3 dans les équations (39) et éliminons le 
; ce ST . . . 
temps, c’est-à-dire : divisons les deux membres de 
chaque équation par le facteur 
exp.127( vs ? + 05 — Sm) 
Nous obtenons les équations séculaires (°) 
w? € = (I(S) + As(S)) (51) 
(5) A partir de ce point, je ne suis plus la méthode clas- 
sique de l'analyse mathématique. Je prends comme incon- 
nues dans les équations (39) les coordonnées vectorielles re 
el les amplitudes 4$+ (24) de Loutes les oscillations fonda- 
mentales, donc les coordonnées vectorielles et les amplitudes 


m 
qui figurent dans les expressions des forces 37 Au 


SC 
contraire, la méthode classique consiste en la résolution des 


m 1 
équations (39) où les forces fai sont déterminées par les 


ST 
coordonnées vectorielles &°$- et les amplitudes Ag+ (24), 
données par les équations séculaires (27) de la première 
analyse. Ainsi définies, les équations (39) admettent la 
solution : 

A$z Cja,sr X57 ! COS 2n{vsr t + Ogx + Pia,sr + Asr — Sm) 
+ Agr D 57 COS 27 (vsx t + Dia sr + 057 — SM) 
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28: est un coefficient déterminé, et 27xA$. une phase déter- 
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où la matrice AA($), tient compte des forces de 


rappel fi (a): 


ÉLTE ST mn 2 î le j le + 
ABUS) = D (As)? (Anb.s+sr + Anog,s sx 
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 k Fe 
+ Anañs.sr)] (52) 
Il vient des expressions (48), (49) et (50) 
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ce QU Ce EXP.— 127 S'ig—v)]] exp.127 S(m—p). 
(55) 


Au cours du temps, toute oscillation fonda- 
mentale, S+, perd ou gagne de l’énergie par quan- 
tum vs. Mais à température constante l’agitation 
thermique d’un cristal conserve la même énergie : 
les quantums d’énergie perdus par des oscillations 
fondamentales sont repris par d’autres oscillations 
fondamentales. Ainsi, les carrés des ampli- 
tudes A5, (24) fluctuent, mais, à chaque instant, 


minée, se rapportant l’un et l’autre à l’oscillation fonda- 
mentale ST. Au contraire, pour déterminer les coef- 
ficients b%" s- et les phases 276/ju.8 (j — 1,2, ...,g, v — 1, 
2, 3), il est nécessaire de se donner l’un de ces coefficients et 
l’une de ces phases. L'introduction des « oscillations anor- 
males », définies par l’expression précédente, entraîne des 
calculs inextricables. Mais le premier terme de cette expres- 
sion, celui qui représente une oscillation dont l’amplitude 
croit avec le temps, révèle que les fréquences déterminées 
par les équations séculaires (27) de la première analyse ne 
sont pas les fréquences véritables des oscillations fonda- 
mentales. : 
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leurs fluctuations se compensent mutuellement. 
Elles n’entraînent pas une variation sensible des 
composantes tensorielles A3 di(S) (2). Et, pour 
évaluer ces composantes, il convient de prendre en 
compte dans l'expression (52) les amplitudes qua- 
dratiques moyennes, As. 

J'ai fait état d’oscillations fondamentales qui 
forment un seul train d’onde emplissant tout le 


cristal ; de la sorte, les oscillations w,5 (23) faites 
par tous les atomes ont, au même instant, la même 
phase temporelle 2r04.. En fait, les trains d’ondes 
élastiques sont de volume borné et leurs dimensions 
diminuent quand la température s’élève [7-8]. Les 
différences entre les phases temporelles, 2r04., des 


. . m ; A . 
oscillations w,i,. transportées par des trains 


d'ondes distincts et la perturbation de ces oscil- 
lations aux frontières qui séparent les trains d’onde, 


influent sur l'intensité des forces de rappel fs 5 tte 


49-50). Ces forces ne sont pas strictement définies 
par les expressions (48-49-50). Maïs les trains 
d'ondes formés par une même oscillation fonda- 
mentale, S+, conservent à toute température des 
dimensions très grandes par rapport à la longueur 
d’onde élémentaire À, celle qui se rapporte au 
vecteur d'onde fondamentale $S : À — 1/|S]. Entre 


m PARC m 
les forces fu véritables, et les forces fs 


définies par les expressions (48-49-50), l'écart reste 
toujours minime. D’autre part, à température 
constante, les trains d’ondes élastique gardent, 
dans l’ensemble, les mêmes dimensions. De même 
que les échanges d’énergie entre les oscillations 


Q 2 6 . , m 
fondamentales, l’incohérence des oscillations 4 RER 


transportées par des trains d’ondes différents, ne 
provoque pas au cours du temps, une fluctuation 
sensible des composantes tensorielles À, Z(S). A 
température constante, ces composantes peuvent 
être tenues pour constantes (°). 

La matrice A,(S) (52-53, 54, 55) est hermitienne, 
etA3(— S) — Af(S). D'autre part, les composantes 
tensorielles A, Z£(S), comparées aux composantes 
tensorielles TË(S), ne sont qu’une correction. Les 
équations séculaires (51) peuvent être résolues par 
itération. On porte dans la matrice A,(S), au lieu 


: RG <. - 
des inconnues €, et Ag., les coordonnées vecto- 


rielles s. et les amplitudes quadratiques 


moyennes A3; (7) déterminées par les équations 


(*) Notons encore que les composantes tensorielles 


A3 «8 (18) ne varient pas avec le volume du cristal, sauf si ce 
volume devenait extrêmement petit. Pour un cristal 
de T atomes, la somme (52) compte 39T termes ; mais les 
carrés des amplitudes quadratiques moyennes As% étant en 
raison inverse de la masse du cristal (24), chaque terme est 


inversement proportionnel à 9T; donc la somme reste 
constante quand varie. 


(7) L’amplitude quadratique moyenne Ag, est déter- 
minée par la pulsation og; (56). 
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séculaires (27) de la première analyse. Ainsi 
modifiées, les équations séculaires (51) peuvent être 
résolues ; elles donnent une première solution 
approximative. Puis on inscrit les coordonnées vec- 


torielles ei s- et les amplitudes quadratiques 


moyennes À,4. fournies par cette première solution 
dans la matrice A,(S), et l’on résout de nouveau les 
équations obtenues par cette seconde modification 
des équations séculaires (51). Cela donne une deu- 
xième solution approximative moins erronée que la 
première, etc. On poursuit le calcul par ce procédé 


jusqu’à ce que les coordonnées vectorielles CP s. et 


les amplitudes quadratiques moyennes À,4.calculées 
se confondent sensiblement avec les coordonnées 
vectorielles et les amplitudes quadratiques 
moyennes portées dans la matrice A(S). 

Les oscillations fondamentales S+ déterminées 
par les valeurs caractéristiques et les vecteurs 
propres de la matrice de Fourier (27-28) sont indé- 
pendantes les unes des autres. Chacune a l’énergie 
d’un oscillateur harmonique et l’amplitude quadra- 
tique moyenne 


Hiu2a 1 tie) 
#) 


M ñ os 
L ST exp ST 


NI 


Les oscillations fondamentales, déterminées par 
les équations séculaires (51) sont liées entre elles 


par les forces fai (ss (48). Elles peuvent être encore 


quantifiées (approximativement) par la méthode 
des perturbations. Mais le calcul des amplitudes 


quadratiques moyennes À,, recèle ensuite des diffi- 
cultés. Présentement il est impossible de le mener à 


bien faute de connaître les grandeurs véritales des 
dérivées troisièmes, Dÿki, et quatrièmes, Eïkie, de 
l’énergie potentielle [1-2] (5). Toutefois, les théories 
classiques (Debye, Max Born), qui font seulement 
état des oscillations fondamentales, et les sup- 
posent indépendantes, retombent sensiblement, 
sauf aux hautes températures, sur les chaleurs spé- 


ASC Aosz 


cifiques mesurées. L'écart relatif 
ST 
reste sans doute minime et peut être négligé en 
première étude. 


20 Retranchons de la force P32 (34) toutes les 
forces élémentaires périodiques Étèpe qui ont les 


1 : . . M 
fréquences des oscillations fondamentales. Soit 9 
la force restante. 


(#) Ce calcul porte sur des séries. Et pour définir, même 
A$- ES 


avec une précision modeste l'écart | ==" |, il serait 


nécessaire d'effectuer de longs calculs numériques. 
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Les équations 
M M—D 6 °m 
tj Us? E 2 C5 ê UÊY? = Pa? 


sont décomposables, comme les équations (42), en 
systèmes de 3g équations de d’Alembert, qui se 
rapportent à des forces élémentaires far 

" : 3@ ,STLS/T' +87" 
toutes de la même fréquence [vg + ve + vor. 
De sorte que chaque système de 3g équations déter- 
mine une oscillation secondaire S7+ + Sr + S’+”, 
d’élongation 


As Age Agrer Ci srEs rs" e" 
X cos 27 [(Vsr + Vs + VS”x") t 
+ Vs Æ Ogre Æ Ogre + piagrtsrt sr" — (S + S’ + S”) m]. 
(57) 


Les coefficients Cérsrss et les phases 
TO sr18r25 SOnt donnés, comme les coefficients 
grise (44) et les phases, 27@s15 (44), par la 
solution des 3g équations séculaires que l’on obtient 
en éliminant le temps des équations de d’Alembert. 

Les coefficients ©” ses. sont du second degré 

F0 » +. MDI 
par rapport aux dérivées troisièmes D dt et du 
premier degré par rapport aux dérivées qua- 
trièmes E his de l’énergie potentielle [1-2]. 

Ainsi, la seconde analyse décompose l’agitation 
thermique en oscillations fondamentales, S+, déter- 
minées par les équations séculaires (51) et en oscil- 
lations secondaires Sr + S’+' (44) de degré 2, et 
Sr + S'r + S’r" (57) de degré 3, respectivement 
de fréquence |vs, + ve et [ve ÆE Ve Æ ver. 


VI. Troisième analyse. 


La troisième analyse (équations (38)) fait état 
d’oscillations fondamentales déterminées par les 
équations 


D pô m mm 
ë L . — 1 1 
f Uk ,sT Îs œ,ST $ 5 GrST * (58) 


on 


Uj U, is = a « 


Elle comprend les deux opérations de la seconde 
analyse et se prolonge par deux opérations nouvelles 
qui peuvent être effectuées selon la méthode suivie 


en seconde analyse. 
Je donnerai seulement les résultats. 


Quatrième opération. — Classées selon leur 
degré, les forces élémentaires comprises dans la 
force ai se rangent en 3 catégories. Ce sont : 

49 Des forces sensiblement constantes au cours 
du temps, développées par les oscillations fonda- 
mentales +, et par les premières composantes des 


2 . . al <a 1 DB À Vif * EE 
oscillations secondaires ST + ST, celles qui sont 
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de degré 2 ; 20 des forces périodiques dont la fré- 
quence est [vs + vsl. Les oscillations qu’elles 
entretiennent, de degré 4, forment les secondes 
composantes des oscillations secondaires 87 + Dé 
30 ce sont enfin des forces périodiques de fré- 
quences [ve + ... + v»»| qui entretiennent des 
oscillations ayant la même fréquence, pilotées par 
les vecteurs + (S + S’ + S" + S”). 


Cinquième opération. — 1° En éliminant le temps 
des équations de d’Alembert (58) on trouve les 
équations séculaires 


t = (T(S) + A, (S) + As) € 


(59) 


où la matrice A,(S) tient compte des forces 1: VX 
Ces équations peuvent aussi être résolues par ité- 
ration. On porte dans les matrices A,(S) et A;(S) 
les coordonnées vectorielles aps et les amplitudes 
quadratiques moyennes A4, déterminées par les 
équations séculaires (51) de la seconde analyse. En 
résolvant les équations (59) modifiées de la sorte, 
on obtient une solution approximative qui peut être 
prise comme point de départ-du calcul itératif. 


29 Retranchons de la force P5i toutes ses com- 
Dosantess see (rte ON tt NE 
3g). La force restante est la résultante des forces 


élémentaires qui entretiennent des oscillations de 
fréquences 


Mae Var l: 


LVsr Æ Vos ÆE vevl et [ve + 
Les oscillations de fréquence [vs + vs + vire 
forment les secondes composantes des oscillations 
secondaires ST + S’x’ + Sr”, 

Aiïnsi, de proche en proche, l'agitation thermique 
peut être résolue en composantes harmoniques. Et 
le nombre de ces composantes et la précision de 
leur définition croissent avec le rang de l’analyse 
effectuée. 


VIT. Conclusion. 


Les forces de rappel développées par l'agitation 
thermique n’obéissent pas à la loi de Hooke. Cela 
entraine les conséquences suivantes : 

19 Mus par l'agitation thermique, les atomes 
exercent les uns sur les autres des forces répulsives 
sensiblement constantes dans le temps, qui crois- 
sent avec la température. Ces forces constituent les 
tensions thermiques ; elles dilatent le milieu eris- 
tallin. 

Si l’agitation thermique développait des forces 
de rappel conformes à la loi de Hooke, elle ne pro- 
duirait pas de tensions thermiques, done pas de 


(*) En toute rigueur il serait nécessaire de redéfinir les 
amplitudes quadratiques moyennes, Agz,à la précision 
atteinte en troisième analyse. 
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dilatation. Or, une élévation de la température pro- 
voque toujours une dilatation du milieu cristallin 
libre de toute contrainte, même si la température 
initiale est très proche du zéro absolu. Donc, même 
aux basses températures, l’agitation thermique 
développe des forces de rappel qui ne restent pas 
proportionnelles aux élongations. 

20 L’agitation thermique n’est pas analysable en 
oscillations harmoniques indépendantes, au nombre 
de 39T dans un cristal formé de 9T atomes. ‘Elle 
comprend des oscillations fondamentales, au 
nombre de 39, et des oscillations secondaires qui 
sont, au sens large, les harmoniques des oscillations 
fondamentales. Les oscillations secondaires sont 
entretenues par des forces de rappel qui sont du 
second, du troisième degré, etc. par rapport aux 
amplitudes des oscillations fondamentales. Aux 
basses températures, voisines du zéro absolu, les 
oscillations secondaires sont infimes ; mais aux 
hautes températures, elles cessent d’être négli- 
geables. À proximité du point de fusion, leur ampli- 
tude quadratique moyenne, totale, excède proba- 


blement chez de nombreux cristaux, le dixième de- 


l'amplitude quadratique moyenne atteinte par 
l'agitation thermique globale. C’est laraison princi- 
pale de l’écart entre les chaleurs spécifiques prévues 
par les théories classiques et les chaleurs spécifiques 
mesurées sur des cristaux portés aux températures 
élevées. 

Comme les oscillations fondamentales, les oscil- 


lations secondaires produisent dans le milieu cris- 


tallin des trains d’ondes de densité électronique et 
des trains d’ondes de potentiel électrique. Elles con- 
tribuent de la sorte à la diffusion des radiations 
électromagnétiques et des électrons. En particulier, 
elles renforcent la diffusion des rayons X [10], et 
augmentent la résistivité des métaux. 

30 La fréquence d’une oscillation fondamentale 
est fonction des puissances paires des amplitudes 
quadratiques moyennes de toutes les oscillations 
fondamentales. 

À élongations égales, la force de rappel véritable 
exercée sur un atome est en général moindre que la 


M . 
force de rappel f, j (20) conforme à la loi de Hooke ; 
elle est semblable en cela au couple de rappel 
exercé sur le pendule gravifique. De même, la force 
m . . . È 
de rappel f5 s- (39), qui entretient une oscillation 
fondamentale est, en général, plus faible que la 
£ m . 
force de rappel, hi s- (26), proportionnelle à l’élon- 
gation. En conséquence, les fréquences vs, des oscil- 
lations fondamentales sont, dans leur ensemble, 
moins élevées que les fréquences vs. déterminées 
par les équations séculaires (27) prises en compte 
par les théories classiques. Et, comme la diffé- 
rence [fi 5 ge "15 seb l'écart [vs — vos.l croît avec 


Pamplitude des oscillations fondamentales, donc 
avec la température. 


N°04 | 


1 


milieu cristallin et toutes ses dérivées CE 

mpq 
D jkl,...,(6 bis) restent constantes quand la tem- 
pérature s'élève, ce qui a sensiblement lieu dans un 
cristal maintenu à volume constant, si les réseaux 
formés par les atomes en même position dans la 
maille ne se déplacent par les uns par rapport aux 
autres et si les transitions des électrons faiblement 
liés vers des états d’énergie plus élevée sont en 
nombre réduit. Cela étant, lorsque la température 
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Supposons que l'énergie potentielle du s’élève, les oscillations fondamentales s’amplifient, 


donc leurs fréquences diminuent. Mais, pour une 
même élévation de la température, cet abaissement 
de fréquence, dû exclusivement au fait que les 
forces de rappel n’obéissent pas à la loi de Hooke, 
reste certainement petit par rapport à l’abaissement 
de fréquence, considérable, qui est entraîné par la 
dilatation thermique du milieu cristallin libre de 
toute contrainte. 


Manuscrit reçu le 14 octobre 1958. 
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REVUE DES LIVRES 


OrTust (J.), Propagation des ondes électromagnétiques de 
haute fréquence (1 vol. relié, vi + 320 pages, 17 X 25 cm, 
Société Française de Documentation Électronique, 
Paris, 4957, 3 100 F). 

Il s’agit du premier ouvrage de la « Collection des Annales 
de Radioélectricité », dirigée par M. Maurice Ponte. Le but 
de cette collection est de faire le point des chapitres les 
plus évolués de l’électronique et de ses applications, en dé- 
gageant les traits essentiels de la technique exposée, sans 
entrer dans les détails. A ce titre le livre de M. Ortusi ne 
peut manquer d’intéresser bien des physiciens, ingénieurs 
ou chercheurs. On notera cependant que les présccupations 


pratiques et utilitaires ne sont pas absentes. Comme l’auteur 


nous le dit lui-même : « Cet ouvrage est rédigé de façon 
à servir de cadre aux ingénieurs chargés de résoudre cer- 
tains problèmes de canalisation de l'énergie pour obtenir 
une utilisation plus rationnelle de ses propriétés. » 

La table des matières indique cinq parties : 1) L'énergie 
électromagnétique ; 2) Les ondes guidées progressives ; 
3) La propagation dans les milieux anisotropes ; 4) Les 
branchements de guides, et 5) La propagation autour de 
la surface terrestre. On remarquera entre autres le chapitre 
consacré à la propagation dans les guides renfermant des 
ferrites (Chap. III de la 3° partie) et le chapitre concernant 
certains branchements de guides : filtres de fréquence, 
gyrateurs, isolateurs, jonctions hybrides et circulateurs 
(Chap. IV de la 4° partie). 

Malgré l’étendue du sujet, l’auteur ne s’est pas contenté 
d’une simple compilation des résultats éparpillés dans la 
littérature technique. Mais il a édifié une théorie simple des 
transformations de lénergie électromagnétique, ce qui 


donne de l’unité à son ouvrage tout en lui fournissant bien 
des aperçus très originaux. 
M. JEsser. 


Hertz (G.), Traité de physique nucléaire I. Méthodes 
expérimentales. (Lehrbuch der Kernphysik, I. Experi- 
mentelle Verfahren) (en allemand). Un vol. de 227 p., 
relié, 16,5 x 23 cm, Teubner, Leipzig, 1958, prix 
DM°#15: 

Il s’agit d’un ouvrage collectif, chaque chapitre est d’un 
auteur différent, et qui s’adresse aux expérimentateurs et 
aux ingénieurs. Les connaissances mathématiques requises 
sont simples. Deux autres volumes doivent suivre celui-ci : 
physique nucléaire proprement dite, et applications. 

Introduction : Généralités sur les photons, les électrons, 
les atomes et les noyaux. À) I. Détermination de la charge 
nucléaire. II. Détermination du rayon nucléaire. III. Masses 
nucléaires. IV. Moments nucléaires. 

B) Preuves et détermination des rayonnements de 
haute énergie. V. Courants d’ionisation, mesure. VI. Par- 
ticules isolées : rayons gammas, interaction avec la matière, 
observation, comptage. Appareils électroniques. Spec- 
tromètres pour rayons alpha, bêta, gamma. 

C) Accélérateurs de particules. VIT. Accélérateurs li- 
néaires. VIII. Accélérateurs circulaires (bétatron, cyclotron, . 
synchrotron pour protons et électrons, et avec champs 
alternés. IX. Développement probable. 

Cet ouvrage contient 139 figures et photos d’une grande 


clarté. 
J. WINTER. 


LE JOURNAL DE PHYSIQUE ET LE RADIUM 


à 


L’ANALYSE DES MÉTHODES DE MESURE DU COEFFICIENT y 
DES GRANDES GERBES DU RAYONNEMENT COSMIQUE 


Par R. FIRKOWSKI et A. ZAWADZKI, 


Institut des Recherches Nucléaires de l’Académie Polonaise des Sciences. 


Résumé. — L’analyse formelle des deux méthodes de mesure de l’exposant Y (A : variation du 
rang de coïncidences, B : variation de la surface des bancs de compteurs) montre qu’elles ne sont 
pas complémentaires indépendamment de l’existence du gradient de densité dans la section hori- 
zontale d’une grande gerbe. La présence du gradient de densité entraîne un accroissement apparent 
de la valeur de y déterminée par la méthode A. : se 

En utilisant un dispositif hodoscopique qui permet l'enregistrement des coïncidences triples, 
quadruples, quintuples et sextuples en même temps, on démontre que la valeur de y obtenue par 
la méthode A est supérieure à la valeur donnée par la méthode B, et que cette valeur croît avec 
l’ordre de coïncidence. Toutes les conclusions ont été tirées d’une série unique de mesures compor- 
tant plus de 124 000 enregistrements des coïncidences. Par une méthode appropriée de calcul on a 
mis à profit toutes les combinaisons d’enregistrements, grâce à la possibilité de passer du nombre 
des coïncidences de l’ordre déterminé au nombre des coïncidences qui sont au moins de l’ordre 
déterminé (p. ex. des coïncidences triples aux coïncidences au moins triples). Cette circonstance 
diminue d’une façon appréciable la valeur des erreurs statistiques. 


Abstract. — The formal analysis of the two measurement methods of the y exponent (A : varia- 
tion of the order of the coincidence, B : variation of the area of the counter trays) shows that they 
are not complementary independently of the existence of the density gradient in the horizontal 
section of the extensive air shower. The presence of the density gradient causes an apparent 
increase of the y value as determined by the method A. By using a hodoscopic arrangement witch 
allows one to record simultaneously threefold, fourfold, fivefold and sixfold coincidences we can 
show that the Y value obtained by method A is higher than the value obtained by method B 
and that the value increases with the order of coincidence. All the conclusions have been based 
on a single set of measurements of 124 000 coincidence records. By using an appropriate comput- 
ing method we have taken advantage of all possible recording combinations ; this was possible 
due to the transition from coincidences of a given order to coincidences which are at least of the 
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given order (e.g. from threefold coincidences to at least threefold coincidences). 


derably reduces the statistical errors. 


1. Introduetion. — On peut représenter le spectre 
des densités des grandes gerbes sous forme de 
Pexpression différentielle suivante : 


N{e) de = Kp= (rt dp (1) 


où W est le nombre des gerbes de densité comprise 
entre p et p + de. 

La fréquence d’enregistrement des gerbes par un 
système de coïncidences de » bancs identiques de 
compteurs est égale à : 


I = KSt Le CU e—a)e av 40 de @) 
où æ — PS est le produit de la densité p des parti- 


cules par la surface S de chaque banc de compteurs. 
On obtient pour la valeur J, : 


= KT D 6 1H (7) er. (3) 


En se servant de la valeur J, on peut déterminer. 


la valeur y ; à cet effet il faut effectuer au moins 
deux mesures en vue d'éliminer la constante K. 

Le but du présent travail est l’analyse de deux 
méthodes de mesure de l’exposant y. Les deux 


This fact consi- 


méthodes s'appellent dans la nomenclature de 
Cocconi (1949) A et B. 

À. La méthode de variation du rang des coïnci- 
dences qui consiste à comparer la fréquence des 
coïncidences de n°% ordre avec la fréquence des 
coïncidences de (n — 1) ordre : 


Jai) 2 (1) (or 


TS) a D pe 1) ( na i “ 


(4) 


B. La méthode de variation de la surface des 
compteurs qui consiste à comparer les fréquences 
des coïncidences du même ordre entre les bancs de 
compteurs dont chacun a changé la surface dans la 
même proportion o : 


In Jn(oS) 
Jalos) , JS) 
SAT) = G QUE ES (5) 


La méthode A exige la supposition d’une densité 
constante moyenne dans toute l’envergure des 
compteurs et par là même néglige l’existence du 
gradient de la densité dans une section horizontale 
d’une gerbe. 


N° % 


L’emploi de la méthode B n’exige cependant pas 
la supposition de la constance de la densité parce 
que la fonction de la répartition des densités des 
particules dans une gerbe [o(r), où r est la distance 
de l’axe] est à peu près la même pour toutes les 
grandes gerbes. 

On peut démontrer que les valeurs de y obtenues 
au moyen de la méthode A sont pour cette raison 
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faussement augmentées parce que le gradient de la 


densité dans la section horizontale d’une grande 


gerbe cause une différence mesurable des densités 
entre deux bancs de compteurs éloignés de quelques 
mètres l’un de l’autre (Broadbent, 1950). 

Nous présentons ci-dessous une comparaison de 
quelques résultats obtenus par différents auteurs 
sur le niveau de la mer (approximativement). 


Ne a MÉTHODE À Mérxone B 
Auger, Daudin (1945) rm 70 1,66 
Cocconi et coll. (1946) 10 — 1 000 1,46 
Latsepine et coll. (1947) DOS) 1,43 
Maze et coll. (1948) > 100 1,67 
Cocconi, Cocconi-Tongiorgi (1949) 5 — 300 1,42 17976) 
Broadbent et coll. (1950) 5 500 4,425 
Hodson (1953) 5 — 500 1,58 1,45 
Murdoch (1958) 1 1,34 


(*) Calculé par nous comme moyenne des résultats obtenus séparément pour les coïncidences triples ‘et quadruples. 


On voit qu’il existe des différences entre les résul- 


tats obtenus par la méthode A et B, mais cepen- 
dant les divergences entre les données des différents 
auteurs rendent impossible une conclusion défi- 
nitive au sujet de l’applicabilité des deux méthodes. 
Ces divergences découlent de l’utilisation des diffé- 
rents types de dispositifs expérimentaux, de diffé- 
rents degrés d’exactitude et de divers autres fac- 
teurs qu’on ne peut pas éliminer. Il faut pourtant 
remarquer que la question de la différence entre les 
deux méthodes n’est pas définitivement éclaircie 
dans la littérature. Dans ses travaux Daudin (1948 
et 1949) a essayé de justifier d’une façon théorique 
et expérimentale l'identité des deux méthodes. 
Dans son rapport monographique Deutschmann 
(1953) n’attache pas d'importance à la différence 
mentionnée. La même chose concerne aussi les 
travaux de Zatsepine (1950) et de Dobrotine et 
coll. (4953). Dans son étude approfondie au sujet 
des grandes gerbes Greisen (1956) suggère qu’on 
peut éviter la différence entre les méthodes À et B 
en utilisant un dispositif expérimental d’une géo- 
métrie appropriée comme celui utilisé dans le 
travail de Coecconi (1949). À notre avis (comme 
nous essayerons de le démontrer dans la suite) cette 
différence peut être diminuée mais ne peut pas être 
complètement éliminée. | 

Le but du présent travail est l’analyse de deux 
méthodes À et B en vertu d’une série unique de 
mesures, ce qui permet d'éviter des facteurs pro- 
venant de la différence qui est liée toujours d’une 
façon plus ou moins expressive avec deux séries de 
mesures indépendantes. 


2. L’analyse du rapport formel des méthodes A 
et B. — Daudin, dans ses travaux (1948, 1949), 
conclut que les deux méthodes ne doivent pas 
fournir de résultats différents parce qu’elles sont 
complémentaires. La preuve fournie par Daudin est 
formelle et ne doit pas dépendre de la présence ou 
de l’absence du gradient de densité. Son raison- 
nement est le suivant : Soit trois bancs de comp- 
teurs des surfaces égales S dont un est séparé des 
deux autres (qui sont juxtaposés) par une distance 
de quelques mètres. Les bancs des compteurs 
peuvent être mis en coïncidences doubles (S, S) et 
(S + 8,85) — (28, 8) ou triples (S, 5, S). La coïnci- 
dence (45,5) (des deux bancs séparés évidemment) 
se compose de la coïncidence (9, S, S) et de la 
coïncidence (S, S$ non S) ce que nous écrivons en 
abrégé : 

(S, 8) = (S, 5, $) + (5, 5,5) ; 


d'autre part : 
SAS) ES SE CESSE) 


En éliminant le terme contenant la négation $ 
on obtient : 


DÉS) = 8 25h (SAS). 


Ce que l’auteur appelle une dépendance complé- 
mentaire, c’est qu’en connaissant les valeurs (S, S) 
et (IS, #, S) utilisées dans la méthode A, on connaît 
aussi la valeur ($, 25). Si la valeur (S, 2.5) pouvait 
vraiment, comme le juge l’auteur, servir dans les 
calculs de la méthode B, sa preuve formelle (dans 
laquelle on ne fait aucune supposition au sujet du 
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gradient de densité) suffirait à démontrer l'identité 
des deux méthodes. La valeur (S, 285) ne peut 
pourtant pas s'appliquer à la méthode B, parce 
qu’elle ne permet pas d'utiliser la formule (5). [Les 
deux intégrales (2) s’éliminent de l'expression (5) à 
condition de changer de toutes les surfaces #S° par 
le même facteur o.] 

On peut poser la question si cette insuffisance ne 
s'ensuit pas du caractère incomplet du dispositif 
admis. Le dispositif (S, 5, S) ne peut pas, notam- 
ment, servir aux mesures par la méthode B. Pour 
les mesures de ce type il faut avoir au moins deux 
banes séparés de compteurs : (5, 5) et (5, 5). Nous 
avons étendu le raisonnement de Daudin à tel sys- 
tème de compteurs, qui peut être mis soit en coïnei- 
dences doubles ($, $) (1), soit aussi en doubles 
(S + S,S + S) = (25,25), soit en triples (5,5, 5) 
soit enfin en quadruples (IS, 15,5, 5). 

On voit facilement que : 


LS}, (ST =, 15, 8,18) + (5,25, 15,18) 

CSS) CSS) I 
EMA AS RE RICA ET PE) 

LaTCS, 5), CS, 8)] TI 
RSS) US RS 40019) LES ERS 20 00) III 


On élimine les négations doubles des égalités I 
et IT : 
AS), (5) — (2, 25) = 3(S2, 15, S) + 4ÛS, 5, 5,15): IV 


On élimine les négations simples des égalités TIT 
et IV : 
4 [US'), (5) Er à (25, 25) = n (US, 5’, S) ER (Se 15% 5’, S). 


En mesurant [(S), (S)]et (25, 28) c’est-à-dire en 
se servant de la méthode B on connaît la différence 
des valeurs qui apparaissent dans la méthode A 
(triples ou quadruples), mais on ne connaît pas les 
deux valeurs séparément, ce qui est nécessaire pour 
appliquer la méthode A. Dans ce but il fallait effec- 
tuer encore une mesure : (9,5, 5) ou (8, 45, 8, 8). 
On voit alors que dans ce cas également les deux 
méthodes ne sont pas complémentaires. 

Nous avons poursuivi le même raisonnement en 
lappliquant au système des compteurs utilisés dans 
le présent travail et notamment : trois groupes 
séparés l’un de l’autre (IS, 5), (S, S), (S, S) et 
coïncidences au moins triples actionnées par au 
moins une impulsion provenant de chaque groupe. 
On a alors : 


DNS AS ST Sen CES te IR SES ES) 
SSSR LOIS) PRES IE AE AR RCE AG I 


(2) Il est évident que dans ce cas les deux ,$ appar- 
tiennent aux différents bancs, [on peut éventuellement, 
pour souligner cette évidence, utiliser la notation sui- 
vante : [(S), (S)]. 
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No [A 
(25, 25, 25) = (S, 5, 15,.5, 8, 5) + 6(S, 5, 5, 5,8, 5) 
Le ASS Sr SSI EE CSSS ASE 


IT 
(5, 5”, ISA 1} = US, #, 5: 1, Ne 5) 
n9(5, SSL SMS) CSS MINES ESRI 
LS SA SS, S) HS US SN SNS) ICS SRE 


IV 


On élimine les négations triples des égalités I et 
IT 


BUS 5 8) = (22 SSL ESS ETES 


LABS SSL S) MOSS SSSR DES 


On élimine les négations doubles des égalités III 
et V 
42(5,.5, 8, 5) —8(5,15, S) + 4(28,26,28) 


= 58,8, 85,188) 6(S, SAS AS ESS 
VI 


On élimine les négations simples des égalités IV 
et VI 


COS EE DA OMC D CC 


ps à G(S, 0 5; 195 S) +- (Ss, 15 1”, 15 S’, 5). 


Ce qui montre aussi que les deux méthodes ne 
sont pas complémentaires. En utilisant la 
méthode B il faut en plus connaître deux des coinci- 
dences : (S, 5, 5, S), (S, 5, 5, De S), (5, S, S, 5, 5,5), 
pour pouvoir appliquer la méthode A. 

On peut conclure d’une façon générale que, indé- 
pendamment de la question du gradient de densité, 
les méthodes A et B fournissent des résultats qui ne 
sont pas complémentaires. 


3. Le rôle du gradient de densité. — L'existence 
du gradient de densité dans la section horizontale 
d’une grande gerbe introduit une erreur systéma- 
tique dans la détermination de Pexposant + par la 
méthode A. Cet effet était déjà mentionné par 
Broadbent et coll. (1950) ; regardons le de plus 
près. Soit deux bancs juxtaposés de compteurs de 
surfaces identiques #. Une gerbe de densité p parti- 
cules par m? (ou x particules par un banc) actionne 
les deux bancs avec une probabilité (1 — e=*)2. Si 
l’on éloigne un banc de autre et si Pon admet que 
la densité x reste la même sur un banc, la densité 
sur l’autre banc est x + x ou x — 5x. La proba- 
bilité d’actionnement de ce banc est dans le premier 
cas [1 — e®+8] dans le deuxième cas 


___ p—{&—0x) 

FLE ET 

en moyenne alors À — e?.(e7 + ex) J2 pourtant 
3 


Ga? (Bx)t 
SE NN 


| 


(et et 


« 


[Ne 
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par conséquent la probabilité d'enregistrement de 
la gerbe par le deuxième banc est moindre si l’on 
admet la présence du gradient de densité (2). 

Au lieu de deux bancs de compteurs on peut, 
de la même façon examiner n bancs, on voit alors 


que lerapport J,(S)/J,-,(S) est toujours diminué 


par la présence du gradient de densité et que la 
valeur déterminée de y est en apparence augmentée. 
Il est évident que la différence de densité sur les 
différents bancs ne dépend pas directement de la 
distance entre les bancs mais de la distance comptée 
radialement à partir de l’axe de la gerbe ; autre- 
ment dit c’est la différence des distances de l’axe 
qui décide. Il est justifié d'admettre que durant 
une longue série d’enregistrement les axes des 
gerbes tombent d’une façon isotrope autour du 
dispositif expérimental. Dans ce cas il est facile de 
démontrer que les « distances radiales » entre les 
différents bancs de compteurs sont égales aux dis- 
tances réelles multipliées par le facteur 2/x. 


4. Dispositif expérimental. — L’hodoscope dont 
on se servait dans le présent travail, était composé 
de six bancs de compteurs GM, chacun comportant 
neuf compteurs de surfaces s identiques ; la surface 
effective s — 130 cm?, la surface d’un banc S — 95, 
= 0,117 m?. Chaque compteur possédait son 
amplificateur individuel. 

Les bancs de compteurs étaient disposés par deux 
dans les sommets d’un triangle rectangulaire dont 
les côtés étaient égaux respectivement à 5,5 m; 
6,5 m ; 8,5 m. Les compteurs dans chaque banc 
étaient écartés de 3 cm l’un de l’autre et placés 
directement sous un toit d’environ 0,9 g/cm? 
d'épaisseur (0,6 cm de contre-plaqué et 0,3 cm de 
carton bitumé pour toiture). 

L'enregistrement des coïncidences triples exi- 
geait qu’au moins un compteur dans chaque paire 
de banes fût actionné. Le temps de résolution du 
système central des coïncidences était de 0,5 y sec, 
le temps de résolution des coïncidences dans chacun 
des bancs de compteurs était de 2 a sec. ; 

On enregistrait alors au moyen d’un appareil 
cinématographique les allumages des lampes à 
néon qui signalaient les coïncidences triples Ses 
S, $, $, S), quadruples (S, S, S, S, $, S), quin- 
tuples (S,8, 8,9, 8,8) et sextuples (5, 5,8, 8,8, 8) 
dans toutes les combinaisons entre les bancs parti- 
culiers. 


(2) En négligeant les termes d’ordre supérieur de 2 on a 
pour la probabilité d'enregistrement de coïncidence double 
par les gerbes d’une densité donnée : 


L o id 
Ro [4 — eT €) TS e—æ(ôx:) | (1 — 7%) 
et pour le rapport de cette probabilité à la probabilité (P)0 
d'enregistrement s’il n’y avait pas de gradient de densité : 
] . —— 
Pa 1 À (Bx)?. 
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5. Résultats des mesures. — On calcule le 
nombre des coïncidences, au moins triples, dans 
toutes les combinaisons entre les bancs particuliers : 
&(S, 1, S) —= (S, W, 9, Se Se $) ie AS 15), f’, #, S) 

SEA CCS CCC ANS RE PR FE ESS 
le nombre des coïncidences au moins quadruples : 
AIO OA OI A NRA 9 EPA OT CO SR 
HAS S STONE 
le nombre des coïncidences au moins quintuples : 
GLS SES 8 5) CS ESS ES SEE GS SSSR S)e 
le nombre des coïcidences au moinstriples entre les 
paires de bancs juxtaposés : 
(CSA ENCRES ER ER ER EE SC OS) 
= (65 ’, 1’, S" st $) + (S 1’, 155, 5’, ou $) 

(ST SNS) OS ES RSS ES 

On obtient ainsi les valeurs : (S,8, 8), (S, 8,8, S) 
(S,15515, 19,0) (SD ASS Set (25 285225), 


qui servent pour les calculs dans la méthode A en 
utilisant par exemple les rapports : 


APR ST) (SES a Ie pe CSSS) 
RS NS TO PO GES TO GO) 
et dans la méthode B : 
(2,5, 29. 25) 
(S, 8,5) 


En profitant de l’enregistrement hodoscopique au 
lieu du simple système de coïncidence, on diminue 
les erreurs statistiques d’une série de mesures. 
Examinons les erreurs commises dans Ja 
méthode B : le nombre (28, 2$, 2S) est dépourvu 
d’une erreur statistique quelconque, puisque c’est 
le nombre total des enregistrements. Pourtant au 
nombre ($, 8,5) on attribue une erreur calculée en 
partant de (25, 28, 25) parce que c’est le nombre 
comprenant toutes les combinaisons indépen- 
dantes (S, S, $). Autrement dit l'erreur relative 
commise sur (S, S, S) est égale à l’erreur relative 
statistique du nombre total des coïncidences (25, 
28, 28). De la même façon on calcule dans la 
méthode A lerreur sur (S, 8, S, $) en partant de 
la valeur [ (25, 28, 25) — (5,85, 5, 5,8, S)] l'erreur 
sur (S, 9,48, 9,8) en partant de es 
OS I ns ES AS SERGE ST S SMS 
et l'erreur sur (S, $, S, $, $, $) directement de 
cette valeur puisque il existe une seule combinaison 
de ce genre : [(25, 25, 25) — (S, S, S, 8, S, S) 
A D D LS 0 50 Do) = (020, 
SAONE 

C’est ainsi que Perreur sur lexposant 
(HSHAS STE 
RICO IDETTE) 
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calculée par la méthode B est : 


AY A(2S, 28, 2) El 1 
Ter (AS MSALS) Vas ES) 
0] HOMO RS RE D , , 
(28, 25, 28) In Ts) 
comme 
29, 28, 28) 
In = fl 
(858) 


En se servant de la méthode À, il est commode, 
dans le cacul de l'erreur commise sur y, de profiter 
des graphiques des fonctions : 


a Il LS SNS) | 
TPS) 


et 
Ç | LSTSES) | 
or DS HAS SN 


On cherche alors les limites de l'erreur commise 
sur y correspondantes aux limites de l’erreur statis- 
tique du rapport respectif des nombres des coïnci- 
dences. En calculant, par exemple, l’erreur dans le 
cas du rapport «triples aux quadruples », on a, 


l'erreur relative sur ($, 89,9), &; — CS, 25, 2578 
sur (5, S, 5, 151) 
1 
SDS) NS Sos ER 
(S, 8, S) 


l'erreur relative sur le rapport =, 
(ASS ES) 


Ega —| 5% 9e 90 
re ED S2) 


1 
+ : et 
(AS, 28028) ER SES SES) 
et l’erreur absolue correspondante : 


(S3SS) 


(S, S, 5, ) RES 


On trouve ensuite sur le graphique mentionné 
les valeurs de (y + AY) correspondantes aux 
(S, 5, S)J(S, S, S, S) (1 dE Era). 

On procède de même façon dans le cas des 
rapports «triples aux quadruples » et «triples aux 
sextuples ». 

Dans une série des 112 segments de pellicule 
photographique nous avons obtenu les nombres 
suivants d'enregistrement des coïncidences de dif- 
férents types : 


JOURNAL DE PHYSIQUE 


N9o% 


) NOMBRE 
TYPE DE COINCIDENCE END EN ETS 
RARES) 39 745 
(516818) 37 619 
ET) 26 484% 
(SSSR ES ES) 20 538 
124 386 


(2.5, 2,5, 2) 
(nombre total) 


On a alors pour les coïncidences au moins triples, 
quadruples et quintuples : 


(S, W, S) Les 
RS 32 501 
ARS SES) 24 952 
(25, 25, 28) 2,583 + 0,007 
(S, 5, S) 
_{S, 8, 8) 1,482 + 0,007 
(S, 5, 8,5) 
_ISSS 1,923 & 0,01 
US, 8,8, #, 6) 
(5, S, $) 


pes 2,345 + 0,018 


(5, Se 15% 15 5, $) 
On obtient alors les valeurs suivantes de . 


MÉTHODE À 


MérHope B 


triples 
sextuples 


triples 
quintuples 


- triples 
quadruples 


1,369 + 0,004 1,398 + 0,010 1,432 + 0,008 1,454 + 0,008 - 


L’ensemble des donnéesstatistiques était examiné 
au moyen du test x? qui confirmait ses propriétés 
normales. 

De l’autre côté il faut introduire une correction 
de la valeur de y due au fait que la dérivée de y 
par rapport à la densité op est finie (Greisen, 1956). 
Pour cette raison le changement de la densité dans 
les méthodes À et B par un facteur donné est 
apparent, en réalité il est un peu plus petit à cause 
de la variation de y, ce qui entraîne la nécessité de 
la correction de la valeur de +. En supposant la 
variation de y en fonction de p proposée par 
Greisen, on obtient pour les valeurs corrigées de + : 
déterminée par la méthode B — y — 1,382, déter- 
minée par la méthode À (coïncidences triples aux 
sextuples) — y = 1,464. 


6. Conclusions. -— Les valeurs de y obtenues au 
moyen de la méthode À dépassent la valeur obtenue 
par la méthode B. Ce fait est.constaté en dehors des 
limites des erreurs statistiques. D'autre part on 


D No4 : 


observe l'accroissement de la valeur déterminée 
de y avec l’ordre de coïncidence employé. A notre 
avis, C’est une preuve de l’augmentation du rôle du 
gradient de densité dans l’enregistrement des 
gerbes par un dispositif de coïncidence d'ordre 
supérieur. Il est évident que la probabilité de 
l'enregistrement, par exemple, d’une coïncidence 
sextuple est diminuée par rapport à celle d’une 
coïncidence quadruple — à cause de la présence du 
gradient de densité. L'influence de ce facteur pro- 
voque un accroissement apparent de la valeur de y 
de 6 % dans le cas du rapport « sextuples aux 
triples », il faut pourtant se rendre compte que la 
disposition géométrique des bancs de compteurs 
réalisée dans cette expérience ne favorisait guère 


l'influence du gradient de densité (bancs juxta- 


posés). Si l’on s’imagine la disposition des bancs 
supplémentaires en dehors du triangle des trois 
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bancs (S, S, S), l'influence du gradient de densité 
doit augmenter. 

La probabilité d'intervention d’autres facteurs 
tel que le changement de densité pouvant causer 
l’accroissement de la valeur de y en passant des 
coïneidences triples aux sextuples est exclue : ce 
changement de densité des gerbes détectées qui 
passe d’environ 8,5/m? (méthode B) à 15,5/m? 
(triples aux sextuples) peut introduire un change- 
ment de y de 1,6 % au maximum. 

Nous remercions M. M. Miesowicz pour ses sug- 
gestions importantes, M. À. Tomaszewski pour les 
discussions, MIIes Z,. Chodkowska, Z. Kedzierska et 
M. T. Szubert pour le déchiffrement de pellicules 
photographiques, MM. C. Lis, P. Sulkowski et 
S. Kowalczyk pour leur aide technique. 


Manuscrit reçu le 12 novembre 1958. 
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Rosix (L.), Fonctions sphériques de Legendre et fonctions 
sphéroïdales. — Tome II (1 vol. 16 X 25 cm, vit + 
384 pages ; Gauthier-Villars, Paris, 1958, 5 000 F). — 
Collection technique et scientifique du C. N. E. T. 
Comme le tome Ï, paru à la fin de 1957, cet ouvrage, 

rédigé en vue des applications, s’adresse aux mathéma- 
ticiens, physiciens, chercheurs et ingénieurs qui ont à se 
servir des fonctions sphériques. On y trouve la plus grande 
variété de formules utiles, avec leurs domaines de vali- 
dité et souvent même leurs démonstrations. 

L'essentiel de ce deuxième volume est consacré aux 
fonctions associées de Legendre (Chap. IV et V) : formules 
de développements en séries, expressions sous forme d’in- 
tégrales définies ou d’intégrales de contour, calcul d’inté- 
grales définies contenant des fonctions de Legendre, re- 
lations entre fonctions de Legendre et fonctions de Bessel ; 
valeurs approchées, développements asymptotiques, iné- 
galités simples vérifiées par les fonctions associées de 
Legendre. Le dernier chapitre (Chap. VI) concerne les dé- 
veloppements en séries de fonctions sphériques : séries de 
polynômes de Legendre, séries de fonctions associées de 
Legendre, séries de polynômes de Laplace, avec leurs 


théorèmes de convergences et de sommabilités. 
M. Jesser. 


Faxo (U.) et Racax (G.), Séries tensorielles irréductibles. 
(Irreductible tensorial sets). Un vol. relié 15,5 X 23,5, 
de 171 p., Academic Press, N. Y., 1959, prix $ 6.80. 
En anglais. 

On appelle série tensorielle un ensemble de quantités 
qui subissent une transformation linéaire, quand on fait 
un changement de coordonnées — (chacune devient une 
fonction linéaire de l’ensemble). Cette notion englobe avec 
les tenseurs proprement dits, les coefficients de décompo- 
sition d’une fonction Ÿ de la mécanique quantique en ses 
composantes si la valeur propre correspondante est dégé- 
nérée, ses composantes elles-mêmes, etc. L'ouvrage com- 
mence par l'étude algébrique de cette notion, et la deuxième 
partie examine les applications à la mécanique quantique, 
puis des appendices. La notion de série tensorielle englobe 
celle de spineur ; la référence aux spineurs n’est faite que 
dans une annexe vers la fin de l’ouvrage. Y a-t-il avantage 
àrendre plus générales et moins précises les notions étudiées ? 
Le livre actuel constitue-t-il un progrès par rapport au 
classique ‘‘ Gruppentheoretische Methode in der Quanten- 
mechanik ” de Van der Waerden ? Un analyste peut seul 


répondre avec autorité. 
J. WINTER. 
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MESURE DU RENDEMENT DE FLUORESCENCE DE LA COUCHE K DANS L’ALUMINIUM 


Par F. SUZOR et G. CHARPAK, 


Laboratoire de Synthèse Atomique, GC. N. R. $. 


Résumé. -— La méthode décrite permet la mesure du rendement de fluorescence de la couche Æ 


des éléments de faible numéro atomique. Pour l’aluminium, on obtient Rx — 


(4,5 + 0,2) 10—2. 


On emploie la technique du compteur proportionnel. 


Abstract. — À method is described which allows precise measurement of Æ-shell fluorescence 


yield of low atomic number elements, 
Proportional counter technique is employed. 


I. Principe de la mesure. — Cette mesure est 
effectuée à l’aide du spectromètre nucléaire à 
compteurs proportionnels [1] que nous avons cons- 
truit. Dans le cas présent, seul un des deux comp- 
teurs est utilisé, et le principe de la mesure est le 
suivant : une source de fer 55 émettant des 
rayons X de 5,9 keV est placée sur la paroi du 
compteur et séparée du gaz uniquement par une 
feuille d'aluminium d’épaisseur a. Les dimensions 
du compteur sont assez grandes pour que le rayon- 
nement X soit pratiquement compté dans un angle 
solide 27. L’absorption des photons de 5,9 keV 
dans l'écran d’aluminium se produit principa- 


lement par effet photoélectrique dans la couche K 


et les rayons X de 1,5 keV qui en résultent sont 
aussi comptés dans le compteur avec un angle 
solide 2x. Le rapport y des intensités des photons 
de 1,5 keV et de 5,9 keV en fonction de l’épaisseur a 
de l’écran permet une mesure directe du rendement 
de fluorescence À4 dans l’aluminium. 

IT. Calcul théorique du rapport Y. — Traitons 
d’abord le cas de Pabsorption du rayonnement X 
dans une géométrie 2x. La source $ a une inten- 
sité Z dans un angle solide 4x ; nous allons calculer 


Pret 


l'intensité P(a) dans l'angle solide 27 après tra- 
versée de l’écran plan indéfini d'épaisseur a. u, étant 
le coefficient d'absorption du rayonnement, le flux 
sortant dans l’angle solide do compris entre les 


For aluminium the value found is Rr — (4,5 + 0,2) 10—?, 


cônes d'angle au somimel 0 et 0 + d0 est, 
traversée de lPécran 


après 


: sin Ô 6 ta4/cos0 0) 


4 


d’où en intégrant dans l’angle solide 2x. 
en posant 


Un photon primaire de 5,9 keV étant absorbé 
en M, dans l’écran d'aluminium, il existe une proba- 
bilité ? d'émission d’un photon X de 1,5 keV, dont 
le coeflicient d'absorption est 1,. Dans le volume 
élémentaire contenu dans l’angle solide do et com- 
pris entre les plans Let k + dA, il y a émission de 


7 sin 0 
E ;—lh/cos0 
PA LMD Ceres ON CU 


photons secondaires X de 1,5 keV. Compte tenu de 
Pabsorption de ces photons en géométrie 2% à 
travers l'écran d’épaisseur (4 — h), on obtient 
après intégration l’intensité S(a) des photons secon- 
daires sortant de l’écran d'aluminium 


S{ah = gun 6 [° Hluste— n) eus 2) an | 


£(z) f exe (dx fx). 


Ayant par définition y(a) = S(a)[P(a), il est 
alors facile de calculer, point par point, la courbe 
y(a)[e en fonction de 4, sachant que , —120 em?/g 
pour les photons de 5,9 keV et que us — 330 em?/g 
pour les photons de 1,5 keV ; la forme de cette 
courbe dépend seulement de y, et de u. 


en posant 


III. Détermination expérimentale du rapport y. 
— Les sources de fer 55 sur les différents écrans 
d'aluminium ont été obtenues par vaporisation 


1 


No 


sous vide. Ces sources sont réparties uniformément 
au centre de l'écran d'aluminium et contiennent 
une quantité assez faible de matière pour que les 
éiectrons Auger de 5,3 keV émis par le fer 55 
apparaissent sous forme d’une raie peu élargie et 
peu abaissée en énergie, comme le montre le spectre 
reproduit sur la figure 2 ; celui-ci a été obtenu 


Electrons Auger = 


5,3 keV 


50 Remplissage du compteur 

2 Argon : 60 cm de mercure 

© 
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avec une source directement en contact avec le gaz 


- du compteur. 


Dans les expériences destinées à la mesure du 
coefficient de fluorescence, ces électrons Auger 
n'apparaissent plus, étant absorbés par l’écran 
d'aluminium. Les impulsions données par le comp- 
teur sont, après amplification, analysées dans un 
sélecteur à 50 canaux ; la forme du spectre obtenu 
pour un écran d’aluminium, a — 0,78 mg/em? est 
donnée à titre d’exemple sur la figure 3. Ce spectre 
correspond à un remplissage du compteur par du 
propane sous une pression de 57 cm de mercure ; 


Nombre d'impulsions 


RENDEMENT DE FLUORESCENCE DE 


LA COUCHE. K AG 
leflicacité pour les photons de 1,5 keV est alors 
égale à l’unité et celle pour les photons de 5,9 keV 
est égale à 0,142 ; cette valeur a été mesurée par 
comparaison avec le cas où le remplissage du comp- 
teur était fait avec de l’argon sous 60 cm de mer- 
cure et du propane sous 6 em, l'efficacité des pho- 
tons de 5,9 keV étant alors égale à l’unité ; en outre 
cette valeur expérimentale est en accord avec la 
valeur calculée connaissant les dimensions du comp- 
teur. La valeur de P(a) est déduite de l’aire de la 
raie de 5,9 keV sur le spectre de gauche, compte 
tenu de l'efficacité indiquée ; le spectre de droite 
avec une amplification double fournit avec pré- 
cision Ja valeur de S(a), déduite de Paire de la raie 
de 1,5 keV. La présence d’un fond n’est pas gênante 
et ne nuit pas à la précision de la détermination 
de S(a), l’interpolation du fond étant choisie de 
façon à obtenir une raie symétrique et de résolution 
connue ; la forme de ce fond a été contrôlée en 
remplaçant l’écran d'aluminium par un écran de 
béryllium. 

La mesure de y(a) a été faite pour plusieurs 
valeurs de a. À ce sujet, nous devons mentionner 
que des résultats cohérents et reproductibles n’ont 
pu être obtenus que si l’écran d'épaisseur a était 
constitué d’une seule feuille d'aluminium ; la pré- 
sence de plusieurs feuilles d'aluminium, entraînant 
des boursouflures et des vides inévitables entre les 
feuilles, avait pour effet de détruire les conditions 
simples d’une géométrie 2x. 


IV. Détermination du rendement de fluorescence. 


— Les valeurs expérimentale de y sont reportées 


Remplissage du compteur 


Propane : 57cm de mercure 


a = 0,78 mg/cm° d'Al 


sur la figure 4 ainsi*que la courbe théorique corres- 
pondant à p — 4,1.10-? donnant le meilleur accord 
avec les points expérimentaux. 

En possession de cette valeur de p, il nous faut 
maintenant obtenir le coefficient de fluorescence Rx. 


À partir de la discontinuité K de la courbe repré- 
sentant le coeflicient d’absorption des photons 
dans l’aluminium en fonction de l’énergie, nous 
tirons la valeur N£/N — 0,92 pour la probabilité 
d’ionisation dans la couche X. La diffusion con- 


tribue pour 1,3 % dans la section eflicace d’absorp- 
tion en bonne géométrie du rayonnement X du 


10! 7. 10° 


—— Courbe théorique pour CELAP 107% 


° Points expérimentaux 


mg/cm? Al 


Xe. «9: 


fer 55 dans l’aluminium ; comme dans notre cas 
une partie du rayonnement diffusé contribue aussi 
à l’ionisation de la couche Æ de l’aluminium, nous 
négligerons cet effet qui est certainement inférieur 
aux erreurs expérimentales. Nous négligerons aussi, 
parce qu’il est d’un ordre de grandeur très infé- 
rieur, l’ionisation dans la couche Æ par des élec- 


ÉLÉMENT 7 
Û 6 0.09 [6] 
N 7 0,15 [6] 
O 8 0,22 [6] 8,2 [10] 
Ne 10 81 [6] 8,3 [10] 
Mg 1192 
Al 3 
Si 1% 
S 16 145 [8] 
CI 17 4 [7] 
A 18 7 [9] 14,910] 
1 R OS 
+ Heintze rs 
o Rubinstein, Snyder 
10 + Charpak , Suzor 
5 
Z {= 
Û 5 . 20 


5 10 15 
F Ne Na Mg A1 Si P S CL A K 


Fross 


Nous avons représenté sur la figure 5, la courbe 
semi-empirique obtenue avec les valeurs de À, B 
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trons Auger de 5,3 keV ou des photoélectrons 
de 4,3 keV. Nous avons, dans ces conditions, la 


relation : 
Rr = p10,92 


qui nous donne : 
Re = (4,5 + 0,2)40-: 


V. Discussion générale sur le rendement de 
fluorescence Rx.— La valeur théorique du rende- 
ment de fluorescence Rx [2, 3, 4] est de la forme 


( Rx 


1 
Re ru 173 
SN 4 + 82 + C2 


où À correspond à une correction due à l'effet 
d'écran, CZ3 à une correction relativiste, tandis que 
le terme principal BZ traduit la dépendance en Z4 
des transitions dépolaires. Différents auteurs ont 
publié des valeurs des coefficients 4, B et C telles 
que la formule précédente rende le mieux compte 
des valeurs expérimentales connues de Rx. La der- 


nière publication en date est celle de Laberigue- 


Frolow et Radvanyi [5] qui donne : 


A=—0,0217 B= +0,03318 C—=—1,14.107$ 


l’accord avec l'expérience est satisfaisant pour 
Z->02p: 

Pour les valeurs de Z inférieures à 20, les résul- 
tats expérimentaux sont très dispersés comme le 
montre le tableau suivant : 


Rx .10? 
1,8 [15] 
15341017] 0,9 [17] 
0,8 [17] 
3,8 [11] 
8,3 [11] 
10,8 [11] 
TNT 6 [13] 12,3 M%] 42;92F5] 


et C précédemment citées et nous y avons reporté 
en plus de notre détermination de Rx pour l’alumi- 
nium, les valeurs récentes obtenues par Heintze[15] 
qui a employé comme nous la technique du comp- 
teur proportionnel ; nous avons ajouté aussi la 
valeur obtenue pour l’argon par Rubinstein et 
Snyder [16] à l’aide d’un calcul théorique rigou- 
reux. Tous ces résultats sont en accord pour indi- 
quer des valeurs supérieures à celles de la courbe 
semi-empirique, et s’en écartant d’autant plus en 
valeur relative que Z est plus petit. 

Ce travail a été accompli grâce aux moyens mis à 
notre disposition par le Centre National de la 
Recherche Scientifique et à une subvention du 
Commissariat à l'Énergie Atomique. 


Manuscrit reçu le 29 octobre 1958. 
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REVUE DES LIVRES 


MacConwezz (J.), Dynamique des particules quantiques 
(Quantum particle dynamics). Un vol. relié de 252 p., 
16 X 28 cm, North-Holland Publ, Company, Ams- 
terdam, 1958. En anglais, prix F1. 20. 

Un excellent ouvrage d’enseignement, clair et de lecture 
très agréable. Il couvre la mécanique quantique ordinaire, 
la théorie de Dirac, l’application aux phénomènes de 
rayonnement, les théories de l’effet photoélectrique et de 
lPémission des rayons X par freinage, et enfin aboutit 
à la théorie des forces nucléaires, et au calcul des phéno- 
mènes de haute énergie et à la théorie des mésons. Le der- 
nier chapitre est consacré à la diffusion des mésons, pro- 
duction et annihilation des pions, des protons négatifs. 
Les éléments de matrice donnant les probabilités de ces 
phénomènes sont calculés. Tout cela en 250 pages, en 
partant de la relativité restreinte et de la théorie de l’onde 
de Broglie. Nous ne saurions assez recommander la lecture 
d’un tel ouvrage, qui rendra service à tous, professeurs 
et étudiants. 

J. WiINTER. 


Doremus (R. H.), RoserTs (B. W.) et TurnBuLz (D.), 
Croissance et perfection des cristaux. (Growth and Per- 
fection of Crystals.) 1 vol. 21 x 28 cm, 609 pages, 
John Wiley and Sons, Inc., New-York, 1958, $ 12,50. 
À une époque où les cristaux jouent un rôle de plus en 

plus grand dans les diverses sciences, toute conférence 

internationale sur leur croissance et leur perfection est un 
événement important. Le présent recueil est la transcrip- 
tion des communications et discussions présentées à la 

Conférence internationale de Cooperstown (N. Y.) du 

27 au 29 août 1958 par les spécialistes des États-Unis, du 

Royaume-Uni, d'Allemagne et de Belgique. 

Un grand nombre de papiers concernent le problème de 
l’origine des imperfections linéaires dans la croissance 
cristalline. À ce propos des résultats récents ont été pré- 
sentés sur les « poils » cristallins (whiskers) dont les pro- 
priétés élastiques sont « parfaites ». 

D’autres communications décrivent les mécanismes par 
lesquels des dislocations peuvent être introduites dans les 
cristaux, soit par précipitation de vacances, soit durant 
la cristallisation même. 

Certains auteurs indiquent comment on peut faire 
croître des grands cristaux apparemment dénués de dis- 
location. 

Les conférenciers remarquent que pour comprendre le 
problème de la croissance cristalline, de gros efforts doivent 
encore être portés sur l’étude de l’interface liquide-solide 


et de la nature de l’état liquide. Quelques remarques sont 
présentées à ce sujet. 

De nouvelles explications sont présentées sur le problème 
de la croissance dendritique, cependant certains méca- 
nismes restent encore obscurs. 

Quelques exposés commencent à jeter une lumière sur - 
le rôle obscur joué par les impuretés sur la vitesse de 
croissance et sur la morphologie du cristal. 

Enfin les cristallographes des hauts polymères ont pu 
confronter leur expérience avec leurs collègues des cris- 
taux plus usuels. 

Signalons que les participants ont tenu à présenter leurs 
plus beaux clichés, les non spécialistes pourront les ad- 
mirer avec profit dans cette édition remarquablement 
présentée. | 

Un seul reproche technique: quelques défauts de brochage 
(dans l’exemplaire qui nous a été confié, une dizaine de 
pages étaient mal découpées, dans un autre, un cahier de 
trente-deux pages manquait tandis qu’un autre était en 
double). Admirons pour terminer avec quelle rapidité la 
General Electric s’est chargée de publier cette impor- 
tante conférence. 

Y. LE Core. 


HerzBercer (M.), Optique géométrique moderne, (Moderne 
geometrical Optics) (en anglais). Un vol. relié de 504 p., 
15 x 24 cm, Interscience publishers, New-York, 1958, 
prix $ 15. 

Cet ouvrage est un traité conçu pour les réalisateurs. 
Le formalisme mathématique est simple, et expliqué en 
fin d'ouvrage par un appendice. Les formules mathéma- 
tiques employées sont calculables par machine électronique 
à grande vitesse. Partie I : Marche d’un rayon à travers 
un système optique. Partie IT : Optique de Gauss. Appro- 
ximation du premier ordre. On introduit un nouveau 
symbole : Les crochets de Gauss. Partie III : Lois générales, 
principes de Hamilton et Lagrange. Partie IV : Systèmes 
concentriques. Partie V : Systèmes ayant la symétrie de 
rotation autour d’un axe. Partie VI: Théorie approxi- 
mative pour un système normal. Partie VIT : Théorie des 
erreurs sur l’image des troisième et cinquième ordre. 
Partie VIII : On développe un modèle mathématique ; 
On trace certains rayons, méridiens et obliques, ce qui donne 
une formule valable au 5° ordre, puis on l'utilise comme 
formule d’interpolation. Partie IX: Milieux non homogènes. 
Partie X : Appendice mathématique. La présentation 


est claire. 
J. WINTER. 


30 


LE JOURNAL DE PHYSIQUE ET LE RADIUM 


TomE 20, AVRIL 1959, PAGE 466. 


DÉPLACEMENTS ISOTOPIQUES DANS LE SPECTRE D’ARC DU TUNGSTÈNE 


Par JEAN BLAISE et GILBERT GLUCK, 
Laboratoire Aimé-Cotton, CG. N. R. $., Bellevue (Seine-et-Oise). 


Résumé. — De la raie 5 642 À (dis? 5), — d* sp *F$), on a déduit le déplacement isotopique 
relatif des 4 isotopes abondants de W. La position relative de T$0W par rapport à TVR A LEE 
a été mesurée dans 4 raies intenses du type ds? — d# sp. Les positions relatives des isotopes sont : 


(186) M OEAMSZL)ENTE 


(183) :1,54 + 0,02 ; (182) : 2,13 + 0,02 ; (180) : 2,92 + 0,05. La position 


anormale de 182W entraîne l’inversion du phénomène de « odd-even staggering ». | 

La mesure du déplacement isotopique de plusieurs centaines de raies permet de poursuivre la 
classification du spectre d’arc où 6 types de transitions ont déjà été reconnus. 

La valeur de la constante de déplacement isotopique pour la paire d’isotopes 18#W — 188W 


CMOS Mig SE né 


Abstract. — From the line 5 642 À (dä5s2 1), — d# sp *F$), the relative isotope shift of the 


four abundant isotopes of W has been determined. 
to 186W and 181W has been measured in four strong lines of the type dis? — d sp. 
(182) EHESS) EM EN 002482) AS ER 02 
The abnormal position of 18W implies the inversion of the phenomenon 


positions of the isotopes are : (186) : 0 ; 
(180) : 2,92 + 0,05. 
of odd-even staggering. 


The relative position of 1#0W with respect 
The relative 


The measurement of the isotope shift of several hundred lines allows a further classification of 
the arc spectrum where 6 types of transitions are already recognized. 
The value of the isotope shift constant for the pair of isotopes 184W — 18W js 


B Cox = 117 + 20 mK. 


I. Introduction. — Le déplacement isotopique 
dans le spectre d’arc du tungstène a été étudié dans 
plusieurs laboratoires depuis une douzaine 
d'années, en utilisant soit du tungstène naturel, 
soit du tungstène enrichi en 188W ou en 18W. 
D’après Williams et Yuster [1], les abondances 
relatives des isotopes du tungstène naturel sont les 
suivantes : 


TABLEAU I 
Isotopes 180 182 183 184 186 
Abondances 0,135 26,4 13,4 30,6 28,4 


Les abondances des isotopes pairs 182, 184 et 
186, apparaissent très favorables pour faire des 
mesures de déplacements isotopiques. Malheu- 
reusement, l’isotope impair a un spin 1/2 et se 
sépare donc en deux, trois ou quatre composantes 
suivant la transition étudiée, et dans la grande 
majorité des cas, les deux composantes les plus 
intenses coïincident approximativement l’une avec 
la composante de 182W, l’autre avec celle de 184W. 
Il en résulte que la mesure de la séparation (186- 
184) et surtout celle de la séparation (184-182) 
devient très difficile et que le déplacement isoto- 
pique relatif mesuré entre les trois isotopes risque 
donc fort de varier avec la raie étudiée. Quant à 
lisotope 180, deux cent fois moins intense que son 
voisin immédiat, il pose lui aussi un problème 
délicat si l’on ne dispose pas d’isotope enrichi ou 
d’un dispositif interférométrique à grand contraste, 


Toutes les raies dont le déplacement isotopique 
a été mesuré avec précision correspondent à des 
transitions 

5d* 6s? — 5d46s6p ou 5d5 6s — 5d8 6s? Gp, 
ce dernier type de transition ayant été reconnu par 
Murakawa [21] en raison de la grandeur et du sens 
du déplacement isotopique observé. 

Tandis que Kopfermann et Meyer [3] et Mura- 
kawa [2], [4] ont utilisé du tungstène naturel, 
Vreeland et Murakawa, dans un travail effectué 
en 1951 [5] et repris en 1958 [6], utilisèrent des 
échantillons enrichis en 188W et en 18W et Barr 
et Jenkins [7] employèrent un échantillon enrichi 
en 18W, Nous avons rassemblé dans le tableau II 
les déplacements isotopiques relatifs mesurés par 
ces auteurs, en prenant comme unité la sépa- 
ration (184-186) ; la précision a étéindiquée chaque 
fois que les données des auteurs le permettaient. 


TABLEAU II 


AUTEURS 186 184 183 182 180 
K.M. (4947) © 1 D A1 14 
V. M. (1951) O0 1 2,13 3,15 
M. (1953) 0 1 2,206 + 0,010 
M. (1956) 0 . ‘1 2,206 + 0,010 3,08 + 0,07 
BJ.(1956) 0 1 249% 2,91 
V. M. (1958) O0 1 2,20, + 0,03 2,98 + 0,40 


On voit que la moitié des résultats indique que la 
séparation (184-182) est supérieure de 11 à 14 VA 
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à la séparation (186-184), tandis que l’autre moitié 
donne un écart de 21 % entre les deux séparations. 
Pour la position relative de l’isotope 180, la dis- 
persion des résultats est du même ordre. Quant à 
lisotope 183, Kopfermann et Meyer pensent qu’en 
raison du phénomène de « odd-even staggering », 
il doit se trouver plus près de l’isotope 182 que de 
lisotope 184 ; Murakawa se basant sur son analyse 
de la raie 4 269 À (d5s 7%, — d8 s2 p 263) estime 
qu'il est équidistant des deux isotopes pairs qui 
l’encadrent. 

Dans ces conditions, il nous a paru intéressant 
d'étudier une raie du spectre d’arc dans laquelle 
l’isotope impair ne présente pas de structure hyper- 
fine. 11 est bien connu que les termes spectraux 
appartenant à des configurations d®s?2 ont des 
structures hyperfines négligeables et que dans le 
cas du tungstène, les structures hyperfines obser- 
vées doivent être attribuées aux termes des confi- 
gurations d* sp ou d5s. Il suffit donc d’étudier une 
transition entre un terme de la configuration d' sp 
pour lequel J —0 et un terme d4 s? ayant un J —1. 
D’après la classification de Laporte et Mack [8], la 
raie la plus intense de ce type est la raie 5 642 À 
(d4 82 5D, — d* sp 7F5) : c’est donc celle que nous 
avons étudiée. 


IT. Résultats expérimentaux. — 1° DÉpLa- 
CEMENTS ISOTOPIQUES RELATIFS. — Nous avons 
utilisé une cathode creuse de Schüler démontable. 
On a introduit une mince feuille de tungstène 
enroulée sur elle-même dans le pot cathodique en 
duralumin ; en se déroulant, la feuille s’applique 
sur la paroi interne du tube assurant un bon 
contact thermique. Un petit disque de tungstène 
a été ensuite poussé à force jusqu’au fond du pot 
cathodique. En utilisant le néon comme gaz support 
de la décharge, une intensité de courant de 40 mA 
suffit pour exciter la quasi totalité des raies classées 
par Laporte et Mack, la cathode étant refroidie à 
l’azote liquide. 

Les enregistrements de la raie ont été effectués à 
l’aide du spectromètre Fabry-Perot [9]. Nous 
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tenons à remercier M. le P* D. A. Jackson qui nous 
a aimablement prêté un étalon d’une qualité 
optique remarquable. 

La figure 1 montre un enregistrement obtenu 
avec une constante de temps de 3,3 sec, la durée de 
balayage d’un ordre de 36 minutes. 


BECOME 


La largeur à mi-hauteur des quatre composantes 
visibles est identique et égale à 18,0 mK 
(1 mK = 1075 cm1); quant à leurs intensités 
relatives, elles sont exactement celles indiquées par 
Williams et Yuster [1], ce qui permet de conclure 
que les deux composantes de l’isotope 183 coïn- 
cident et que, comme prévu, le terme 5d4 6s? 5D, a 
une structure hyperfine nulle. 

Dans ces conditions, la mesure des distances 
entre composantes fournit directement les dépla- 
cements isotopiques relatifs des quatre isotopes les 
plus importants du tungstène. En raison de sa 
faible intensité, il n’est évidemment pas question 
de déceler l’isotope 180 qui, avec l’épaisseur d’éta- 
lon choisie, tombe sous la composante 18W d’ordre 
d’interférence inférieur. Le tableau IIT indique les 
résultats obtenus. 


TABLEAU TITI 


Isotopes 186 
Position des composantes (en mK) 0 
Déplacements isotopiques relatifs 0 


Ces résultats ont été confirmés par les mesures 
que nous avons faites sur 22 raies intenses appar- 
tenant aux multiplets 5D — 7F0, 5D — TD° 
et 5D — 7P° et dont la figure 2 montre un 
exemple. 

Nous avons en outre étudié les raies 4 983 À, 


184 183 182 
BE01200 6 082,50 2005 117,0 + 0,3 
1 1,54 + 0,02 2,13 + 0,02 


5 015 À, 5 069 À et 5 225 À avec un étalon donnant 
un intervalle entre ordres de 240 mK. Avec une 
intensité de courant de 15 mA, ces raies deviennent 
si intenses que la distance (186-180) devient faci- 
lement mesurable avec une précision de + 2 mK 
(la hauteur de la composante de 1#W est de 8 mm 
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sur notre enregistrement de la raie 4 983 effectuée 
dans ces conditions). Pour les quatre raies le rap- 
port (486-180) : (186-184) — 2,92 + 0,05. 
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Pour l’ensemble des isotopes stables du tungs- 
tène, les déplacements isotopiques relatifs sont donc 
les suivants : 


TABLEAU IV 


Isotopes 186 


Déplacements isotopiques relatifs 0 


La comparaison avec le tableau I montre le 
complet accord avec les résultats de Barr et Jenkins 
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pour tous les isotopes pairs et avec Kopfermann et 
Meyer et les premiers résultats publiés par Vreeland 
et Murakawa pour les isotopes 182, 184 et 186. 


29 POSITION RELATIVE DE L’ISOTOPE IMPAIR PAR 
RAPPORT AUX ISOTOPES PAIRS. — La séparation 
(182-183) est supérieure de 10 % à la sépara- 
tion (183-184), alors que dans tous les éléments 
étudiés jusqu’à présent, on constate le déplacement 
de l’isotope impair vers l’isotope pair léger. En fait, 
c’est la position relative de lisotope 182 par 
rapport aux autres 1sotopes pairs qui est anormale. 
Elle indique que le carré du paramètre de défor- 
mation du noyau décroît beaucoup plus rapidement 
entre les isotopes 180 et 182 qu'entre les iso- 
topes 182 et 184 ou 184 et 186. De plus, comme 


(183-184) > à (184-186), cela confirme bien la posi- 


tion anormale de l’isotope 182 par rapport à tous 
les autres isotopes. 


ITT. Déplacements isotopiques des termes spec- 
traux du spectre d’arc : Résultats préliminaires. — 
Nous avons entrepris en 1957 de contrôler et 


184 


180 
9,92 + 0,05 


182 
2134 0,02 


183 
1,5% + 0,02 


d'étendre la classification du spectre d’arc du tung- 
stène en utilisant les résultats fournis par l’étude 
des déplacements isotopiques. Dans ce but, il suffit 
de s'intéresser au déplacement isotopique des deux 
isotopes les plus abondants 184W et 18W. Nous 
nous contenterons de donner dans les paragraphes 
suivants, un résumé des résultats déjà obtenus, les 
résultats complets devant être publiés ultérieu- 
rement. 

Nous avons déjà mesuré le déplacement isoto- 
pique de plusieurs centaines de raies entre 3 700 et 
6 000 À. Parmi ces raies, les deux tiers environ ont 
été classées par Laporte et Mack [8] et leur classi- 
fication a été généralement confirmée. Quand plu- 
sieurs classifications sont proposées pour la même 
raie la mesure du déplacement isotopique permet 
de choisir celle qui convient. Les nouveaux niveaux 
indiqués dans Atomic Energy Levels [10] ont per- 
mis de retrouver la classification d’un certain 
nombre de raies. Enfin nous avons trouvé de 
nouveaux termes spectraux 1inpairs, mais Jusqu'à 
présent, aucun nouveau terme appartenant aux 
configurations 5d4 65? ou 5d5 6s. 

Dans les dix raies mesurées du multiplet 
5d® 6s? 5D — 5d*6s6p F0, le déplacement isoto- 
pique ne varie que de — 52 à — 55 mK, ce qui 
permet de conclure à l’absence de perturbations du 
multiplet fondamental et du multiplet impair le 
plus profond par des termes d’autres configurations 
électroniques et nous considérons donc la valeur 


00(184-186) — — 55 mK comme caractéristique 
des transitions ds? — d# sp. Des déplacements 
isotopiques — 45 mK mesuré dans la raie 


4074 À (dts25D, — dtsp?P;) et +25 mK 
mesuré dans la raie 4 295 À (d5s 75, — q4 sp *P5), 
on déduit la relation ST(d4 52) — ST(d7 s) — 69 mK. 
Si l’on suppose en outre que le rapport 

dT(d* sp) 


DT(d5 s) 


comme dans les éléments voisins [11], [12], on est 
conduit à attribuer à la configuration 5d5 6s un 
déplacement isotopique de + 70 mK par rapport 
à la configuration 545. 

Nous avons indiqué dans le tableau V les dépla- 
cements 1sotopiques caractéristiques des différents 
types des transitions étudiés. 


r 1,20 
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TABLEAU V 
TRANSITION Ôo (184-186) 
5d4 682? — 544 6s6p "i 55 
5d5 6s — 5d#* 6s6p + 14 
5d5 6s — 5d$ 6s? 6p + 8 
544 68? — 54365? 6p + 12 
5d4 6s6p — 5d46s7s mm + 10 
54? 65? 6p — 5d4 6s7s m — &0 


IV. Attribution de configurations électroniques 
aux termes speetraux. — Toutes les configurations 
électroniques indiquées dans Atomic Energy 
Levels [10] sont confirmées à l'exception des sui- 
vantes : 

1° Les termes 19 535,04 K, et 19 826,04 K clas- 


_ sés respectivement d5 s5 G;et f4s? 3G; doivent être 


permutés et l’ordre des niveaux du quintuplet est 
maintenant 5G:, 5°G:, 5Gx, 5G, et 5G; au lieu de 
l’ordre normal. Les valeurs trouvées pour les dépla- 
cements isotopiques des termes d4 s? 3G;et d5 s 5G, 
confirment leur perturbation mutuelle. 

2° Les raies ultimes du spectre d’arc ne sont pas 
des transitions d5s %S, — d5 p ?P534, mais bien 
des transitions d5 s — d sp, comme le pensaient 
Laporte et Mack [8]. Le terme 7P% est même légè- 
rement perturbé par un terme de la configu- 
ration dÿ 5? p. 

30 La configuration dÿ s?p fait son apparition 
avec le terme 25983,60 K, classé antérieu- 
rement d*sp 5F,. De nombreux autres termes de 
cette configuration ont été identifiés. Par contre, 
aucun terme de la configuration d5 p n’a été trouvé 
jusqu’à présent. 


V. Constante de déplacement isotopique. — Le 
déplacement isotopique de 70 mK que nous avons 
attribué à la configuration d5s pour le couple 
d’isotopes 184-186 est précisément la valeur 
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admise antérieurement par Brix et Kopfer- 
mann [13] pour le couplage d’isotopes 182-184. 
En utilisant la valeur de 
1 d'a 


né(6s) dn = 0,596 + 0,06 


évaluée par ces auteurs et les déplacements isoto- 
piques relatifs déterminés au début de cet article, 
on peut calculer la constante de déplacement isoto- 
pique expérimentale pour chaque couple d’isotopes 
pairs définie par la formule 


B 01e, 
PSN 
n3(6s) dn 


B Cexn 


et la comparer à la constante de déplacement isoto- 
pique dû à l'effet de volume, Cw, calculée pour 
R5 —1,20.107 18 em, 


TABLEAU VI 
COUPLES 
D’ISOTOPES 180-182 182-184 184-186 
B Cexp (MK) 93 133 117 
Cin (mK) 310 308 306 
Cexp /Cth 0,30 0,48 0,38 


Il semble que les valeurs de 8 Ce, puissent être 
considérées comme correctes à mieux de 20 près. 

Nous tenons à exprimer notre reconnaissance à 
M. le PT P. Jacquinot sous la direction duquel ce 
travail a été effectué. Nos remerciements s’adres 
sent également à Mme M. Diringer-Louvegnies et 
à M. S. Gerstenkorn pour le concours qu’ils nous 
ont apporté dans la partie expérimentale de ce 
travail et à MM. P. Giacomo et S. Roizenqui ont 
préparé les couches diélectriques des étalons que 
nous avons utilisés. 


Manuscrit reçu le 27 décembre 1958. 
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ÉTATS D’UNE PAIRE D’ISOBARES 


Par R: NATAF, 


Collège de France, Orsay. 


Résumé. — On étudie la dégénérescence des états de deux particules isobares et de leurs anti- 
particules, en l’absence de champ électromagnétique : 
1. Lorsque ces particules sont des Fermions (nucléons). 


2. Lorsque ce sont des Bosons (mésons Æ). 


La quantification des champs associés résultant de règles de commutation différentes, la dégéné- 
rescence est d’ordre plus élevé pour les Bosons ; elle porte sur les quatre états de particules et 
antiparticules. Les transformations possibles entre ces 4 états comprennent un sous-groupe homo- 
morphe des rotations dans un espace à 4 dimensions ayant la métrique de Minkowski. 

3. On compare enfin ces dégénérescences à celle qui résulte de l'ignorance de l’état de « charge » 
pour la paire neutrino-antineutrino à quatre états. 


Abstract. — The degeneracy of the states of two isobaric particles (with their antiparticles), in the 
absence of any electromagnetic field, is analyzed within the frame of present field theory : 
1. When these particles are Fermions (nucleons). 


2. When they are Bosons (Æ mesons). 


The quantization of the associated fields being then derived from different commutation rules, 
the degeneracy for Bosons has a higher order: it involves the four states of particles and antiparticles. 
The group of transformations between these states contains a sub-group homomorphic to rota- 
tions in a four dimension space with Minkowski metric. 

3. These denegeracies are compared with that resulting from ‘* charge ”’ uncertainty on the four 


states of the neutrino-antineutrino pair. 


A. États d’une paire de fermions isobares. — 


19 Supposons que deux nucléons soient placés dans 
une enceinte isolée E, où il n’y a pas de champ 
électromagnétique appliqué. Ces deux nucléons 
semblent identiques : leur différence de masse est 
inférieure à l'erreur des mesures que l’on peut 
effectuer dans E. Cependant, l'étude expérimen- 
tale de leurs états possibles (au cours de leurs colli- 
sions, par exemple) montre qu’ils peuvent se 
trouver tous deux dans le même état quantique 
d'espace et de spin. Il faut donc admettre qu’ils 
correspondent à des états «intrinsèques » différents, 
soient p et n, orthonormaux. Comme nous ne 
savons pas en quoi diffèrent ces états, nous pou- 
vons les remplacer par : 

p° = ap + bn 


1.6 


n = Cp + dn (1) 


eux-mêmes orthonormaux. 

z > Het oo 

Si nous représentons p par Lo , An par (4 ,» et 
posons en général : 


C=tp+tn(s) 


2 
l'orthogonalité de €, n se traduit par : Ci n —0, 
C étant normé par Cf € — 1. 
On voit alors que : 


doit être unitaire, 


| [p'12 = 4 imposant 24] ER = (2) 
In'|. =1 imposant tel al (3) 
p'fn’ = 0 imposant : ca* + db* — 0 (4) 


(1) correspond à une dégénérescence des états de la 


\ 


paire isobare (?) en l’absence de champ électro- 


magnétique, qui est levée par l'application d’un tel 
champ : on distingue alors sans ambiguïté l’état p 
chargé et l’état n non chargé. 

L'intérêt de (1) provient, on le sait, du fait que la 
dégénérescence n’est pas levée par les forces 
nucléaires ; ainsi, les nucléons d’un noyau léger 
sont pratiquement dans les conditions de E. 

Mais on peut considérer les transformations uni- 
taires (1), par le même raisonnement, pour une 
paire quelconque d’isobares qu’il est impossible de 
distinguer dans des conditions données — en parti- 
culier pour une paire de Bosons isobares de charges 
différentes, placés dans E : la fonction d’onde de 
l’ensemble n'étant pas toujours symétrique (d’es- 
pace), on pourrait voir encore qu’il y a deux états 
intrinsèques différents. 


29 Le groupe des transformations unitaires U de 
(1) définit le spin isobarique. — En effet, substi- 
tuant : d = — c(ab—1)* (5) de (4) dans (3), et 
tenant compte de (2), [| —{b| ; le groupe des 
U;, unitaires à deux dimensions, dépend donc de 


No 4 


4 paramètres réels, (2, 3, 4) étant vérifiées pour : 


a — eu cos À ; b — eiv sin À; c — — ew sin à, et (5). 


PONT Er 


Avec ces paramètres : 


x es ( eiu cos À eivsin À ] 

È — eiw sin À  eilw+v—u) cos À 

x 

3 .U+W i(u— iT ù 

4 rs sa cos À e cos À | 

“4 S . W—v . (V+W 

E De i —ù 

; 7e 2wsinx ( # ) cos À 
, Dre D + w 

… c’est-à-dire, en posant : 0 — = : 

| p' = eï{a; p + bin) 

i , 

n' = ei8(— pr p aï n) 

4 

£ ) Tv b « 

D où : = | : = a un déterminant 1. 


Les matrices V, de déterminant 1 forment, on le 
sait, le sous-groupe unitaire unimodulaire de U,, 
à 3 paramètres réels, homomorphe du groupe À; 
des rotations autour d’un point. 

V prend la forme classique : 


nl alé 5 oi || de à 


52 . 


its Ê 
ONCE 


(EL 
ep (6) 
- en fonction des matrices +; de Pauli et des «angles 
. de rotation » «fy tels que : 
D EG BEN ED NET VIS 
DIR TU: TR ER Te 


u 


De et : 
U — ei0 p. (7) 

30 Cependant, à chacune des particules de la 
paire isobare, on peut associer son antiparticule, et, 
pour bien préciser la relation de « conjugaison de 
charge », nous nous restreindrons au Cas d’une 
paire de Fermions isobares. 

Dans E, en l’absence de conditions de seuil 
d'énergie, on ne pourrait pas savoir en quoi dif- 
fèrent les 4 états : p, n, p°, n°. Mais ces 4 états ne 
sont pas tous indépendants, et pour bien voir leurs 
relations, il faut examiner les transformation du 
type (1) dans le cadre de la théorie quantique des 
champs. 

a) On peut toujours effectuer (1) sur les 
champs p, n initiaux que nous savons indépen- 
dants (sans savoir en quoi ils diffèrent) ; nous pou- 
vons supposer que ces champs anticommutent : 


ns, p}= {tn} = 0 
tandis que : 
Epi,mi=de, non} 


‘(F fonction des coordonnées) 


{pi ,pi= {pri} = 0 
et de même pour n. (8) 
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En imposant les mêmes relations d’anticom- 
mutation aux champs p', n’ on retrouve alors les 
conditions (2) : relation (8) pour p’, (3) : relation (8) 
pour n’, et (4) : anticommutation de p’et n’. 

b) On ne peut pas effectuer de transformation 
du type : p” — Ap + Bp°. Pauli [1] a montré que 
seule une transformation : 


p" = Ap + B y;p° (9) 


pourrait conserver les relations d’anticommutation, 
et celle-ci n’est pas compatible avec les équations 
d'évolution de Dirac : 


(Yu du — Mo) p = 0 (Yu du — Mo) n = 0 (10) 


que vérifient les champs p, n, ainsi que p°, n°, en 
général. 
Par contre, des transformations : 


ap + bnc 
n'© — cp + dn° 


TS 
I 


de matrice Ü, unitaire sont possibles, les champs p 
et n° étant indépendants. 
En posant : 


ile) 


on a donc les transformations possibles : 


(11) É=AURE n = Un (12) 
d’où : 
(13) Ee = U*Ee n° = Ur, (14) 
(7) donne : 

U* = e—i0 p* (15) 


et l’on voit que le paramètre 6 est associé au 
nombre de nucléons : si le Lagrangien est invariant 
par la transformation de jauge 4 — eï 4, 
* — de, celui-ci est constant (nombre de 
nucléons diminué du nombre d’antinucléons, 
comme le montre (13)). 

Avec V* seulement, (13) s’écrit : 


pe=amp+bhn 


DE b, pt +an 
OU : 
ne = g; n° + bil— p°) (16) 
(pe = — 6 n° + ai (— po). (16) 


Si p, n correspondent aux valeurs propres + 1, 
De 
—Aders — Lo 
et p° à — 1. 
Remarquons enfin que la phase 0, de U, n’a pas 
de signification physique : elle est associée au 
nombre de protons et antineutrons, diminué du 
nombre d’antiprotons et neutrons. : 
On peut mettre en évidence les transformations 
simultanées V et V* d’une part, V, et VS d’autre 


de i) , n° correspond donc à + 1, 
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part, en disposant les 4 états possibles de parti- 
cules et antiparticules dans la matrice colonne : 


p 
210 

17 
Fe (17) 


— p° 


V, V* sont obtenues à l’aide des matrices à 4 lignes 
et 4 colonnes : 
Ti 0 
ie fe # 


V,, V* de la même manière, à l’aide des matrices 
à 4 lignes, 4 colonnes : 

ss É 0 ) 

HR TAN Er 


fo 1 20 ai 
AE DL RO E\CT (PM 


où L est la matrice unité à 2 lignes, 2 colonnes. 
Évidemment, les deux ensembles s’excluent 
mutuellement puisque le produit d’une transfor- 
mation V et d’une transformation V, permettrait 
de combiner net n°pet p',ce qui est exclu (schéma 


fig. 1). 
V'(r;) 


Cependant, on peut associer x3, diagonale, aux ;, 
ce qui donne : U = ef“ V dans la représentant (17) 
De même, on peut associer 7; aux x;, d’où : 


Us = eita0c (Fe 


B. Cas d’une paire de bosons isobares. — La 
paire X*, K°en est le seul exemple connu. 

Dans l’enceinte E considérée en A, on ne dis- 
tingue pas en quoi diffèrent ces deux états ortho- 
normaux et on peut encore les remplacer par Æ'*, 
K'® qui en sont déduits par une transformation 
unitaire U, : les relations d’anticommutation (8) 
sont seulement remplacées par des relations de 
commutation. 


+ 
Posant : Ë = (ro): on retrouve (11), ainsi 
; KT 
que (13) pour £° — (ai L 


Ici, la phase 0 de U ne peut pas être associée au 


N°4 


nombre de mésons À puisqu'ils n’y a pas conser- 
vations des Bosons. V unitaire unimodulaire définit 
toujours le spin isobarique. 

À côté de transformations (11) et (13), on peut 
ici effectuer des transformations : 


KT AR CRE (19) 
KT = B* K+ + A*K+ 
ete 
OCR? DK 
.. Fe (20) 
K9— D* K° + C* KO 


puisque, si XÀ* correspond au champ œ(w), K+ 
correspond à p*(— w), les relations de commu- 
tation s’écrivant : 


LA+, At] =[Xr, KF]= 0, (AFF, KE] = Fo SA) 


Il semble donc que les transformations V et W, 
soient simultanément possibles. En fait, par des 
produits VV, on obtiendrait des transforma- 
tions unitaires unimodulaires analogues à (19-20). 
Mais celles-ci ont une structure différente. 

En effet, pour que Æ”’+, K'+ satisfassent aux 
mêmes relations que À+, X+, À, B doivent satis- 
faire à : 


EEE (22) 


puisque : 
CAF, k'+]= |APTKX+, K+]+ |B[X+, K+] 
—(|4P — 18/) C4, At] 
tandis que : [K°+, K°+] = [K'+, K'+] = 0 résul- 
tent de (21) quels que soient À et B. 
De même avec (20) on obtient (1) : 
ICÉID E (28) 


Les transformations (19) satisfaisant à (22) 
dépendent en fait de 3 paramètres réels : 


A = ei% ch :. B = ei sh ne (22b) 
De même (23) donne : 
C = ei ch - D = ei sh : (236) 


Ainsi, si lon retient les transformations V entre 
K*, K9 (et V* entre KT, KÔ), il faut leur associer 


non pas V(V*) mais des transformations W 


(1) On admet évidemment que Æ°et Ksont discernables, 
c’est-à-dire que le champ Æ° n’est pas « absolument neutre » 
(hermitique). 


* 


TOMENIET À 


SENTE 


DV. 


Lila) ER 


be de 


at 2e, É 


a tt De nd de dt le nn à; Le Ré 
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définies par (19-20) avec les relations (22-23), sui- 
vant le schéma figure 2. 


Ho Va(c:;) K° 
W(xT) Wixiti) 
de V(rt:) (Ye 
A H° 


Toutes les transformations (19) et (20) peuvent 
être obtenues en associant une transformation (20) 
arbitraire W, à un ensemble (19-20) où les cons- 
tantes C, D sont liées à À, B, soit W;. Prenons par 
exemple : 

ASC ch 2)E PERD SnNUu 2) 
c’est-à-dire : 
RU NE 


MRC Dr =?) 


et comme au $ A, disposons les 4 états du méson K 
dans la matrice colonne : 


FH 
K° 
K5 


re 
Nous obtenons alors par (19-20) : 
RE 
K'9 
K°9 
ch{u/2) . sh?) K+ 
ch (u/2) ish (u/2) fe k° 
% . —ish (uf2) ch(uf2) | 1 Lo 
: sh (u/2) ch (1/2) — KT 


cette matrice étant : 
ch® IRC sh® .Xy To — 17/2 


où les +, x, à 4 lignes, 4 colonnes sont définies 
comme au $ A. Elle est analogue à celle qui effectue 
une transformation de Lorentz le long de Oy sur 
un spineur de Dirac : 


eYGU2 — € P102L/2 


nu 
K=A—(, sa PE 


D'autre part, les transformation #, V* s’obtien- 


en prenant 
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nent à l’aide des +, à 4 lignes, 4 colonnes comme 
au $ À, une transformation V (et V*) infinitésimale 
pouvant s’écrire : 


V=1 + (8m + PT: + ts) 


homomorphe d’une rotation spatiale sur un spineur 
de Dirac. 

En mettant de côté les transformations W, por- 
tant sur Æ9, X° seuls, on doit donc pouvoir écrire 
l’ensemble des transformations VW, comme homo- 
morphes des transformations de Lorentz homo- 
gènes sur un spineur de Dirac : 


PIE = : (ÔT, + Pts + œt:) 
d 
de OV T1 + We To + Vas Ta). 
En posant 9’ — 3/2 etc…., ceci correspond aux 


transformations infinitésimales à 6 paramètres 
réels : 


K+ = (1 +ix) At + (840) K° 
EVA +(N — in) (KT) 
ROSE PA 0) AT LE ix) A? 
STI) A ER) 
COURT = NV ART UN ip) A (24) 
| + (1 + io) Æ9 + (B + id) (— A+) 
TRE ST core 
+ (8 + 8) Kô + (1 — in’) (— KT) 
qu’il faut comparer à : 
ER GR Clos HS ane 
RE GS AR CEE PT Ce PT NC 


RARE * * * TT dE 
KT = a3 KT + as K° Lay KT + a K° 


RO CR DE I NRC AN LA 


où les 8 constantes complexes a,;, b, doivent vérifier 
les 6 conditions : - 


(264) 
la:l2 + la2l2= |a:l? + [a]? + 1 pour que [K+,K+]l=F 
(26b) 
[8,12 + 1h21? = ]b3l? + 1b4l? + 1 pour que [KS, x°l=r 


ai ba + c> D = a3 b3 + a b, (et C. C.) pour que : 
[Re Ke] =TA TA N|= 0 (6?) 
Gi bg +. Ge ba = @3 b1 + à b, (et C. C.) pour que : 
(A HO LAS, 020 NS) 


Pour une transformation infiniment voisine de 


l'identité, les relations (26) impliquent que 
a — À, b, — 1 soient imaginaires purs : 
a —=41+ix; by =1—1ie, 


Elles sont alors vérifiées au 22 ordre près. (27) 
donne, au 4er ordre : b, =— aÿ vérifiée dans (24); 
et (28) : 


— / e / 
b3 = 4 = V — 10. 


Les transformations (24) où €’ —«’, 0” —0 et 
b, = — af forment donc un groupe plus restreint 
que (25-28), ce qui est évident puisque les pre- 
mières dépendent de 6 paramètres réels, et les 
secondes de 16—6 — 10 paramètres réels. 

Cependant, si nous effectuons après (24), une 
transformation W, infinitésimale arbitraire (qui 
dépend de'3}paramètres réels) : 


ET RAT NOR 


KE0 = v* KO (1 — 10) K° 


nous remplaçons Æ”?, K'9 par : 
RD EC Lo) Rte Herald ue 
+ (X + zu + v) K° + v KH+ 
FAR Ge AO ee ed 
\ AE Te At (Pau ct 


c’est-à-dire que, dans W,5 VW, : 
Se a du et ol bi ATEN VU UP 


est sans relation avec a3. Enfin, effectuons 
après (24) une transformation e°’* infinitésimale, 
soit À + 10’ x. Nous obtenons : 


RP EG a et 00 Run ID RES AE OO 
DNA AN CRC ET em Ter 
RON NE ED 7 AR A a "Di Te 
(N + iu i 0’) K° + (v — à 07) Fes 


et les complexes conjugués, où 0’ 0. 

Les transformations (25-28) les plus générales 
sont donc bien obtenues comme produits des VW, 
homomorphes des rotations dans un espace à 4 dimen- 
sions à métrique de Minkowski (à 6 paramètres), 
par ef°°% (1 paramètre) et par les W,sur X°, K0 seuls 
(3 paramètres). 

ei9"%s VW, est l’analogue de ex} du $ A. 

+3 est l’analogue de y, pour les spineurs de Dirac 
dans la représentation Ys = p,, Ya = p3. On peut 
donc dire qu’elle effectue une « symétrie spatiale » 
dans cet espace des états de Minkowski : comme 4, 
x, cCommute avec les «rotations spatiales », c’est-à- 
dire avec les r;, mais non avec les «rotations » dans 
un plan défini par un «axe d’espace » et un «axe de 
temps », c’est-à-dire avec les x, +. 


+ K0 
Les deux doublets (Æ | et _ sont les 
re oh? 


analogues de deux spineurs de Pauli de parités 
opposées. En revenant à la signification physique 
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de +3, #3 on a le tableau des valeurs propres sui- 
vantes : 


T3 X3 
p + 1 + 1 
n — 1 + 1 
ne + 1 — 
p° — 1 — 1 
1} Qu + 1 + 1 
K9 — 1 + 1 
K° ESA M 
K+ 21 —4 


Si l’on considère l’ensemble des baryons, on sait 
que l’ensemble des valeurs propres de T3 = T3/2, 
T (ici T — 1/2), x3 (particule ou anti-particule) ne 
suffit pas à les classer. 

En se limitant aux doublets d’isobares (7 = 1/2), 
il faut une autre « variable dichotomique » pour dis- 


tinguer le doublet (? 


4 (avec ses antiparticules) et 


Can 


le doublet (es) (avec ses antiparticules). 


Ce nouveau nombre quantique est lié à l’étran- 
geté S qui indique aussi le décalage des valeurs de 
la charge pour un même 7, : 


Q—=T:+ + D = gro + % À + 6) avec À = 1. 
| (22) 
S —0 pour les nucléons, S$ —=—2 pour les Æ 


(et + 2 pour E). 

Pour les mésons X, on ne connaît qu’une paire 
d’isobares et les valeurs propres de +3, x3 suff- 
raient à distinguer les différents états de particules 
et antiparticules. Cependant, comme À —0, la 
relation (29) conduit à leur attribuer une valeur S 
(+ 1 pour les X, —1 pour les Æ). Celle-ci n’a pas de 
relation avec »x3. D'ailleurs, la valeur de x, est 
purement conventionnelle pour les K puisque le 
nombre de ces particules n’est pas conservé : on 
pourrait échanger les dénominations Æ° et ÆX9 
(ce qui n’est pas possible pour », n° une fois p, 
p° fixés). 

Au contraire, les valeurs de S (ou de +3) peuvent 
être déterminées par les règles de sélection dans les 
processus de création associée hypéron + X. 

Ainsi, aux groupes de transformations que nous 
avons étudiés, se superpose toujours un autre 
groupe définissant les valeurs de S. 


_C. États du neutrino. — 10 Le champ neutri- 
nique vérifie les équations de Dirac (10) où 
Mo — 0. On sait que les combinaisons : 


L + y, (l 
(Is = D) L W, da = % Fe ” 
orthogonales, vérifient alors chacune (10) et sont 


#, 


À 
ps 
: 
fi 


PSE F5. TNT — 
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indépendantes ; elles ne sont couplées ni par les 
équations d'évolution (10), ni par les relations 
d’anticommutation : dr, Qu | = | VX, ds | = 
comme on le voit en prenant une représentation 
où y; est diagonale : 


|: )= : . = (| à 
LS NS EE RE À o)° 


En effet, dans cette représentation : 


UE et Da = 
Ÿ9 ( Ya de 


0 
(Ÿ%,, Ÿ, étant des spineurs à 2 composantes) corres- 
pondent à des indices différents. 
On se trouve apparemment dans un cas ana- 
logue à celui du $ À pour une paire d’isobares, 
avec les états : 


| (in os. (31) 


Avec la représentation (30), ces 4 états sont 
représentés par les 4 spineurs de Pauli : 


(a 


D’après (31) Ÿ£ est, comme d,,un état propre 
de y; avec la valeur + 1, 4£ avec la valeur — 1, 
comme , (?). 

20 Comme au $ A, il semble que l’on puisse 
effectuer les transformations (11-13) sur les dou- 


blets ” 
Rae a (0) 


c’est-à-dire en particulier les transformations V 


sur 
V9 “ à 
(5) “ és ds) 
Cependant, les équations d'évolution s’écrivent 
avec la représentation (30) : 


Ro pop 0. (ut) 


db, et 4, représentent des états symétriques l’un 
de l’autre. Dans un système d’axes droit, 4, (ou 4) 
représente un état gauche d’hélicité — 1, 4, (ou à) 
un état droit d’hélicité + 1. 


(2) On peut dire aussi que 4, Ÿ5 ont la « chiralité » + 1, 
Va, VS la chiralité — 1. 
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Une transformation V qui combine ces deux 
états ne conserve donc pas strictement les équa- 
tions d'évolution et doit être exclue (tout comme (9) 
au $ A (5)). 


On peut aussi, comme au $ A, effectuer sur les 


doublets 
(D « (à 


les transformations (12-14), en particulier les V.. 
Ici Yet d£ d’une part, déet 4, de l’autre obéissent 
bien aux mêmes équations d'évolution (10 b). 

Nous avons alors les V, associées aux x, de $ A 
et eirûe — bide +, correspondant ici à Ys. 

Nous pouvons associer ici V, aux matrices 7, en 
inversant les rôles des 7; et x;, c’est-à-dire en écri- 
vant les 4 états du neutrino et de l’anti-neutrino 
dans l’ordre : 


REA (32) 


Les valeurs propres de 3, x, sont alors associées 
aux différents états suivant le tableau : 


va OI AG 
Vo A dl 
Vg — 1 2e Al 
La ul =" 
28 A1 1 


x, Correspondant maintenant à y;, T, au Caractère 
neutrino, ou anti-neutrino. 

Ainsi, le spin isobarique T = 1/2 se trouve 
attribué à la paire 


C) 


Les deux états de la paire diffèrent par des pro- 
priétés «intrinsèques », comme deux états isobares, 
l’un étant un état neutrino, l’autre un état anti- 
neutrino de même hélicité. L’analogie obtenue 
pourrait donc avoir une signification physique. 
C’est ce que Heisenberg et Pauli ont admis dans 
leur tentative de théorie unitaire, en suivant la 
remarque de Gürsey. 

Remarquons que l’on peut faire porter la trans- 


formation : 
> a b 
Poe 4 : 


sur les champs : 


y correspondant à T3 = + 1. 
ve correspondant à T3 = — 1. 


(3) Ceci n'empêche pas que la différence entre les états 
isobares puisse éventuellement résulter d’une sorte d’«asy- 
métrie interne » ne correspondant à aucune asymétrie Spa- 
tiale ordinaire. 
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dans leur ensemble, en combinant les relations (x), 


(y) et (8), (à) de : 


| Yo — a Vo + b Vo (x) 
À Se = — 0 ve + at de (f) 
base () 
Va = — b* Va + a* da (È) 
où l’on a posé d, = dj et da = — Y$ pour mettre 
en évidence les hélicités de chacune des paires : 
Ÿ Va 
() « (à 
On obtient ainsi la transformation (9) de 


Pauli [1] : 
l=aÿ+bysd, de a* 4e—8* y à, 
avec : [a[? + |b[? 
qui est compatible avec les équations d'évolution 
(celles-c1 étant strictement conservées) dans le cas 
actuel mo —0. On doit y ajouter les transfor- 
mations à 1 paramètre : 
= eve 4, 


introduisant le déphasage 20, entre d, et d,. 
Gursey [2] a précisé l’analogie en posant : 


d'e == eye be 


A + y; de 
Pi=dXo ia = EE 
(34) 
1 à L— 
Eee nr x — — te 


En considérant aussi les conjugués de charge 
de (31), on en tire : 


9 4 + 1 — y; 

Rene do 
Le Learn 
yen ae 


Ces champs y, € définis par Gürsey satisfont à 
des équations couplées, invariantes par (9°) sur y, C 
ce qui correspond bien à une transformation V 
pour p, n, d’après (30) et (31). Elles sont aussi 
invariantes par 

Anne r Men ED AUS DU a der CR 
d’où la transformation de jauge en 0,. 

Evidemment (33) implique l'existence dans la 
nature des 4 états d,, V6, 5, d, du neutrino. Dans 
une théorie à 2 étais du neutrino, on n’a que 2 états 
d’hélicités opposées qui ne peuvent être combinés 
par une transformation V, : 4, et 45 dans la théorie 
de Weyl et Lee-Yang, à, et 4, dans la théorie de 
Majorana. | 

Si le spin isobarique Ta un sens, dans ce 
dernier cas 1l devrait être pris nul, comme pour A° 
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(états A°et A). Il en serait de même pour tous les 
leptons autres que le neutrino qui n’ont qu’un état 
de particule (et l’état antiparticule correspondant), 

Âvec la première hypothèse, les désinté- 
grations É : 


n—>p+eT +vw; pn+et+v 


1 4 1 1 1er 1 
a Vip = re RE EN || = 
ee 2° T2 ù 2 RES 2 


> 

vérifieraient la règle AT] — 17/2 des interactions 
faibles (en admettant que 7 = 0 ait un sens pour 
l’électron). 

Avec la seconde (7 — 0 pour v, v°) cette règle ne 
serait pas vérifiée (4). 

D'autre part, l'incertitude sur l’état (4, ou dé) 
d’un neutrino d’hélicité déterminée ne peut pas 
être levée aussi facilement que pour l’état de 
charge du nucléon, par l’application d’un champ 
électromagnétique. Ici, un processus de désinté- 
gration $ inverse pourrait seul, en principe, lever 
cette incertitude, et fait disparaître le neutrino. 


Conelusion. —— 1° Dans le cadre de la théorie des 


champs actuelle, l’analogie suggérée par Gürsey 
entre le doublet proton-neutron et le doublet 
neutrino-antineutrino à quatre états, conduirait à 
attribuer le spin isobarique 7 — 1/2 à ce dernier. 
Pour un neutrino à deux états au contraire, comme 
pour les autres leptons, on serait conduit à prendre 
T = 0, si la notion de spin isobarique a encore un 
sens. 

20 À cause des relations de commutation de 
champs différentes, la dégénérescence des états 
isobares en l’absence de champ électromagnétique 
est d’ordre plus élevé pour les Bosons de T = 1/2 
(mésons À) : elle porte sur les 4 états de particules 
et antiparticules. Parmi les transformations alors 
possibles entre ces 4 états, certaines sont d’un type 
essentiellement différent de celles que lon a pour les 
états isobares de Fermions, ou les états du neu- 
trimo : le groupe de toutes les transformations 
contient un sous-groupe homomorphe des rota- 
tions dans un espace quadri-dimensionnel à 
métrique de Minkowski. 


Nous avons une grande dette de reconnaissance à 
l'égard du Pr Frédéric Joliot-Curie pour les encou- 
ragements qu’il nous a toujours prodigués, pendant 
plus de dix ans, au cours de nos travaux. 


(*} Quelle que soit l'hypothèse, elle ne serait pas vérifiée 
pour certaines désintégrations. Ainsi : WE — e£ + y + ve 
n’y satisfait pas ; T pour y + ve étant de toute façon 
entier avec TE") 


Manuscrit reçu le 31 octobre 1958. 
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PRÉVISION DE LA SUSCEPTIBILITÉ PARAMAGNÉTIQUE DES CORPS FERROMAGNÉTIQUES. 


APPLICATION AU COBALT. 


Par D. TAUPIN et G. FOURNET, 


École Supérieure de Physique et de Chimie de la Ville de Paris. 


Résumé. —— Au moyen de la théorie des phénomènes coopératifs d’Yvon, nous avons obtenu 
l'expression de la susceptibilité paramagnétique des corps ferromagnétiques. Dans le cadre du 
modèle d’Ising, cette expression nous a permis de voir l'influence des termes à deux nœuds, trois 
nœuds, quatre nœuds, etc. Nous pouvons également prendre en compte, non seulement les énergies 
d'interaction entre premiers voisins, mais aussi entre deuxièmes, troisièmes voisins, etc. Pour le 
co nous avons obtenu un très bon accord entre notre courbe théorique et les résultats expéri- 
mentaux. 


Abstract. — By means of the Yvon cooperative phenomena theory, we have obtained the para- 
magnetic susceptibility of ferromagnetic materials. If we consider the Ising model we may study 
the influence of the terms of two sites, three sites, four sites, etc. We may also introduce the 
Our theoretical 
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interaction energy between first, second, third, etc, nearest neighbours. 


curve y(7') for cobalt is in quite good agreement with observation. 


1. Introduction. — De nombreuses études ont 
été consacrées à la statistique des corps ferro- 
magnétiques ; outre les travaux de Néel [1] consi- 
dérant les fluctuations du champ moléculaire, la 
plupart des auteurs utilisent le modèle simplifié 
d’Ising [2] où les interactions entre atomes sont 
supposées se limiter aux premiers voisins. On cons- 
tate alors que les prévisions coïncident assez mal 
avec l'expérience, surtout au voisinage de la tempé- 
rature de Curie. Toutes les méthodes de calcul 
employées sont soit des méthodes approchées — 
telles celles de Bethe [3], Kikuchi [4], Zernike [5], 
Ter Haar [6] — soit exactes, mais donnant toujours 
le résultat sous forme d’un développement en série 
(Trefftz [7]). On pourrait imputer les divergences 
entre les prévisions et l'expérience à l’imprécision 
des méthodes de calcul, mais il est beaucoup plus 
probable que les défauts sont dus au modèle utilisé, 
celui d’Ising. 

Nous nous sommes proposés de faire une nouvelle 
comparaison entre la statistique du modèle d’Ising 
et l'expérience en nous intéressant à la suscepti- 
_bilité paramagnétique. 

Nous avons utilisé une méthode de calculexacte 
qui, comme celle de Trefftz, fournit les résultats 
sous forme de série rapidement convergente : la 
théorie d’Yvon [8] et [9]. Yvon avait lui-même 
amorcé le calcul de la susceptibilité, mais il n’avait 
explicité que les deux premiers termes du dévelop- 
pement. Fournet [10] a repris ces calculs de manière 
plus approfondie, afin de calculer les températures 
de Curie dans les divers réseaux simples. Nous 
avons utilisé ses résultats pour le calcul de la 
susceptibilité dans les réseaux cubiques centrés 
(Fer) et à faces centrées (Cobalt et Nickel). 

Nous avons ensuite perfectionné la théorie pour 


tenir compte des interactions avec les voisins 
d’ordre supérieur à 1, et l’avons appliquée au cas 
du cobalt. 


2. Rappel de la théorie d’Yvon. —- Rappelons 
brièvement les hypothèses de la théorie d’ Yvon : 

— Le réseau est supposé indéformable. 

— Le moment magnétique de chaque atome peut 
prendre 2 positions, l’une parallèle, l’autre anti- 
parallèle au champ magnétique Hi. 

— L'énergie mutuelle de deux atomes ne dépend 
que de l'orientation relative de leurs moments 
magnétiques et de leur distance. 

— On suppose que les seules énergies qui inter- 
viennent sont les énergies potentielles dues au 
champ macroscopique et les énergies d'interaction 
entre les atomes deux à deux. 

— Il est fait usage de la statistique classique. 

Dans sa théorie générale sur la thermodynamique 
ordre-désordre des cristaux constitués de deux 
sortes d’atomes À et B, Yvon obtient une équation 
de l’ordre à grande distance sous la forme 


Nr T3) Ci = RTS (1) 
dans laquelle 
n; — probabilité d’avoir en i un atome de type A, (2) 
Wis = Waais + Wei — 2Wapi; (3) 


où Waax est l'énergie mutuelle d'interaction de 
deux atomes de type À situés respectivement aux 
nœuds i et J, 

Ui = Uni — Uni (4) 
où UA 1 est l'énergie potentielle d’un atome A situé 
eni. 

Yvon décompose le premier membre de (1) en 
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termes de complexité croissante classés d’après le 
nombre de nœuds distincts qu’ils font intervenir 


f(ni, es W;5) 


log Ass + Dog Ai + 2 log Ask + .… (5) 
E ik 
Le fait remarquable est que les Avi... Sont des 
fonctions de z, W;; et T indépendantes du réseau, 
lequel n'intervient que dans leur existence et leur 
dénombrement. Ils sont calculables 4 priori dans 
un réseau fictif à nombre fini de nœuds. 


3. Application aux corps ferromagnétiques. 
Caleul de la susceptibilité x. — L’équation (1) 
devient, à condition de considérer qu’un atome À 
est, par exemple, un atome de moment + uet B, 
un atome de moment — pi : 


HR EE CREER PETITE (6) 


où À; est le champ agissant sur l’atome situé au 
nœud 1. 

Nous savons qu’au-dessus de la température de 
Curie et en l’absence de champ magnétique l’ordre à 
grande distance a complètement disparu et par 
conséquent 

ni = 1/2 
ce qui se traduit (cf. (6)) 
HOT Mi) = C 

Pour les mêmes conditions de température l’ac- 
tion du champ magnétique ne fait que varier très 
faiblement la probabilité r; et nous pouvons donc 
écrire 


DA ARRETE Î 
ff 7 mi) ie 


+= cet (8) 


Tee (7) 


en ne conservant que les termes du premier ordre 
dans le développement de la fonction f. La compa- 
raison de (7) et (8) fournit 


NU ou 
SL ( T, is): (ru —;) - PE (9) 


Il est commode de poser 
m=>{t+s) (10) 


et d'utiliser les quantités 


1 dAn ij 
He (11) 


déjà introduites dans un précédent article [10] pour 
obtenir expression de la susceptibilité 


2u? N 
X = (12) 
RT (ai + D ai; Don sin He) 


N° 4 


où W est le nombre de nœuds par cm, gramme, ou 
atome-gramme, selon qu’on considère la suscep- 
tibilité par cm, spécifique ou atomique. 


4. Développement de l’expression de x valable 
dans le cadre du modèle d’Ising. — Nous suppo- 
serons que les interactions se limitent aux premiers 
VOISINS : 


pour des premiers voisins : 


Wa = W 4 
DE = pour les autres. . 
Pour simplifier l'écriture nous poserons 
x = eWil2RT, (14) 


Fournet a calculé et dénombré les termes de la 
parenthèse de l’équation (12) et nous nous bor- 
nerons à rappeler ses résultats : 


RÉSEAU CUBIQUE CENTRÉ : 
1 EL 8 as LES aus 04 a Sa RE 


(La figure 1 indique les interactions représentées 
par ces termes.) 


À 2 À 3 A 41 À 42 
A 43 A4 A 5,4 
Ass À 61 À 62 

À6s 
Eric de 


RÉSEAU CUBIQUE A FACES CENTRÉES : 


41 + 19 2 + 24 az + 180 4, ; 


+ 120 das + 120 Gus + 24 Gus (16) 


N°4 


Au lieu de considérer, comme Fournet le faisait, 
un seul terme à 4 nœuds, pour chaque figure de 
4 atomes, nous avons préféré scinder ce terme en 


plusieurs de complexité croissante. Par exemple 
nous posons 


da,2 — y 


Gas = Ga, + ae. 


Ces nouveaux paramètres sont intéressants pour 
deux raisons : 

Le terme ü@y,2 par exemple diffère peu du 
terme &4,, et dans une première approximation 
leur différence 44, pourra être négligée ; nous 
aurons donc moins de termes à calculer. 

Les formes mathématiques de nos à, sont plus 
simples que celles relatives aux 4, ; nous avons en 
effet défini les & de telle façon qu’ils n'existent 
que si toutes les liaisons elles-mêmes existent et 
par suite l'annulation de chacune d’elles isolément 
correspond au moins à une racine simple æ = 1. 
Les &, contiennent donc toujours un fac- 
teur (x — 1)", n étant au moins égal au nombre de 
liaisons que le a, fait intervenir. 

Nous avons posé de même : 


aa — Aaa — Ua,1 (17c) 
Gas — Aa,s 3 G4,a 3 Ga,2 — 3 Gy,1 (174) 
dos — gs (17e) 
Gi = A1 3 (17) 
Ugo = go — 6,3. (178) 


Nous donnons ci-dessous les expressions des a 
déjà calculées par Fournet [9] ainsi que les nou- 
velles pour les à : 


CEE | (18a) 
= 
HUE _—— : à (188) 
7. 
Bai = Gas1 = Se . (18c) 
> (22 \ 
Ug,2 = 5 RAT 
DU + 3) @ +1) Fe 
TT GE ni) 
a Dr El} (&r 3) 
E (DA 2 3x x 5) 
Va DE re 3 
Dire A) (re 8) 8e) 
(os re 2re 1) 
da,a — (x? + 1) (a? + 3) | 
n 3{x —1)6 (x$ — x? — 2) 
U4,3 — € TS) 2 2 
2(x? 2 le 3) (a NE ) (18) 
3(x — 1)4 (3x — 1) 
Aa,3 Da? — x Eu 1) (x? 2e 3) 
e ed). 
DIRE 
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Nous corrigerons ici une expression donnée dans 
l’article déjà cité [10] 


D) 5 


6,3 le,3 — (2x? + 1) &? +3) (18h) 

EI 9f2 + 1) (x? + 3) (3æ° + 1) 

; z—1)S (x +1)? 

D NOR ET (re CENT AC ET) 

de AN (at D Sale re 2) 

fe. (2 Dre 3) (5x? + 1)? 8x) 
(m —1)6 (78 + 2 67 Er 

ere } 


(PCT TES) NOTE En) E 


Nous avons calculé numériquement ces diverses 
fonctions en fonction d’une variable réduite m, pro- 
portionnelle à la température ; nous avons posé : 


m = —2kTIW, (19) 


== CON (20) 
Pour les corps ferromagnétiques W,; est négatif 
et la variable m est donc positive. 
Pour fixer les idées sur les ordres de grandeurs, 
nous préciserons que, selon les approximations 
cf. [1ÔJ), pour T = T,, température de Curie, 


Y= m[2+8as+...] 
8 


RS US r 6.585 8e 
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m, est voisin de 3,5 pour le réseau cubique centré 
et de 5,25 pour le réseau cubique à faces centrées. 

Il est commode d’étudier, non la susceptibilité %, 
mais son inverse 


Des 2 7, man DD dit Das Li) NeNT) 


1 /x est proportionnel à une fonction réduite Y, 
laquelle ne dépend plus que du réseau et peut être 
calculée numériquement & priori : 


T= more LOUE UE 


(22) 


Ve mfo;+2a,+...] 


Y= m[2+#2a+2403+180 à,,] 


Y= le lee 110042 a 
IE (2 


l /) 
3) À) 


Séut24Q4s] 


Y2mf2#2 a+ 24 


7 


F1G. 3. — Réseau cubique à faces centrées. 


Nous avons indiqué (fig. 2 et 3) les courbes pour 
les diverses approximations suivantes : 


RÉSEAU CUBIQUE CENTRÉ : 


approximation à 2 termes (série limitée 
termes à 2 nœuds), 

approximation à 3 termes (série limitée 
termes à 4 nœuds), 

approximation à 4 termes (série limitée 
termes à 5 nœuds). 


aux 
aux 


aux 


RÉSEAU CUBIQUE A FACES CENTRÉES : 


approximation à 2 termes (série limitée 
termes à 2 nœuds), 


aux 
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approximation à 3 termes (série limitée aux 
termes à 3 nœuds), 

approximation à 4 termes (série limitée aux 
termes à 4 nœuds), les &y,9 dy, da,a étant négligés), 

approximation à 7 termes (série limitée aux 
termes à 4 nœuds, les Gyo, Qu,3) Gasa n'étant plus 
négligés). 


5. Examen des courbes. — Il apparaît que le fait 
de négliger les Gy,9, 4,3 Et Ga,a à peu d'importance, 
d’ailleurs, parmi eux, seul 44,4 à une valeur notable. 

On peut constater que la convergence est beau- 
coup plus rapide pour le réseau cubique centré. 

Un autre fait remarquable est la’faible courbure 
des courbes Y(m") ; en particulier, et contrairement 
à l’expérience [11], on ne constate pas d’accrois- 
sement notable de cette courbure au voisinage du 
point de Curie. (La courbe en pointillé sur la 
figure 3 représente les résultats expérimentaux 
relatifs au cobalt.) 

Enfin, la position de l’asymptote, très écartée de 
la partie accessible à l’expérience, semble rendre 
illusoire toute détermination expérimentale [11] 
des températures paramagnétiques de Curie. Pour 
améliorer la coïncidence entre les courbes expéri- 
mentales et théoriques, on pourrait évidemment 
pousser plus loin les développements en série, mais 
il est bien clair que cela ne changera pas l’aspect 
général des courbes, et en particulier la courbure au 
voisinage du point de Curie. Il faut donc perfec- 
tionner le modèle d’Ising. Nous avons donc fait des 
essais en considérant les interactions non seulement 
entre premiers voisins mais aussi entre voisins 
d’ordre supérieur. 


6. Caleul de la susceptibilité du cobalt, en tenant 
compte des interactions avec les voisins d’ordre 
supérieur à 1. — Lorsque nous supposons que les 
interactions ne se limitent plus aux premiers voi- 
sins, d'énergie mutuelle W;, mais ont lieu éga- 
lement entre voisins plus éloignés, d’énergies 
mutuelles W,, W3 etc, l’équation (12) devient : 


x = QU NIRTa; + p1 ao(Wi) + po Go(Wa) + ps 41(Wa) 
onde dr an Je, 7) 

+ ris As(Wa, Wa, We) + r199 Ga(Wi Wo W3s) “578 

(23) 


Les p, r ... sont des coefficients entiers qui 
dépendent du réseau considéré. 

Etant donné que les interactions entre voisins 
d’ordre supérieur à 1 sont faibles, nous nous bor- 
nerons à les faire intervenir dans les termes 
en dyWo), aW3) ... et nous continuerons à 
négliger les a;(W, W, W.), as(W,, W:, W,) 
Les calculs précédents nous ont également montré 
que les Gy,2, Ayg Axa pouvaient raisonnablement 
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être négligés ; notre équation (23) devient donc : 


= UNE NIRAT[a + 12a:(W;) + 24a3(Wa) + 180aa,1(W 1) 
+ 6 a(W, JP=P28 aolWa) + ... Pn Ga(Wa) + .] (24) 


Nous ferons encore une approximation : D’après 
(14) et (18a) : 


ao(Wn) = eWnl2RT 1, (25) 
Nous admettrons, et nous vérifierons à poste- 
| riori Que pour n > à W,lkT est petit, et que par 
suite 

Wa 


(26) 
Il vient alors 


KT[a1 + 12a2{W;) + 24 a3(W;) + 180 aa 1lW;)] 


FU 


1, [ee] 
T3 à Pn Wa 


Nous avons déterminé les constantes u, W, 
et E Pn W, de façon à ce que la courbe théorique 


Gi. %) coïncide au mieux avec la courbe expéri- 
mentale pour Co [11]. 

Nous avons trouvé : 
spécifique) 


1 [x = 28 850 + 40 T[2 + 12 a:( Wa) + 


(il s’agit de susceptibilité 


“icosaure-tmne 


(28) 
Soit 
uw = 1,40 magnéton de Bobr (29) 
; W; = —1,068.10—18 erg — —— 0,0667 eV . (30) 
4 CO 
2, 2, Pn Wn = — 0,995 1018 erg — + 0,0624 eV. (31) 


Rappelons que, pour deux nœuds en positions 
de n° voisins : 


Wan = Waais + Waris — 2 Wapii. (32) 

Remarquons que p, est au moins égal à p, — 6 
ce qui fait que les W,/2 (tous positifs) sont infé- 
rieurs à 0,01 eV. Par suite W,/24T est de l’ordre 
de 0,08 et ‘notre approximation (26) était justifiée. 


PRÉVISION DE LA SUSCEPTIBILITÉ PARAMAGNÉTIQUE 481 


L’examen des courbes (fig. 4) montre que, dans 
la limite des erreurs Fan AE la fonction BE) 


15000 TN SR a = Sete 
k spécifique ee 


14000 
15000 


PASS ee Sen FT 


FOOD ONE re 


10000 


9000 | = ———+—— 


8000 


7000 |———— 


Se = 


EXPÉRIMENTAL à) 
THÉORIQUE 


T (°H) 
1650 700 


1450 1500 1550 
Free 


1350 1400 1600 


représente de manière très satisfaisante les varia- 
tions de la susceptibilité en fonction de la tempé- 
rature. 


7. Conclusion. — La présente théorie ne peut 
avoir la prétention d’être une théorie exacte, 
puisqu'elle fait intervenir séparément le cortège 
électronique de chaque atome ; une théorie exacte 
devrait être au moins en partie collective et déduite 
de la théorie des bandes. Nous croyons cependant 
que notre travail permet d’obtenirune bonne appro- 
ximation des faits expérimentaux. 

Nous sommes heureux de remercier M. Herpin 
pour nous avoir donné l’idée initiale d’effectuer ce 


travail. 
Manuscrit reçu le 12 décembre 1958. 
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COUCHES NUCLÉAIRES ET ÉNERGIES DE DÉSINTÉGRATION DES ÉMETTEURS 


Par Maurice-E. NAHMIAS (*), 


Laboratoire de Physique et Chimie Nucléaires du Collège de France. 


Résumé. -- Le premier article publié sur ce sujet dans le Journal de Physique et le Radium en 
avril 1958 a été remis à jour au moyen de la nouvelle Table des Isotopes, parue en avril 1958, dans 
la Review of Modern Physics. Il subsiste toujours des anomalies qui ne sont pas communes aux 
neutrons et aux protons dans le remplissage des couches nucléaires. 


Abstract. — The plots of Qg versus the number of protons or neutrons published in the first 
paper on this subject, in Journal de Physique et le Radium of April 1958, have been revised with 
the help of the new Table of Isotopes — Review of Modern Physics, April 1958 —. There are still 
discontinuities in the two kinds of plots which do not coincide. 


Protons (I) 


La DR 
| Î mc? Eee 


RÉCRSREE 
1 
ARE EE 


Fc. TL A XIII. 


En abscisses : 


+ D 1] + A a Ua ere > x a] d ‘ j ( l fi ñ 1 
EU : Pour les émetteurs 8 le nombre N de neutrons du noyau émetteur est donné en abscisses. Pour les émetteurs 
Ê : Les sen ue des eus zéro, les abscisses donnent le nombre N de neutrons du noyau final 
aut donc retrancher un neutron à lindication de l’abscisse pour avoir le nombr 6 ins- 
û e de neutror - 
crit sur les courbes. D 


Protons : Pour les émetteurs B— l : 6 
“ DR teu Le les abscisses donnent le nombre de protons Z du noyau émetteur. Pour les émetteurs 
. On sur l’abscisse le Z du noyau final. Il faut donc ajouter un proton à l’indication de l’abscisse pour avoir 
le Z du noyau émetteur inscrit sur les courbes. 
En ordonnées : 


Les énergies totales Æ de désintégration B + Y des radio-éléments B sont représentées en unités MeV. en ordonnées dont 
la partie supérieure est relative aux radio-éléments f— et la partie inférieure aux GA 


(*) Adresse actuelle : $, H. A. P. E. Technical Center. La Haye. Hollande. 


No 


Le premier article sur ce sujet a été envoyé au 
Journal de Physique et le Radium en décembre 1957 
et fut publié [1] en avril 1958. Les derniers docu- 
ments utilisés dataient de fin 1956. Une nouvelle 
Table des Isotopes est parue [2] dans le numéro 
d'avril 1958 de la Review of Modern Physics et m’a 
servi à remanier complètement les graphiques du 
premier article. 
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Voici ce qu’on peut déduire des 43 figures qui 
suivent et pour l'établissement desquelles on 
voudra bien se référer au premier article cité. 

Sur 24 anomalies, 10 sont communes aux neu- 
trons et aux protons ; 8 sont uniquement visibles 
sur les graphiques relatifs aux neutrons, dont 4 
pour un nombre de neutrons supérieur à 100 et par 
conséquent non contrôlables sur des graphiques en 


Protons (II) 


Neutrons (IV) 


protons. Enfin 7 discontinuités sont propres uni- 
quement aux protons. 

Le tableau I résume ces constatations. Lorsque 
les anomalies sont voisines, comme à 24 et 26, à 50 
et 52, ou à 58 + 2 nucléons, protons ou neutrons, 
j'ai désigné la discontinuité par la même lettre : 
a et a’, Cet C’ ou f, f’, f”, la majuscule faisant 
ressortir un nombre magique du modèle nucléaire à 
couplage spin-orbite. 

Lorsqu'une discontinuité est fondée sur un seul 
point expérimental la droite à été hachurée et la 
position de l’anomalie est dans ce cas douteuse, 
exemple f” à 60 neutrons sur la droite 19-17 de la 


figure VIII, bien que cette discontinuité soit con- 
firmée sur la même figure sur la droite 17-15. 

Ces graphiques illustrent, une fois de plus, lexis- 
tence de nombres magiques correspondants à des 
couches saturées avec 14, 20, 50, 82 et 126 parti- 
cules et des sous-couches à 38, 58 et 68 nucléons 
Ils montrent aussi qu’il existerait d’autres confi- 
gurations particulièrement stables à 24, 26, 30 et 
56 nucléons (protons ou neutrons), 36, 42, 62, 72, 
84 et 88 protons, 52, 60 et 74 ainsi que 120, 122 et 
134 neutrons. Certaines de ces configurations sont 
proches, à une paire de nucléons près, de sous- 
couches prévues par le modèle de Mayer-Jensen. 
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L'absence de tout indice d’une couche à 28 nu- 
cléons est pourtant assez énigmatique. On observe 
une pénurie de radioéléments dans cette région — 
voir droite 5-3, figure II —. Le Ti49 avec 27 neu- 
trons est stable — 5 % — ainsi que le Cr53 avec 


N°4 


29 neutrons —9,5 %— Par contre le #Ni56 dou- 
blement magique n’est pas stable pas plus que 
le ?7Co55, ni le $Mn53, ni le #Sc49, alors que sont 
stables le %Fe54(6 %), le ZÆCr52(84 %), le 
28 V51,(99,8 %) et le ATI50(5 %). Enfin le ÆCa48 


Protons (W) 


| 
ml147 IL Eu 155 

U x: 
C2 Gd” (27-25) 


at iot 


Neutrons (V1) 


n’est pas particulièrement stable puisque sa con- 
centration isotopique n’est que de 0,18 % et qu’il 
est radioactif $, avec, il est vrai une période 
de 108 années. Pour le nombre magique 20 on 
observe au contraire avec les neutrons la grande 
stabilité de #K39(93 %) et de 170137(24,5 9). Der- 
nière remarque, parmi les 10 isotopes connus à 
28 protons du nickel, 5 sont stables, dont 4 pair- 
pair et 5 instables, dont 4 impair-pair. 
Remarquons en outre que ces graphiques 


peuvent si bien servir à interpoler ou à extrapoler 
des valeurs de Qs que plusieurs points anormaux 
des figures du premier article [1] se placent à 
présent beaucoup mieux sur les droites, avec les 
nouvelles valeurs publiées. C’est ainsi qu’on peut 
par exemple déduire de la droite 19-17 de la 
figure VIII que ©» du Pd111 doit être proche 
de 2,6 MeV alors que les valeurs expérimentales 
données jusqu'ici sont : 2,15 par spectromètre [3] 
et 3,5 par absorption [41. 


r 
N°4 COUCHES NUCLÉAIRES ET ÉNERGIES DE DÉSINTÉGRATION 485 


îh 


; Neutrons (IX) f$nis 
| Le 


ce 


rni par Au191 dont Q droite 43-M de la figure VII bis — ce qui confirme 
he de 1,9 MeV, voir la coïncidence f-Y de 4,6 + 1,8 MeV jusqu'ici 
considérée comme douteuse [5]. 


Un autre exemple est fou 
(C. E., B*) doit être proc 
droite 35-33 de la figure XI. 

Enfin le @e de 1r197 se place bien à 3,4 MeV pu Manuscrit reçu le 3 novembre 1 958. 
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ce 


(55-58) 


Protons (XIII) 
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Neutrons (XII) 


Pa 22 (H7:4S) | 
(45-43) 1 ll 
m 
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Ê Ô ° At !211 
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5 (87-35) —TTLISS l (43-41) 
Hp + +697 
35-33) | pres 
ie | Bi 203 
TT Ie 
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TABLEAU I 


Coucnr Ne ANOMALIE EN NET P EN N NES EN P 
INDex Ale SUR FIGURES : SUR FIGURES : 
DOI E SE eee Se A SR A EU AR Re LATE 
RÉNES TRS 20 RAA TA ner tac Len 
24 es ne er tb re ne one RTL 
26 RU See On TO NOT DATI 
28 :11f7/2 
30 VAE SE NP REA ee À Je LHLRIONE 
34 : 1f5/2 
36 CRT ee ren Se tout PR PROS ie me ME VU Rio IIL 
DOME RER de 38 L'TE OmR ed IT, IV 
40 : 2p1/2 
42 CRT Ie DD me nie del dede La pal dot er ne TT PE INPALVE 
ROSES LE MP SAUT 50 CNE A Ce ira EP HNPAV AVI 
52 CRE AN ENE ET RS ORNE D ARR CH MAR RUN ON ONE VI 
56 UT DATI EE PE V bis, VI 
PEN PSN RP EE 58 EE PS BR OI AN OT PE RE NE RE A DO ES OS Om OP SI V bis, VII 
60 AR Sn ni me nee ee ns Bts ST Une EUR VIII 
62 D Le CEE non NU DE nl À En ne OLe no el a CR V bis 
64 : 245/2 
DOP MERE E RES 68 CRE MMS Rs VII, VIII 
70 : 3s1/2 
72 le A SA OP D TE D Er RAC EE D AN ee VII 
74 M RP A M NET A D NO CR VIII 
1h11/2 COLE CT RS Éd IX, VII bis 
84 DÉOU ere e DR ER MAR AS REA ne eh en ERA RE XIII 
88 RE ES PA RE le mèe Tite DR RS A EN D Ce XIII 
92 : 1h9/2 
100 : 2f7/2 
106 : 2f5/2 
410 : 3p3/2 
412 : 3p1/2 
120 1e Se At dE MAO ne D ER DIE NN SN DRE XI 
122 RSR RSR A DRASS. dre el de nn XI 
A NE 0 CE Ne ER RE D RACE RS RE CE XT 
13% m, à moins que Po217 ait : Qg — 1,8 MeV XII 
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CORRÉLATION ANGULAIRE BÊTA-GAMMA A L'ÉNERGIE DE RÉSONANCE 


Par Maurice SPIGHEL, 
Laboratoire de Synthèse Atomique, G. N. R. $., Ivry-sur-Seine. 


Résumé. — Étude des conditions à remplir pour réaliser une expérience de corrélation ne 
laire bêta-gamma à l’énergie de résonance. Les calculs antérieurs d’autres auteurs sont complé- 
tés, en insistant sur le calcul de grandeurs directement accessibles à l'expérience. 


Abstract. — The author studies the conditions for the success of an experiment of beta-gamma 
angular correlation, with resonance fluorescence of the gamma-ray. Additions are made to 


previous calculations of other authors. 
accessible to experiment. 


Introduction. — L'étude de la corrélation angu- 
laire entre un rayon bêta et un rayon gamma émis 
en cascade et diffusé par résonance a déjà fait 
l’objet de divers calculs théoriques [1-21, quand la 
transition bêta est permise. Une expérience de ce 
type a l’intérêt de fournir des renseignements sur 
le type de l'interaction responsable de la désinté- 
gration bêta. Elle permet en principe de distinguer 
entre une interaction de type V et À ou Set T. 
Les difficultés expérimentales d’une telle mesure 
sont assez considérables. Cet article analyse les 
relations entre les grandeurs mesurables et les 
grandeurs théoriques qui peuvent être atteintes ; 
il montre que des additions doivent être apportées 
aux calculs cités [1-2]. 

Une telle expérience apporte des renseignements 
théoriques à peu près identiques à ceux fournis par 
une expérience de corrélation bêta-noyau de recul. 
La diffusion du noyau de recul par la matière per- 
turbe sa direction et son énergie initiales et fausse 
les conclusions que l’on peut en tirer [3] ; cette dif- 
fusion est sans importance dans une expérience du 
type que nous considérons si on détecte le gamma 
à la résonance lorsque la condition 7, Te est 
satisfaite, 7, et T° étant respectivement le temps 
de vie du noyau émetteur et le temps de collision 
de ce noyau dans la matière de la source. 

Les contradictions [4] entre les différentes expé- 
riences relatives à la corrélation angulaire bêta- 
noyau de recul et à la non-conservation de la 
parité sur n?, ‘He, li, PNe, Ne, A -152Eu, 
semblent se lever petit à petit [5]. Il est admis à 
l'heure actuelle que l’interaction bêta est du type V 
et À. Il serait cependant intéressant d’avoir par la 
méthode de corrélation angulaire bêta-gamma à la 
résonance, une indication sur le type d'interaction. 

Cet article analyse aussi certaines conditions qui 
doivent être remplies pour la réussite d’une telle 
expérience. Le calcul de la fonction de corrélation 
angulaire bêta-gamma à la résonance a été fait à 
un facteur de proportionnalité près [1-2]. L'expé- 


Emphasis is put on the calculation of quantities directly 


rience peut se schématiser par la figure 1. Il est 
nécessaire de connaître les deux grandeurs sui- 


les 4e 


vantes, qui correspondent aux évolutions succes- 
sives du gamma avant sa détection : 


a) La probabilité P,,5(0, W) dE qu’aun gamma 


d’être à l'énergie de résonance Ex à dE près quand . 


il est émis simultanément avec un bêta d’énergie W 
et quand les directions du gamma et du bêta font 
entre elles un angle 6 ; le taux de coïncidence est 
proportionnel à cette grandeur ; 

b) Le taux de coïncidences bêta-gamma diffusé 
par résonance : pour que l’expérience soit possible, 
il faut que ce taux ne soit pas trop faible. 


Caleul de P,,5,(0, W). — Les notations employées 
sont les mêmes que dans la référence [1]. Pour une 
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transition permise, Lewis et Curtis ont calculé 
la forme de la fonction de corrélation angulaire 


ce qui donne pour P,,»r(0, W) 


M QE [: apE 


P 0 — 


xp? ; 
Wa cos 0 — Wa cos o] 


(1) 
dans le système % —=m, =C —1:;q —W,—W : 
M —= masse du noyau. 
W, et E — énergies totales des transitions bêta 
. ét gamma. 
.__ Wet p — énergie totale et quantité de mouve- 
ment de l’électron. 


x = (B + 2/3.C)JA 


À, B, C'sont calculés en fonction des paramètres 
nucléaires des transitions bêta et gamma en se 
servant des références [1] et [6] par exemple. Ils 
dépendent : des constantes de couplage de la radio- 
activité bêta : Cy, Cs, Ca, Cr ; des matrices de 
transition Gamow-Teller et Fermi Mor et Mr; 
des Spins J1, Jo, Js des noyaux initiaux, inter- 
médiaires et finaux ; des amplitudes M, et E, de 
la transition gamma. 

Il faut ajouter à l’équation (1) les inégalités 


E < VWè — 1 (2) 
qg >|E + cos 0.VW?—1| (3) 
La formule (1) ne se différencie de la formule (A) 
que par des considérations cinématiques sans faire 
intervenir les caractéristiques nucléaires des émis- 
sions bêta et gamma. Ilsuffit de démontrer (1) dans 
un cas particulier 0 = x où nous faisons arbitrai- 
rement « = 0 : pour cela, nous posons B —0 ce 
qui entraîne l’émission isotrope du neutrino par 
rapport au bêta et € — 0, ce qui entraîne l’émis- 
sion isotrope du gamma par rapport à l’émission 
du bêta et du neutrino. 


Remarques.— Nous aurions pu démontrer la for- 
mule (1) en cherchant la valeur de Ce qu’il fallait 
pour normaliser l'intégrale 
SE (a—pcos0—Z) 

RIM. (—a—pcos0 -E) 


xp 
[1 Wa 


Cte, 
(£ “ie e) cos 0 Wa cos? o] dx 
x est le déplacement d’énergie du gamma à partir 
de la valeur Æ?/M nécessaire pour qu’il soit à 
l'énergie de résonance et causé par l’émission du 
bêta et du neutrino. Les limites de l’intégrale défi- 
nissent le domaine de variation x. La fonction à 
intégrer représente la fonction de corrélation entre 
un bêta et un gamma porté à l’énergie (E?/M + x). 
Quand qg tend vers zéro, d’après l’inégalité (3) la 
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résonance n’a lieu que pour une valeur déterminée 
de 0 ; à cet angle déterminé, dÆ tend vers zéro. 

Jusqu'à présent, il n’a pas été tenu compte pour 
établir (1) de la largeur naturelle de la raie et de la 
largeur Dôppler d’agitation thermique. 

Le terme remplaçant Cfe fait intervenir le fac- 
teur g7? qui ne figurait pas dans la formule (A). 
Le taux de coïncidence bêta-gamma àla résonance, 
à énergie W fixe, mesuré en fonction de 6 
suit toujours la loi donnée par la formule (A); à 
angle fixe en fonction de l’énergie W de l’électron, 
il ne suit pas la loi indiquée par (A) car il faut 
tenir compte en plus de ce facteur g—1 variable. 


Probabilité du processus de résonance. — La 
probabilité y9 qu’a un rayon gamma à l’angle 6 en 
coïncidence avec un rayon bêta d'énergie W d’être 
diffusé par le processus de résonance est, en tenant 
compte d’une largeur de raie quelle que soit sa forme 
et son origine, naturelle ou thermique : 
ap 


os 0 
Wa Wa 


ap? 


Me In 
4 = Mr | 


SE cos? 0]. (4) 


La section de diffusion au voisinage de la réso- 
nance étant : 


o= gXT2.(8x)-1.[(E — Ep)? + T?f4]— 
g = (1 + 2Jexe) (1 + 25) 


— facteur statistique, 

— longueur d’onde du gamma, 

— nombre d’atomes diffuseurs par cm?, 

— largeur naturelle de la raie, 

Jexe — spin du niveau émetteur gamma, 

J; —spin du niveau fondamental du noyau 
diffuseur où aboutit la transition 
gamma. 


HS >0 


l'est liée à la vie moyenne 7, de la transition 
gamma par la relation : ['.7, — 1. La formule (4) 
n’est valable que si la largeur de la raie gamma 
d’origine naturelle ou thermique est petite devant 
le déplacement d’énergie £?/M du gamma à la 
résonance. On peut alors considérer P,,:8(0, W) 


constant sur toute la largeur de la raie. 

Pour étudier l’ordre de grandeur du taux de 
coïncidence bêta-gamma à la résonance, il ne sera 
considéré que xp. # doit être choisi d’après l’équa- 
tion (4) pour que y, ne soit ni trop petit ni trop 
grand. 

— Si xp est trop petit, l'expérience sera impos- 
sible à réaliser vu le faible taux de coïncidences. 

— Si xp est trop grand, il n’y a plus de dis- 
tinction entre les gamma qui sont portés à la réso- 
nance et ceux qui n’y sont pas ; si Z est choisi 
trop grand, le diffuseur peut aussi diffuser les 
gamma qui sont assez loin de l'énergie de résonance. 

{x doit être ajusté par le choix de 7 à quelques 
centièmes. ; 

Si le processus de diffusion élastique prépon- 
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dérant nonrésonnant en compétition est la diffusion 
Rayleigh, estimons d’après Storruste [7] la section 
efficace 68 


or < 10—%5.(0,8/Æ).(Z/82) 


pour un angle de diffusion plus grand que x/4. 
. La condition pour que la diffusion par résonance 
soit prépondérante est 


q Z\? 1 
IS SAOE 7e PAR 5 
JR TT Te (5) À (5) 


Aw et Z étant le nombre de masse et le nombre 
atomique du noyau, À étant l’abondance isoto- 
pique du noyau considéré dans le diffuseur. 
Écrivons la condition pour que la largeur de raie 
Dôppler d’origine thermique, relative aux noyaux 
émetteur et diffuseur, L'% soit petite par rapport au 
déplacement d’énergie du gamma £?/M pour être 
à la résonance. C’est aussi la condition pour que 
l'énergie de recul du noyau après émission bêta 
nécessaire au gamma pour être à la résonance, soit 
grande devant l’énergie d’agitation thermique. 


Tr est calculé en tenant compte de la vitesse rela-. 


tive d’agitation thermique des noyaux émetteur et 
diffuseur gamma en les supposant à la même tempé- 
rature (1) 


Dre = E.Amt?.1,4.10—5 


à la température ambiante normale. La condition 
s'écrit : 
E > A°.2,8.10—2 (6) 


Si la condition (6) n’est pas réalisée, la for- 
mule (4) n’est plus valable pour interpréter l’expé- 
rience. 


Remarques : 


1) Dans Particle [1], la largeur de raie due à 
l'agitation thermique est 0,02 eV pour la raie 
de 436 keV du ?Ne alors que nous trouvons 
Vyr — 1,25 eV en accord avec la référence [8]. 

2) yrp permet de déterminer F, donc 7. 

3) 1, € Test généralement une condition facile 
à satisfaire pour une source gazeuse. Elle conduit à 
l'inégalité suivante si'L est le libre parcours moyen : 

W? —1\1? 
> re FLO 

Taux de coïncidences bêta-gamma à la réso- 
nance. — Soient «, l’angle solide du diffuseur 
gamma par rapport à la source, «, l’angle solide 
du compteur gamma par rapport au diffuseur, 
e l’efficacité photoélectrique du cristal détectant 
les rayons gamma, a la proportion de bêta pris dans 


(:) Cette façon de traiter le problème est une simpli- 
fication. En toute rigueur, il faut faire intervenir la tempé- 
rature de Debye dans le solide constitué par le diffuseur, 
Voir Lame (W. E.), Phys. Rev., 1939, 55, 190. 
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la bande d’énergie étudiée par rapport à tout le 
spectre, et b la proportion de l’embranchement 
bêta-gamma considéré par rapport à tous les 
bêta, le taux de coïncidence bêta-gamma à la 
résonance { est : 


le =. 01 2 € 4.070 (7) 


k est un coefficient voisin de 1 dû à la corrélation 
angulaire gamma à la résonance-gamma diffusé par 
résonance [8]. k est le rapport de la probabilité que 
le gamma soit diffusé sous l’angle solide «2 à 
l’angle considéré, et de la probabilité moyenne 
qu'un gamma soit diffusé sous un angle solide &2. 
À la suite d’une transition bêta permise, le gamma 
n’est pas polarisé rectilignement : la corrélation 
angulaire entre un gamma à la résonance et le 
gamma diffusé par résonance ne change pas avec 
l’angle 8 entre les directions d'émission du bêta et 
du gamma (9). 

Dans les meilleures conditions, estimons la valeur 
de -# maxima., Si l’on "choisit co 10 
Gy 410%; a =10 Pb = AO00 EE 
Y —50/ ontrouvet A0 

Regardons grossièrement l'influence de l’éner- 
gie Æ du gamma sur &. y,» Varie comme ÆT#, 
À variant comme Æ£71. Pour un M, pur, l varie 
comme Æ*%{[101] : y.x est indépendant de l’énergie 
du gamma pour un M,. Pour un ÆE;,, l' varie 
comme Æ5 : y, Varie comme Æ?. € varie à peu 
près entre 100 keV et 1 MeV comme Æ164 
dans INa ; donc pour un £, pur il y a intérêt à 
RE E assez grand et pour un M, pur, E assez 
faible. 


Influence de l’état chimique de la source. — 
Il n’y a pas de difficultés si la source gazeuse est 
monoatomique, mais on est limité aux gaz rares 
radioactifs. Si la source gazeuse est polyatomique, 
il faut tenir compte des processus de dissociation 
et de choc entre le noyau de recul et le reste de la 
molécule, qui perturbent le recul du noyau ; T, est 
souvent plus grand que le temps nécessaire à 
rompre les liaisons chimiques (& 1071 sec) [11]: 
entre l’émission bêta et l'émission gamma, le 
noyau parcourt une distance plus grande que les 
dimensions moléculaires. L’énergie de liaison Or 
du noyau de recul dans la molécule excitée après 
radioactivité bêta avec changement de Z d’un des 
constituants de la molécule, est très faible. Après 
le choc avec le reste de la molécule, la direction 
et la grandeur de la vitesse du noyau de recul n’est 
sensiblement modifié que si la masse du reste de la 
molécule est grande vis-à-vis de la masse du noyau 
de recul. Ces quelques considérations limitent 
pratiquement le choix des radioéléments pour 
l'expérience corrélation angulaire bêta-gamma à 
la résonance aux deux cas suivants : à) gaz rares 
radioactifs ; D) radioéléments qui peuvent former 
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un hydrure gazeux à basse pression aux tempé- 
ratures ordinaires. 

Comme x doit être faible devant l’éner- 
gie E?JM, il importe peu de savoir s’il y a ou non 
dissociation de la molécule. 


Remarque. — Dans l'expérience de Goldhaber et 
autres [12] avec une source solide ou liquide 
sur Eu, 7, & 3.107 sec, temps plus petit ou 
de l’ordre du temps nécessaires au noyau de recul 
pour parcourir les distances intermoléculaires et 
intramoléculaires. Le raisonnement précédent n’est 
donc plus valable. 


Conclusion. — Le choix des émetteurs bêta à 
transition permise qui peuvent servir à une expé- 
rience de corrélation angulaire bêta-gamma à la 
résonance est très limité. Les inégalités rencontrées 
montrent que l’énergie du gamma ne doit être ni 
trop faible, vu l’agitation thermique, ni trop forte 

comparée à l’énergie de la transition bêta. La vie 

moyenne du noyau intermédiaire émetteur gamma 
doit être suffisamment faible pour que les processus 
de diffusion en compétition ne soient pas prépon- 
dérants. L’abondance isotopique du noyau final 
qui servira de diffuseur ne doit pas être trop faible. 
La proportion de l’embranchement bêta-gamma 
considéré doit être grande. Les spins des noyaux 
émetteurs bêta, émetteur gamma, et diffuseur 
doivent être connus aïnsi que, dans une moindre 
mesure, la proportion de caractères £, et M; de la 
transition gamma. De plus l’état chimique de la 
source limite le choix des radioéléments possibles. 
Les quelques cas possibles sont ?F, ?%Ne, SKr, 
123$n. Le plus favorable est ?3Ne bien que ne soit 
pas connu le rapport des matrices Fermi et Gamow- 
Teller : il peut y avoir ambiguïté car l’interpré- 
tation peut dépendre de ce rapport. #F nécessite 
un diffuseur constitué par du néon liquide. L’utili- 
sation de %Kr et 12%Sn dépend de la valeur de la 
vie moyenne du niveau intermédiaire émetteur 
gamma : cette valeur n’est pas encore connue. 
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Curie pour sa bienveillance et ses encouragements. 
Je remercie M. F. Suzor, Sous-Directeur du Labo- 
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Appendice. — Quand la condition (6) n’est pas 
respectée, 1l faut tenir compte de l'agitation ther- 
mique dans la source gazeuse à la température 7, 
et dans le diffuseur solide à la température de 
Debye Th. Dans ces conditions 4e a l'expression 
suivante qui remplace la formule (4) : 


ve Me 22 ln 


8£q ; 
ap CADET l À 
. Li re 0 me cos 0-5 (ert x, + erf *;) 
« p cos 0 [2k(Ts + Th M\ 24 , : 
en posant 
PR COS UN L q p cos 0 —Æ 
 LERUTSE Tp).M]72 D (9X(Ts To) M1 


Manuscrit reçu le 5 novembre 1958. 


Note ajoutée à la correction des épreuves. Nous avons 
ris connaissance de l’article théorique de N. A. Burgov 
13]. La première partie de cet article traite le même pro- 

blème. La formule (11) de N. A. Burgov est identique à 
notre formule (4’), mises à part les notations, à un coefñ- 
cient deux près. Sa formule (11) donne pour la section 
efficace X9/n une valeur deux fois trop grande car sa for- 
mule (3) contient un facteur deux en trop. 

Il faut d’autre part remarquer que son terme À est uni- 
quement lié à la corrélation angulaire bêta-neutrino et est 
indépendant du caractère E, ou M, de la transition gamma 
qui suit : N. A. Burgov ne tient pas compte de la corréla- 
tion angulaire entre la direction de l’émission du gamma 
et les directions des émissions bêta-neutrino. Cette corré- 
lation angulaire existe, mais intégrée sur toutes les direc- 
tions du neutrino la corrélation angulaire résultante bêta- 
gamma pour une transition bêta permise est nulle. 
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CALCUL DE LA DIFFUSION ÉLASTIQUE SANS POLARISATION ET SANS ÉCHANGE 
D'UN ÉLECTRON DE VITESSE NULLE 
PAR UN ATOME NEUTRE DE SODIUM DANS L'ÉTAT FONDAMENTAL 


Par AzBerT SALMONA et PrEeRRE FRENKIEL, 


Institut Henri-Poincaré, Paris. 


Résumé. — Les auteurs calculent la section efficace de choc d’un électron très lent contre un 
atome de sodium dans l’état fondamental sans tenir compte des effets de polarisation et d’échange. 


Abstract. — The cross section of the collision between a slow electron and a sodium atom in 
the ground state is computed without taking into account the effect of polarisation and of exchange. 


Introduction. — Nous nous proposons de calcu- 
ler la section efficace de choc entre un électron 
lent et un atome de sodium. On peut montrer en 
se fondant sur un principe variationnel, les fonc- 
tions Gr(l'Îr) étant seules « variées », que l’on 
aboutit au système d’équations intégro-différen- 
tielles suivant 


dé Upiip 1) RUES 
Fe - ET à à] Gr(T/r) 
— 9 2 | Vrr” — Woyr] Gr(1”}r) (1) 
avec 
3 1 ,e . LT 
EE Etes AM =; Er +ShF, (2) 


où [” désigne la voie d’entrée, F la voie de sortie, 
Er l'énergie de l’atome dans l’état l”, A$ celle de 
Pélectron diffusé, avec 


% . 
Vrr/(rn+1) = foret 
4 N 
e TN+1 er FN+1,i 


] bre[riN +1] dr + De (3) 


Wrr- Gre(L'frn £a) = TN: 
, : “ea 
X Nf D [r(A ED] [I — F1] bref] — Gre(lfrx] drNH, 
LN 


(1) 


brfr®] désignant le produit de fonctions anti- 
symétrisées de l’atome par des fonctions angulaires 
et de spin de l’électron 1. 4 désigne l’hamiltonien 
du système complet et E l'énergie totale ; 
Tnt = [ni — Ti. 

Les formes asymptotiques des fonctions G sont 
exprimées moyennant la matrice de réactance À. 

Ces équations étant évidemment très difficiles 
à étudier, 1l est intéressant de résoudre des équa- 
tions plus simples — même si les hypothèses phy- 
siques que ces dernières supposent sont physi- 
quement inacceptables — afin de posséder un sys- 
tème de fonctions qui puisse servir de point de 
départ au caleul par approximations successives. 


Cas du sodium neutre entrant en collision avec 
un électron très lent. — La section élastique est 
seule à considérer si l’électron n’a pas suffisamment 
d'énergie pour porter l’atome (initialement à l’état 
fondamental) au premier état excité. La confi- 
guration du sodium neutre à l’état fondamental 
est (152 2s22p6 3s), c’est-à-dire 3 couches complètes 
et un électron de valence. Dans ce cas, seule l’onde 
de diffusion saura une forme asymptotique finie qui, 
dans le cas où Ær — 0, est ar + b, les ondes de 
diffusion p, d'etc., décroissant exponentiellement. 
Si on néglige l'effet de la déformation des ondes p, 
d, etc. sur l’onde s, c’est-à-dire, en particulier, la 
polarisation contenue dans les termes (3), et si l’on 
ne tient pas compte de l’effet d'échange représenté 
par les termes intégraux (4) dans les équations (1), 
le système (1) se réduit à une seule équation diffé- 
rentielle 


qd? lr'Ur- Le 1) | * 272 
= - L D Pareil) E 1. Cr) 0 
où 


le = 0 si Jp 0, 


Les fonctions de l’atome de sodium étant tabulées 
par Hartree, d’une part, par Fock et Pethrashen, 
d'autre part, il est possible de calculer le poten- 
tiel — 2Vrr qui n’est autre ‘que le potentiel 
central de l’atome dans l’état 3s auquel est soumis 
l’électron incident. 

L'intégrale (3) dans la configuration 


(1522522p6 3s) 
donne en effet 
1 
Vrr(r) — Aer [1 — 2%, (15, 1sfr) — 2 Y,(2s, 2s/r) 


700 Yo(2p, 2p|r) Fa (35, 3s|r)] (5) 


N°0 4 
avec 


Folalnllr) =7r 1 P{nlli) > Pulse) de 


Tan Le Pile) . Pro) dt, (6) 
les P(niir|) étant tabulées, 
Yolnt, nll0) = 0,  Yi(ni, nlloo) = 1, (7) 


les fonctions Y (nl, nllr) peuvent être calculées 
numériquement soit directement soit — ce qui est 


préférable — par résolution de l’équation diffé- 
rentielle 
dF P?{(nllr 
Te tr nlir _ (8) 


qui est, moyennant (7), équivalente à (6). Les équa- 
tions (8) ont été résolues numériquement en utili- 
sant la méthode de Gauss-Jackson, à l’aide de 
machines à calculer mécaniques. Ensuite, par la 
même méthode — fondée sur les différences finies 
— nous avons résolu l’équation (1") pour obtenir le 
coefficient D dans la forme asymptotique de 
GrA{l'}r) qui est a(r + D). La section élastique 


CO 
totale étant Q — Fe > (21 + 1) sin? », nous avons, 
0 


avec les hypothèses précisées plus haut, les relations 
suivantes : m— A, k%+1, si k +0, n9 —kD = A, k, 
d’où Q(k — 0) — 47 D?. Les résultats numériques 
donnent D — 2,36, Q # 22 raÿ, a, étant le rayon de 
l'atome d'hydrogène. Les mesures expérimentales 
laissent prévoir une section efficace qui est proba- 
blement bien plus élevée (de l’ordre de 400 x &), ce 
qui implique qu’il est indispensable de tenir compte 
des effets de polarisation et d'échange. Aussi 
nous sommes-nous actuellement attachés à cal- 
culer les corrections dues à l’échange et à la polari- 
sation sur la fonction Gr-(l”|r), solution de (1) qui 
est prise comme première approximation. 
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Notalions utilisées : Représentation 
y = 4LSM; Ms ln ns, 


où aLS M, Ms représentent les nombres quantiques 
de l’atomeet lmy,m, ceux de l’électron diffusé. 

Représentation [' — aZLILT ST MT ME, où L' 
et S7 sont obtenus par couplage de L et de Z d’une 
part, de Set du spin de l’électron diffusé, d’autre 
part ; z désignant le 5% électron r; = (r;, r;), 
où r;, désigne l’ensemble des coordonnées angu- 
lares ; &, = (r;, ©) — (r;, ri, oi), o; étant le spin. 

Si l’on considère un système à N + 1 électron 
(N pour l’atome), le symbole 7° désigne l’ensemble 
(Li, La, Li. Mix, . .. yyx) À (ri, Gi), C'est-à- 
dire toutes les coordonnées sauf r;. dr désigne de 
même l'élément de volume sur toutes les coor- 
données sauf r;. Dans tout ce qui précède, les 
fonctions de l’atome sont supposées exactes. Les 
problèmes de choc se ramènent essentiellement au 
calcul de la ‘fonction de l’électron diffusé Gr(T"|r) 
où [” désigne la voie d’entrée du système, et [°, la 
voie de sortie. 

Les équations sont exprimées dans le système 
d'unités atomiques où e = = m —=c—=1. 

Ce présent travail est le premier résultat obtenu 
par l’équipe de chercheurs français dirigée par 
M. T. Kahan dans le cadre du programme établi 
conjointement par ce dernier et un groupe de 
chercheurs anglais que dirige le D' Seaton de 
l’University College de Londres, groupe qui depuis 
de nombreuses années s’est spécialisé dans les 
calculs des sections de chocs atomiques. 

L’équipe française tient à exprimer sa vive recon- 
naissance à M. Coulomb, Directeur du Centre 
National de la Recherche Scientifique pour l’aide 
et l’intérêt qu'il porte à ses travaux. 


Manuscrit reçu le 16 janvier 1959. 
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Conférence Internationale sur la Physique des Semi- 
conducteurs, Prague, 1960. — L'Académie tchécoslovaque 
des Sciences organise une Conférence Internationale sur la 
Physique des Semiconducteurs qui aura lieu du 29 août 
au 2 septembre à Prague sous le patronage de l’Union 
internationale de Physique Pure et Appliquée, et qui fera 
suite à celles de Reading en 1950, d'Amsterdam en 1954, 


de Garmisch-Partenkirchen en 1956 et de Rochester en 
1958. 

Secrétaire du Comité Organisateur : Dr. Miloÿ Matyéÿ, 
Institut de physique appliquée de l’Académie tchécoslo- 
vaque des Sciences, Prague 5, Cukrovarnickä 10, Tchéco- 
slovaquie. 
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EXPOSÉ ET MISE AU POINT BIBLIOGRAPHIQUE 


L'EFFET POYNTING-ROBERTSON 
par G GUIGAY, 


Observatoire de Marseille. 


Résumé. -- Soit une particule se mouvant dans le champ de gravitalion du Soleil et soumise à 
l’action de son rayonnement. Poynting et Robertson ont montré que la pression de radiation n’est 
pas exactement dirigée suivant la droite Soleil-particule, mais possède une composante opposée 
à la vitesse de cette dernière. Une démonstration de cet effet est donnée pour le cas d’un rayon- 
nement électro-magnétique et pour celui d’un rayonnement corpusculaire. à 

l’étude du mouvement de la particule montre que le demi-grand axe et l’excentricité de son 
orbite décroissent continuellement : la particule tombe en spirale sur le Soleil. Pour un petit corps 
sphérique de 10—? centimètre de rayon et de densité 4 le temps de chute est de l’ordre de 105 ans. 
Pour une densité constante, ce temps est proportionnel au rayon. 

Les conséquences de cet effet sur la disposition des orbites au sein d’un courant météorique, 
sur l’équilibre du nuage zodiacal.et sur la capture de matière météorique par la Terre sont passés 
en revue. 


Abstract. — Let a particle move in the sun’s gravitation field. According to’ Poynting and 
Robertson the radiation pressure is not exactely directed along the straight line Sun-particle, but 
possess a component opposite to the particle’s velocity. A simplified demonstration is given for 
the case of a solar electro-magnetic radiation and for that of a corpuscular radiation. It is shown 
that the semi-major axis and the eccentricity of the orbit of the particle decrease continually : the 
particle falls in spiral on the Sun. For a particle of 10—? centimetre radius and a density of 4, the 
time of fall is about 105 years. For constant density this time is proportional to the radius. 

À study is given of the effects on the distribution of the orbits in a meteoric current, on the 
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zodiacal cloud, and on the capture of meteoric bodies by the earth. 


1. — L'effet Poynting-Robertson joue un rôle 
important dans l’évolution d’un nuage de poussières 
entourant une étoile. Il entraîne la chute en spirale 
des particules sur l'étoile en un temps de l’ordre de 
quelques millions d’années. Cet effet a été signalé pour 
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la première fois/par Poynting, mais c’est Robertson [1] 
qui, en 1937, en a donné la démonstration fondée sur 
la Théorie de la Relativité et le Calcul Tensoriel. Nous 
nous proposons de donner ici une démonstration plus 
simple de cet effet et d'en examiner les conséquences 
astronomiques. 


Considérons une particule P se trouvant dans le 
champ de gravitation d’une étoile E et exposée à son 
rayonnement. Indépendamment de l’attraction new- 
tonienne G, P est soumise à la pression de radiation ®. 
Si P est immobile, ® est opposée à la direction de 
l'étoile (fig. 1 a) ; il n’en est plus de même, comme 
nous allons le montrer, dans le cas où la particule est 
en mouvement : On voit apparaître une composante f 
de la pression de radiation opposée à sa vitesse ®. 
(fig. 1 d). Il en résulte un freinage qui constitue l’effet 
Poynting-Robertson. 


_2. -- Nous admettrons que la particule est sphé- 
rique, de diamètre nettement plus grand que la longueur 
d'onde de la lumière incidente et parfaitement absor- 
bante. 

Considérons deux systèmes de référence : 

19 un système S lié à l'étoile E et constitué d’axes de 
directions fixes ; 

2° un système S', en translation uniforme par 
rapport à S et instantanément au repos par rapport 
à la particule P, c’est-à-dire ayant un vecteur-vitesse 
par rapport à S équipollent à celui de la particule à 
l'instant considéré. 

. Nous conviendrons d’affecter d’un accent les quan- 
tités relatives au système S’ et de ne pas accentuer 
les quantités relatives au système S. 

Pour l'observateur de S' la pression de radiation 
s'exprime par 


DATE (1) 


N°4 


A” désignant la section droite de la particule et£W” la 
densité de l’énergie du rayonnement. W' n’est pas un 
mvariant, sa valeur W dans le système S est reliée à sa 
valeur dans le système S' par la relation : 


où v est la vitesse de la particule pour l'observateur du 
système $ et s la densité du courant d'énergie dans ce 
système ; c représente la vitesse de la lumière [2]. 
La notation [v p] désigne le produit contracté 
VE Dia (3) 
du vecteur-vitesse © par le tenseur des tensions de 
Maxwell dans le système S. 

Mais si, dans le système S, on désigne par 4 le 
vecteur unitaire porté par la direction du rayon lumi- 
neux : 

o=—c.Wl (4) 
d'où 


Pme (r 2 lol: : (ep). (5) 


On voit que le second terme de la parenthèse est de 
l’ordre de v/c alors que le troisième est de l’ordre 
de v?/c?. Nous conviendrons de négliger les termes de cet 
ordre, de sorte que : 


12 
Te l 
re (: ] 


ou, à la même approximation : 
.l\2 
Et des "| (6) 


D’après l'équation (1) la pression de radiation s’ex- 
prime alors par : 


l\2 
D — AW (1 —?) , 


C 


Tel est le module de la pression de radiation pour 
l'observateur du système S”. En réalité cette pression 
est un vecteur porté par la direction du rayon lumineux 
pour ce même observateur ; son expression complète 
est : 


1e 
D — AW (1 —7) L (?) 


l désignant le vecteur unitaire porté par le rayon 
lumineux pour l'observateur de S’. 

Pour passer à la pression de radiation pour l’obser- 
vateur S il est nécessaire d'utiliser le théorème de la 
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Relativité des forces. Considérons pour un instant deux 
systèmes d’axes dirigés comme l’indique la figure 2, le 
système Ox'y'z' étant en translation uniforme par 
rapport à Oxyz. Dans ce cas une même force agissant 
sur un point instantonément au repos dans le système 
O'x'y'z' a, dans les deux systèmes, des composantes 
XYZet X'Y'Z' liées par les relations 

ANRT } 


CORYA 7 UN" 


où 


c’est-à-dire que pour obtenir la force relative au sys- 
tème Oxyz il faut ajouter au vecteur force relatif à 
O’x'y'z" un vecteur supplémentaire de composantes 


LE 


ONE : 

262 2 Ce 

; . : ne ; a 

A l’approximation indiquée ce vecteur supplé- 

mentaire est négligeable. Revenant à notre problème 
nous pouvons écrire : 


2 
D — D — AW ( _ v'. 


Pour achever d’écrire l’expression de ©, il faut 
exprimer l'en fonction de 4. Pour cela appliquons le 
théorème de la composition des vitesses au vecteur- 


vitesse de la lumière 
uw = cl u = c.l. 


La relation de composition des vitesses sous forme 
vectorielle s'écrit [3] : 


u © u Ù 
uw | JE LA r)v EL au ———.Ù 
U” v= 


c’est-à-dire : 


2(4 . v)| (Le De 


ou, en remarquant qu'à l’approximation indiquée 
« se réduit à l’unité : 


Fever led 


(9) 
Il reste à exprimer la force de gravitation de Pétoile 
qui comprend [4]: 
19 la force newtonienne classique 


u Le L 
Fo 


ou m est la masse de la particule, r la distance EP et pr 
le produit F. M de la constante # de la gravitation par 
la masse de E, 
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20 des forces supplémentaires : 
a). une force dirigée suivant le rayon EP et d’expres- 


SION ;: 
( 2u , 20  2r* Te | L 
"pè ep | çä c? CE (2u/r) \ 


ps dr Jde 


où 


b) une force dirigée suivant la tangente à la tra- 
jectoire de P dans le sens du mouvement et de module 
2ur lol 


Fe — (upr)] 


A l’approximation indiquée nous ne devons con- 


server que la force newtonienne et la force 
0) 


Pre L 
c? rè 
toutes les autres étant de l’ordre de v?/c? 
Finalement l’équation du mouvement de la parti- 
eule s'écrit : 


d dv Um Dre 
d R + SR LES 
æ (nv)  m. SN 5 sl 
DNA 
+ AW (e ———— 4): 
C C 
; ue au 
Si nous négligeons encore le terme Ta FL dont je 
rapport à la force newtonienne est de l’ordre de + 
il reste : 
v.l AW 
mer + an (1 nn v 
de 1? c F 
ou en posant 
AW = f (10) 
dv u m v.l Î 
— = — — l NE es l RE Surr Rio 
m U+ ii tie. (1) 


Le deuxième terme du second membre exprime la 
pression de radiation affectée d’un « facteur Doppler ». 
Au dernier terme correspond l'effet Poynting-Robert- 
son. D’après la démonstration précédente où nous 
avons utilisé les directions différentes de la lumière 
dans les deux systèmes de référence, il apparaît comme 
lié à aberration. 


8. — Whipple [5] a étendu l’effet Poynting-Robert- 
son au cas où la particule est exposée, non plus à un 
rayonnement électro-magnétique, mais à un rayon- 
nement corpusculaire. Notons que, d’après les idées 
modernes, le rayonnement corpusculaire des étoiles, 
consistant surtout en atomes et ions d'hydrogène est 
extrêmement important ; il serait même l'élément 
directeur de leur évolution. Par suite de ce rayonne- 
ment corpusculaire la masse d’une étoile diminuerait 
considérablement au cours de son existence : la masse 
du soleil par exémple aurait été à l’origine cinq à huit 
fois plus grande qu’aujourd’hui [6]. 

Conservons les mêmes notations que précédemment, 
remplaçant seulement la vitesse de la lumière par la 
vitesse g des corpuscules incidents (g' dans le sys- 
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tème S' lié à la particule). Le vecteur vitesse d’un 
corpuscule dans les deux systèmes s’exprime par : 


ga=—aql d'= qu. 


Nous admettrons que les vitesses v et q sont assez 
petites par rapport à c pour que les rapports v?[c? et 
qg2[c? soient négligeables. Dans ces conditions la loi de 
composition des vitesses et la loi fondamentale de la 
dynamique se réduisent à leurs expressions non rela- 
tivistes. 

La masse totale des corpuscules qui frappent la 
particule P pendant le temps dt est : 


A'd.g dt 


d désignant la masse totale des corpuscules dans un 
centimètre cube. Leur quantité de mouvement est 


A'd.g'? dt.v'. 


Si les corpuscules incidents restent, après le choc, 
unis à la particule, le théorème de la quantité de mou- 
vement donne : 

dv’ 


je A'd.q°?V. 


Le premier membre représente la force ® exercée | 


par les corpuscules sur la particule P. 
D — Ad.q?l — Ad.qd. 
Mais 
LIT 
g? = qg? + v? —2q.%. 

Admettons encore que la vitesse v de la particule est 
petite par rapport à celle des corpuscules incidents, de 
sorte que l’on peut négliger v?[q?, alors : 

De QUE) 
d’où 


D — Ad.qli — (q.v/q")] (a —®). (13) 


L’équation du mouvement de la particule s’écrit : 


dé q* 
ou 
av Um Te Ÿ 
ou, en posant 
ZCRTÉEE (14) 
av Um lv Î 
mare hf) 2e 


L’analogie des formules (11) et (15) apparaît plei- 
nement si l’on remarque que la loi d'équivalence 
matière-énergie fait correspondre à la densité d’éner- 
gie W une densité matérielle 


de sorte que l'équation (10) s'écrit 


1 AC?) 


N°4 ! 


(2). 


AE PEN TT A 


N° 4 


4. Mouvement d’uné particule dans le champ de 
gravitation du soleil. 


Introduisons les données relatives à la particule et 
au rayonnement. Nous ne tiendrons compte que du 
rayonnement  électro-magnétique. La Constante 
solaire S c’est-à-dire le flux reçu à la distance du 
Soleil où se trouve la Terre vaut : 


1,35 10% ergs par centimètre-carré et par seconde. 


* Si b désigne la distance de la Terre au Soleil l’éclai- 
rement à la distance 7 est : 


156? 
F2; 

Ce nombre exprime l'énergie contenue dans un 
cylindre de 1 em? de base et de longueur c, de sorte que 
la densité de l’énergie est : 

Sb° 
CT 


= 


L’équation (10) donne : 
Hd Some A SL 


j cr? 7? me? 
ou en posant : 
= 6, (16) 
mc” 
mc 
f = HE G. (17) 


Pour une particule sphérique de densité € et de 
rayon s, o est une constante 


EE 
MENTION 
Pour : 


b—/150410=%6èm c — 3.101 cmy/sec 


CR ES ao (18) 


L’équation vectorielle (11) s’écrit en coordonnées 
polaires : 


“(a 
7 l'a dé 


PAITTE r? TE r? dt 
4 (20%) mod? 
r dt dt r dt 
ou 
derr ) u—oc  2cçdûr 
de A que ne er EL rad 
dt dt, 7 (19) 
a ( :T) se 
dt dt dt 


Dans la première équation (19) le terme 


LOC, 


r? 
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est analogue à un terme de gravitation mais cette 
dernière est diminuée par la pression de radiation et 
remplacée par une « gravitation effective » où la quan- 
tité : 


= Ve (F — Gte de la gravitation) 


est remplacée par 


W'=up—oc= Fm—oc. (20) 


Nous considérerons le cas où 
oc < Fm 


de sorte que la gravitation l'emporte sur la pression de 
radiation. La deuxième équation (19) s’intègre immé- 
diatement et donne : 


(21) 


h étant une constante. L'expression au premier membre 
est le double de la vitesse aréolaire de la particule P. 
Cette quantité, qui serait constante dans un mouve- 
ment képlérien (correspondant à o — 0) diminue pro- 
gressivement lorsque croît l’angle © balayé par le 
rayon vecteur. À est la valeur du double de la vitesse 
aréolaire à l’origine des temps. Posons : 


He Res h — 69. (22) 


dt 


Quand Æ s’annulera, sans que r soit nul, nous aurons 


de 
ds 
de sorte que la vitesse de la particule se réduira à sa 
composante suivant le rayon vecteur et il en sera de 
même de son accélération d’après la 2€ équation (19). 
Avec les forces agissantes de telles « conditions ini- 
tiales » ne peuvent entrainer qu’un mouvement recti- 
ligne dirigé vers le Soleil ou dans la direction opposée. 
L’angle y ne peut excéder la valeur : 

k 


re 


ce qui limite le nombre AN des révolutions à 
N =. 


2 Ko 


(23) 


Robertson [1] étudie le cas d’une particule décrivant 
autour du Soleil une orbite circulaire de rayon : 


ro — une unité astronomique, 


sa vitesse angulaire n étant donnée pour la 3€ loi de 
Képler : 


fé 


DAT LA 
Dans ces conditions. 
h = nrè = r$ Vu’frè = Var, 


d’où 


N = Vurol2 To. (24) 


Dans un calcul approximatif on peut prendre pour u' 
32 
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sa valeur newtonienne déduite de l’application à la 
Terre de la 3e loi de Kléper, soit : 


4 T°? bÿ 
æ = "LE 


T désignant la durée de la révolution de la Terre : 
31,6.106 secondes. 
Avec la valeur déjà donnée de b, on trouve 


u—/1393.10%; 
Pour 
SEMI CIN 


c’est-à-dire, d’après (18) 
6 = 4,6.1010 


on trouve. 
N —1,5.108. 


Mais s peut être beaucoup plus faible et N est propor- 
tionnel à s d’après (24) et (18). Avec s — 10-$ centi- 
mèêtre, on trouve seulement : 


NET 1107 


Ce résultat justifie la conclusion de Roberston 
d’après laquelle l'effet étudié peut entraïînerla chute des 
particules sur le Soleil en un temps suffisamment court 
pour avoir une signification cosmogonique. 

L’équation de l'orbite s’obtient en éliminant t entre 
les équations (19) et (22). Posons : 


Re Aloe (25) 
L’équation (22) s'écrit : 
dofdt — Hu? (26) 


et la première équation (19) : 
due (& 


du 
LR ER Nr 4) D RS OT Et GC Es 
uw de | us \d Fes CRT RE 
ou, en remplaçant du/dt et d’u/dt? par leurs valeurs : 
dufdt — (duJde).(do/dt) — Hu?(dufde) 


RE CRUE 


dde dr/Mdr 
= re (a à ane (ES — mu () 
= w | u d à [ d 5 
dard HT 
de? */Hdo + U — HE (27) 


Robertson remarque que o/h est en général une 
petite quantité. Considérons par exemple une particule 
se mouvant à la vitesse de 30 kilomètres par seconde 
sur un orbite circulaire de rayon égal à une unité astro- 


nomique. Si nous prenons pour À la valeur r? LE trop 
faible d’après (22), nous obtenons 
Red) 7015408 240025 100 
d’où, d’après (18) 

cfh — (0,56/s C).10—$8. 
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Pour 
s — 10—3 cm "C—=55 
oJh #10—$ 
On peut donc développer les divers termes de l’équa- 
tion (27) suivant les puissances deT. Négligeant les 


termes en (o/h)?, on obtient : 


L’équation (28) admet l'intégrale générale : 


5 2 - 
= Be? cos (\/1 D 


2u’ © a’ 
HAE +E (1 


B et y étant des constantes 


Développant encore les divers termes suivant les 


puissances de o/h et négligeant les termes en (o/h)?, on 
obtient 


/ 2 
DA D Ut (1 


— (e— »)] 


_ ) cos 
A 
ce que l’on peut écrire : 


BR? [1 — co/2h] cos (® — y) 


1 + 


u' 1 + 26o/h) 
4 n° 1 
u'1 + (20h) 
Bh? C) 2 
+ (1 0) (1% à) cos (9 — 1 
4 h°? 26 
ml FR +) 
BR. 
1 + 7 [1 — 50p/2h)] cos (® — y) 
ne A : (30) 
&  — (26p/h)] 


L’équation (30) peut être ramenée à l’équation 


képlérienne : 
1 e cos Ps 
Re RE ER re) (31) 
P 
en posant : 
en 
19 uw Fi 
(3 
BR? 56 (P2) 
Ne (1 D +) 
u 2h 


Les quantités p et e jouant le rôle d'éléments oscu- 
lateurs de l’orbite de la particule, les expressions : 
__26h9 N __5Bho@ 


her tie ne 
be 2, ue 


(33) 


+ "RON A 


No 4% 


représentent les perturbations du premier ordre en /h. 
On en tire 


Spip = —20p]h ele — — (5/2) (ofh).æ. (34) 
L’équation de définition du paramètre : 


p = a(1 — e?) 


permet de calculer la perturbation du demi-grand axe 


Op — da(1 — e?) — J2ae De 


Ôp 2ae 
Ô = a ——_—_—_— 
& De aqua 


ou, en remplaçant Ôp et de par leurs valeurs 
pure te (35) 
L’angle @ ne croît pas proportionnellement au 
temps mais s'exprime par: 
® — 0 + nt + des termes périodiques. 


Nous obtiendrons les inégalités séculaires de p, a, e 
en réduisant ® à : 


P = Po + nt 
que l’on peut écrire plus simplement 
o —= nt 


en choisissant convenablement la direction de l’axe 
polaire. L'expression des inégalités séculaires est donc: 


Le 2. L 8e? 
= rie 


nt 


5 o 26 
D; ent dpi y pnt TACTON 


2h 
Puisqu'on ne désire que les inégalités séculaires du 
premier ordre en c/h on peut prendre pour 4e, p,nla 
valeur de ces éléments pour { — 0. Une relation connue 
du problème des deux corps donne : 


p = hu 


(puisque k est la constante des aires à l’époque { — 0). 
La 3° loi de Kepler donne : + 


h2 
Dr Ne Mie 
DEN TE ne 
h ; 
ne (36) 
ai VAR? 
d’où > 
5 2 3e 
$e — : = DEQU— = T TE (37) 
2 a? VA — e? ai — e°)# 


Les équations (37) ne sont pas valables pendant un 


: eue (9 
temps très long car certains termes négligés, en _ 


étaient multipliés par ® ; ils risqueraient de ce fait de 
devenir importants au bout d’un temps suffisant. On 
peut, avec Wyatt et Whipple [7] les appliquer pen- 
dant un temps Ôt suffisamment court (temps qui peut 
cependant comprendre plusieurs révolutions). On à 
alors 


Be = —— ———— Ûùt 
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Revenant ensuite aux variations de aet e pendant un 
temps très long on peut identifier les quotients Se/3t 
et Sa]Ôt avec les dérivées de/dt et da/d£. On obtient 


de 5 ce da 


; 6.2 + 3e? 
(OI, 


To (1 — e?)#2 


vas dt ni 
où a et e ne sont plus considérées comme des constantes 
au second membre. Ces formules montrent que aet e 
décroissent continuellement ; la particule tombe en 
spirale sur le Soleil. 

La première équation (38) admet la solution parti- 
culière : 4 


; C0) 
qui entraîne 
dafdt = — 26]a 
‘a da = — 20 dt 
il 
3 (a? ss aë) = — 26. 


Cette solution concerne le cas d’une orbite primitive 
circulaire. La trajectoire demeure circulaire mais son 
rayon décroît continuellement. La particule tombe sur 
le Soleil au bout d’un temps : 


t = aÿl4o 
ou, en remplaçant 6 par sa valeur (18) 


ie TAONÉ 
ii RDS STORE (89) 


Cette formule a été établie dans le système CGS. Si 
l’on exprime 4, en unités astronomiques (15.10!? centi- 
mètres) et ten années (31.6.1056 secondes) on obtient : 


lans nn 7va. S. OPIOS (40) 
s restant exprimé en centimètres et € étant la densité 
de la particule par rapport à l’eau. 

Dans le cas d’une orbite excentrique, on tire des 
équations (38) 


da  2a(2 + 3e!) 
de M5 'e(d 16?) 


, da 2(2 + 3e) ;, &de  5de  5de 
"a  ell — e) E MERE 
d’où 
Kesl5 
CR apr (41) 


K est une constante que l’on peut déterminer en 
appliquant cette équation à l’époque origine. Alors : 


er Keël5 Es doll — ei) = Goll + €) (42) 
MINE e4l5 cale 


G étant la distance périhélie à l’époque t = 0. 
En portant la valeur (41) de a dans la première équa- 
tion (38) on obtient 
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Si l’on calcule des tables de la fonction 
* £ au LL, 
Ge) — À (1 — e?)#/2 de (an) 


la formule (43) permet de calculer le temps nécessaire 
à une diminution donnée de e. Utilisant la valeur (18) 
de o, on obtient dans le système CGS : 


DIHCISE 
5.2,53.101 


l Lo= [G{e) en G(e5)] 


ou, en remplaçant X par sa valeur (42) 


DEEE 


es 1 + en) 60 5 SI Geo) — G{e)]. 


1 — 
Si l’on convient d'exprimer 1-4, en années et qÿ en 
unités astronomiques en trouve : 


(t — to)ans = 1,13.107 va (1 ner ess 
Stem) SLE(E0o) — G(e)] 


e diminue toujours au cours du mouvement et l’équa- 
tion (41) montre que as’annule avec e. Le temps néces- 
saire à la particule pour tomber sur le Soleil s’obtient 
en posant 


(45) 


EN) G{e)=10; 


L’équation (45) donne alors 
(Et — to)ans = 1,13.107 Gualt + 20)? é "5 Sem) € G(e0) 


Wyatt et Whipple donnent pour diverses valeurs 
de e,, les valeurs correspondantes de : 


Er 
107s € qè 


t-t, étant exprimé en années, soit : 


2 tt : t— to 
107 SC qé , 107 SC q$ 
0,00 0,70 0,60 2,49 
0,10 0,86 0,70 3,16 
0,20 1,04 0,80 4,42 
0,30 1,27 0,90 7,29 
0,40 1,55 0,99 28,89 
0,50 1,92 


Considérons une particule telle que : 


SH CmM Ge 


décrivant initialement une orbite circulaire de rayon 


do — o — 1uA 

Alors 
t— to 
107 SC qè 
t— 1 = 2,8 107 ans. 


SA — 0,7 (voir tableau) 


Considérons maintenant la même particule décri- 
vant une orbite originelle de demi-grand axe égal à une 
unité astronomique, mais d’excentricité 0,7 ; le tableau 
donne : 

t— to 


res = 19) 10, 
107.sÈ gè 


N° 4 


PHYSIQUE 
Mais : . 
do = aol —— €0) = se qÿ = 0,36 
119 —6MD:0 0 A0 .107 ans. 


Ainsi, le temps de chute a diminué lorsque, con- 
servant la même valeur du demi-grand axe, nous avons 
augmenté l’excentricité. Si nous avions attribué à la 
particule, dans les deux cas le rayon 


s — 1072 centimètre, 


nous aurions trouvé respectivement pour durée de la 
chute sur le Soleil : 


2,8.105 ans ATOME: 

5. — L’effet Poynting-Robertson et l’Astronomie 
météorique — Wyatt et Whipple [7] ont calculé la 
durée de chute sur le Soleil des particules qui consti- 
tuent les principaux courant météoriques. Les valeurs 


utilisées de 4, et e, sont celles des comètes génératrices. 


Ils trouvent : 


DüuRÉE DE LA CHUTE 


COURANT Es Re 
Géminides 0,143 sC.107 
Taurides Encke 0,605 56.107 
Bielides 1852 III 2 TOMES UE 
Giacobinides 1933 III SOS CT 
Orionides Halley EA0MSS TO 
Léonides 1866 I T2 RRES CO 
Perséides 1862 III 11900 ES TO 
Lyrides 1861 I 11860 CAO! 


D’après ces résultats une météorite de densité 4 et 
de rayon 0,04 cm (pouvant donner naissance à un 
météore de 5° magnitude) se mouvant sur l’orbite de 


la comète de Halley, tomberait en spirale sur le Soleil 


en dix millions d’années. Toutefois les nombreuses 
perturbations auxquelles sont soumises les météo- 


rites de la part des planètes peuvent modifier notable- 


ment leurs éléments dans un délai beaucoup plus court. 
. Revenons à la variation de a et e en fonction du 
temps. L’équation (45) s’écrit : 
esls 


æ 
LL 


G{e) = Geo) 


LOT EN PARC RS 
AE &uA)l1 réa) SR CELUE 


et permet en utilisant la table de G(e) de calculer e pour 
chaque valeur de 


_t—t 
 stC107 


a se déduit ensuite des équations (41) et (42). Effee- 
tuant les calculs pour divers courants, Wyatt et 
Whipple trouvent une décroissance presque linéaire 
de a jusqu’à ce que cette quantité atteigne une valeur 
voisine d’une unité astronomique puis une décrois- 
sance plus rapide ensuite. Par contre, la distance 
périhélie q varie peu. Les courbes de la figure 3 repré- 
sentent ces résultats pour les Léonides (ar, gr) et pour 
les géminides (ag, ge). Par ailleurs les formules (38) 


ibn Edfe ic ct 


MSN 


Re en ES LR NE RL 5 ut 
. < ? 


N°4 


montrent que les variations de a et e et par suite de g 
sont proportionnelles à o, donc inversement propor- 
tionnelles à s et € d’après (18). Il en résulte qu’au bout 
d’un temps suffisamment long, alors que la distance 


a ou q (UA) 


CPE PT ENT 8 
(t _ € )années 
Use 
TGS: 


périhélie aura peu varié pour toutes les particules le 
demi-grand axe aura diminué davantage pour les 
petites particules que pour les grosses. Iltend à s’établir 
la disposition indiquée par la figure 4. Lors de la tra- 


nee 


versée par la Terre d’une partie du courant (sui- 
vant AB), on rencontrera d’abord de grosses parti- 
cules, ensuite de plus petites. Il y aura, en général, une 
corrélation entre le temps écoulé depuis le début de la 
traversée du courant et la grosseur, et par suite l’éclat, 
des météores. Des recherches effectuées dans le but de 
mettre en évidence un tel phénomène n’ont pas conduit 
à des résultats probants. Wyatt et Whipple en con- 


æluent que les courants météoriques sont d’origine 


suffisamment récente pour que la séparation des par- 


. l'équation : 
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ticules n’ait pas eu le temps de se produire. On obtien- 
drait d’après ces auteurs un effet sensible pour les Gémi- 
nides si leur âge atteignait 7.104 ans et pour les 
Taurides si leur âge atteignait 5.105 ans. 

Considérant les météores sporadiques, Wyatt et 
Whipple examinent l'hypothèse d’après laquelle une 
partie de ces corps proviendrait d’une planète gravi- 
tant entre Mars et Jupiter (à une distance du Soleil 
de 3 unités astronomiques environ) et qui aurait éclaté. 
On peut calculer le temps nécessaire à la chute des 
fragments sur le Soleil en attribuant à leurs orbites 
des demi-grands axes inférieurs à celui de l’orbite 
de Jupiter, et ceci pour diverses valeurs de l’excen- 
tricité &. La densité des particules est supposée de 
4 grammes par centimètre cube. Les résultats sont 
donnés dans le tableau ci-dessous. 


TEMPS DE CHUTE POUR DES PARTICULES ASTÉROIDALES 


CARRE THEN (distance aphélie) LiAntets! 
5 OMRNES 5 70.107s (s en em) 
3 0 3 3 25 t00S 
4 02255 5 11=10%S 
2 0,50. 3 2.107s 
3 07071 5 SOS 


Considérons des particules de 4 centimètres de rayon. 
Dans le cas où la chute est la plus lente, ils seront 
tombés sur le soleil Soleil en 


70 10742 8 410%ans 


et des particules plus petites tomberaient en un temps 
moindre. Si, donc, on fait remonter la rupture de la 
planète à 3.10° ans, on ne peut expliquer par cette 
hypothèse que l’existence de météores très lumineux. 
Si l’on suppose avec Bauer [7] une rupture remontant à 
60 millions d’années, la table précédente montre que les 
particules de rayon moindre que : 


8,5 10—? centimètre 


seraient déjà tombées sur le Soleil. En effet, dans les 
conditions de chute les plus lentes, une particule qui 
tombe sur le Soleil en 6.107 ans a un rayon donné par 


CAO 0 100 
S5—18,0 1l0=-CIn: 


Les auteurs, considérant qu'un météorite de 
5.10-2 em de rayon cause un météore de cinquième 
magnitude environ, en concluent que seuls les météores 
de luminosité supérieure à la cinquième magnitude 
peuvent provenir d’une planète éclatée. 

L'effet Poynting-Robertson provoque la chute en 
spirale sur le Soleil des particules qui constituent la 
couronne F et le nuage causant la lumière zodiacale. 
Il en résulte la nécessité d’un apport extérieur de parti- 
cules si le nuage doit se maintenir. Avec les hypothèses 
admises sur l’absorption intégrale du rayonnement 
incident et en admettant une rémission intégrale 
ultérieure, connaissant en outre le rayonnement éner- 
gétique du Soleil, il est possible de calculer la quantité 
totale d'énergie rayonnée au cours de sa chute par une 
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particule de masse m. On en déduit l’énergie rayonnée 
par unité de temps. D'autre part le rayonnement 
global de la lumière zodiacale et de la couronne F se 
monte à 3.10-7 de la quantité d’énergie rayonnée par 
le Soleil. On peut en déduire [5] l’apport nécessaire 
au maintien du nuage dans son éclat actuel, soit : 


une tonne par seconde. 


Ce résultat a été obtenu en tenant compte unique- 
ment du rayonnement électro-magnétique du Soleil. 
Si l’on tient compte de l’effet Poynting-Robertson dû 
au rayonnement corpusculaire et aussi de la diminu- 
tion de masse des particules correspondant aux atomes 
arrachés par les corpuscules incidents (ions et atomes 
d'hydrogène) on trouve que la quantité de matière à 
fournir au nuage pour assurer son maintien est for- 
tement augmentée ; mais les connaissances actuelles 
sur le rayonnement corpusculaire du Soleil ne per- 
mettent pas d’en tenir compte d’une manière sûre. 

D’après Whipple il faut rechercher l’origine de cet 
apport dans les éjections nucléaires des comètes. Ses 
recherches sur la quantité de matière émise par un 
noyau comètaire d’une part, et les statistiques sur le 
nombre des comètes d’autre part l’ont conduit à éva- 
luer à trente tonnes par seconde la quantité totale 
de matière émise par les comètes sous forme de parti- 
cules et à cent tonnes par seconde leur émission ga- 
zeuse. Les particules chassées des noyaux cométaires 
décrivent des orbites peu différentes de celles des 
comètes parentes. Elles tombent ensuite en spirale 
sur le Soleil, mais subissent au cours de cette chute des 
influences qui peuvent les arrêter ou les détruire. 
Ce sont : 

10 les perturbations des planètes, en particulier de 
Jupiter qui tendent à communiquer aux particules une 
distance aphélie voisine de la distance moyenne de 
cette planète, 

20 La possibilité de rencontres destructrices entre 
particules, 

30 L’existence possible d’un courant de matière 
(atomes surtout d'hydrogène) qui traverserait le Sys- 
tème Solaire. 

Whipple analyse ces influences et conclut qu’une 
partie notable des particules peut atteindre les régions 
voisines du Soleil (c’est-à-dire les régions situées à 
moins de 1,66 unité astronomique de cet astre). Les 
résultats conduisent en outre à une prédominance des 
mouvements directs sur les mouvements rétrogrades 
et à une répartition des rayons des particules qui 
montre une brusque chute de fréquence pour des rayons 
supérieurs à 10? centimètre. Il y a sur ce point accord 
avec les résultats de Van de Hulst d’après lesquels les 
-particules du nuage zodiacal montrent une chute de 
fréquence pour des rayons supérieurs à 0,035 centi- 
mètre. 


L'effet Poynting-Robertson permettra peut-être de 
résoudre une difficulté qui a surgi à propos de la den- 
sité y du nuage zodiacal au voisinage de l’orbite terres- 
tre. Utilisant les mesures de luminosité de la lumière 
zodiacale et de la couronne F et moyennant certaines 
hypothèses sur la fonction de fréquence des rayons des 
particules divers auteurs sont parvenus aux résultats 
suivants : 


N° #4 


Allen (1946) u "4 ,1072%81g/cm* 
Fessenkov (1947) u'="64100% 2eme 
Van de Hulst (1947) 9 100) 
Behret Siedentopf (1953) u —  10—?3 g/cm? 


Le travail de Van de Hulst a été critiqué par 
Lévine [8] qui lui reproche d’avoir attribué aux parti- 
cules un albédo «incroyablement faible » (< 0,005) et 
d’avoir utilisé une représentation de la brillance de la 
lumière zodiacale en fonction de l’élongation du Soleil 
qui est hypothétique puisqu'elle résulte d’une inter- 
polation entre des observations intéressant de petites 
élongations moindres que cinq degrés (Baumbach et 
Grotrian) et des observations concernant des élon- 
gations supérieures à trente degrés (Hoffmeister, 
Élvey et Roach). En outre les travaux de Van de 
Hulst sont fondés sur la réflexion et la diffusion de la 
lumière par des particules alors que ceux de Behr et 
Siedentopf ont mis en évidence le rôle fondamental 
joué par les électrons dans la partie interne du nuage 
zodiacal (à moins de 0,6 unité astronomique). 

La densité de la matière au voisinage de l’orbite 
terrestre peut aussi se déduire des comptages de 
météores, de la théorie de leur luminosité, et de la 
répartition de leurs vitesses géocentriques. On trouve : 


LME p)/CmS 


en accord avec les résultats de Fessenkov, Allen, 
Behr et Siedentopf [8]. 

Le résultat discordant de Van de Hulst trouve cepen- 
dans une confirmation dans la détermination de la 
masse de matière reçue par la Terre, détermination 
fondée sur l’étude des limons océaniques. Déjà, en 
1876 Murray et Renard [9] avaient signalé des sphé- 
rules magnétiques de 0,3 millimètres environ de dia- 
mètre qui avaient été draguées au fond des mers. 
Ces sphérules consistaient, comme les fers météo- 
riques en fer allié au nickel et au cobalt. De nom- 
breuses » carottes « de sédiments ont été prélevées ces 
dernières années dans diverses mers profondes. Les 
échantillons soumis à un puissant extracteur ma- 
gnétique ont donné un nombre de sphérules dix 
fois plus grand qu'il n’avait été signalé par Murray et 
Renard. Mais il y a mieux : en déterminant la teneur 
en nickel de ces sédiments et en admettant que ce 
nickel est d’origine cosmique, utilisant d’autre part 
la proportion de nickel contenue dans les météorites, 
Petterson et Rotchi en ont déduit la quantité de fer 
météorique tombant sur la Terre. La proportion des 
fers parmi les météorites étant connue, ces auteurs 
peuvent évaluer à 4 000 à 5 000 tonnes la quantité 
totale de matière météorique tombant chaque jour sur 
la Terre. La collecte de férules magnétiques sur les 
montagnes et dans les régions polaires conduit à des 
résultats de même ordre. On en déduit une valeur de la 
densité de matière au voisinage de la Terre 


u = 102 gramme par centimètre-cube. 


“ Ainsi, même si l’on prend en considération les 
remarques de Levine au sujet du travail de Van de 
Hulst, il semble que la densité de la matière au voi- 
sinage de l'orbite terrestre soit mille fois plus grande 
que ne le font prévoir les comptages de météores 
L’explication ‘de ce’résultat peut se déduire des idée 


No 


actuelles sur les orbites des corps météoriques. D’après 
Opik [10] une partie très importante de ces Corps, 
peut-être 90 % décrivent autour du Soleil des trajec- 
toires peu excentriques et de faible inclinaison ; une 
partie beaucoup moins importante se meut sur des 
orbites de grande excentricité et provient de la désin- 
grégation des comètes ; enfin une faible partie serait 
hyperbolique. Parmi les corps du premier groupe dont 
les vitesses par rapport à la Terre sont relativement 
faibles, un grand nombre ne pourraient de ce fait 
donner un météore lumineux, surtout les plus petits de 
diamètre voisin de dix microns (micrométéorites). Ces 
très petits corps sont fortement influencés par l’effet 


. Poynting-Robertson et tombent*en réalité en spirale 


sur le Soleil. De Jager [11] a montré qu’un grand 
nombre d’entre eux doivent être captés au passage par 
la Terre. Soit un nuage de particules tournant autour 
du Soleil sur des orbites circulaires avec des périodes 
de révolution obéissant à la 3e loi de Kepler. 
Les vitesses de ces corps (et celle de la Terre) 
décroissent quand la distance du Soleil augmente, ceux 
qui sont plus éloignés du Soleil que la Terre ont une 
vitesse moindre que celle-ci, ceux qui sont plus rap- 
prochés une vitesse plus grande. De Jager montre que 
les particules qui passent à une distance de notre pla- 
nète moindre que 


4,5 106 kilomètres 


décrivent par rapport à celle-ci des ellipses et 
deviennent des satellites de la Terre. Mais à ce moment 
(fig. 5) l'effet Poynting-Robertson introduit une force 
jf constamment opposée à la vitesse héliocentrique de 
la particule P, c’est-à-dire perpendiculaire à la direction 
du Soleil. Cette force tend à augmenter continuellement 
l’excentricité de l’orbite par rapport à la Terre, et cette 
orbite s’allongeant de plus en plus finit par rencontrer 
la surface de notre planète. Ainsi, toute particule 
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captée est destinée à tomber sur la Terre. L’effet 
Poynting-Robertson joue ici un rôle double : 

19 Il fait tomber les particules en spirale sur le 
Soleil et les amène ainsi à passer à moins de 4.5.105 
kilomètres de la Terre ce qui entraîne leur capture. 


Prices 


20 Au cours de leur mouvement autour de la Terre 
il allonge leur orbite jusqu’à amener une rencontre 
avec la surface de notre planète. 

Ces particules animées de vitesses géocentriques 
trop faibles pour donner naissance à des météores 
lumineux et d’ailleurs trop petites augmentent cepen- 
dant la quantité de matière météorique reçue par la 
Terre. 

Manuscrit reçu le 10 novembre 1958. 
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LETTRES À LA RÉDACTION 


STRUCTURE ET POUVOIR RÉFLECTEUR 
DES COUCHES DE SILICIUM 
OBTENUES PAR ÉVAPORATION SOUS VIDE 


Par Mme Sonja RoBiN-KaNDaRE et M. H. Damanvy, 
Laboratoire des Hautes Pressions, 


et Mne L. TERTIAN, 
Laboratoire des Rayons X, C. N. R. $., Bellevue. 


Le silicium présente dans l’ultra-violet lointain un 
très bon pouvoir réflecteur et il nous a paru nécessaire 
de compléter nos expériences préliminaires [1] par 
l’étude de couches préparées dans des conditions sensi- 
blement améliorées. 


I. Préparation des couches. -— Toutes les couches 
ont été préparées dans un vide dynamique d’envi- 
ron 3.10-5 mm de mercure, créé par une pompe à 
diffusion d'huile de grand débit. Nous avons utilisé 
trois méthodes d’évaporation thermique : a) Emploi 
d’une gouttière de tantale chauffée par passage de 
courant. b) Emploi d’un creuset de graphite bombardé 
par des électrons. c) Emploi d’un canon à électrons [2]. 

Les deux premières méthodes présentent de gros 
inconvénients. Pour a), le tantale est très fortement 
attaqué par le silicium : il se forme un alliage qui pro- 
voque rapidement la rupture de la gouttière. De plus, 
les couches ainsi préparées contiennent certainement 
des traces non négligeables de tantale. Pour b), le 
silicium forme avec le graphite du carbure de silicium 
et le plus souvent la couche de carbure ne protège pas 
le graphite d’une attaque ultérieure : le creuset se 
fendille et ne garde pas le silicium fondu. Aussi est-il 
difficile d’obtenir des couches épaisses par ce procédé. 

La méthode c) nous a permis la préparation de 
couches épaisses, bien adhérentes et très pures car, 
convenablement appliquée, elle élimine tout problème 
de support. Le dispositif est le suivant : un morceau 
d’un monocristal de silicium très pur repose sur un 


bloc de graphite, en bon contact thermique et élec-. 


trique avec ce dernier. On bombarde la partie supé- 
rieure du cristal avec ‘un faisceau d’électrons fourni 
par un canon ; le silicium fond partiellement et c’est le 
reste du cristal qui sert de creuset ; le bloc de graphite 
joue le rôle de volant thermique et évite la fusion 
totale du silicium. (fig. 1). 

Nous avons évaporé Si sur trois sortes de supports ; 
verre, quartz cristallisé, corindon. Le dépôt se faisait 
sur une face polie plan-optique, soigneusement 
nettoyée et dégraissée avant introduction dans la 
cloche, puis soumise à un bombardement ionique de 
20 minutes. 

L’épaisseur de nos couches n’était pas rigoureu- 
sement uniforme, car l’évaporation donne une symétrie 


semi-circulaire sur une plaque de grandes dimensions ; 
mais les plaquettes présentaient des dimensions assez 
faibles et la différence d'épaisseur n’était pas tres sen- 
sible. Nous avions la possibilité d'améliorer l’uniformité 
en faisant tourner la plaquette au-dessus de l’évapo- 
rateur. L’épaisseur des différentes couches étudiées 
était comprise entre 3 000 et 15 000 À (mesurée par 
interférométrie), les temps d’évaporation variaient de 
40 à 90 minutes. Une évaporation aussi longue élève la 


Basse 
tension 
A 


Canon 


Si fondu 


Si solide 


Graphite | 
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+ 15 kV 


température du support jusqu’au voisinage de 100 oC. 
Pour éviter cet effet, nous avons dans certains cas, 
fractionné l’évaporation : par exemple, évaporation 
pendant 10 minutes, puis arrêt pendant 20 minutes, etc. 

Nous avons soumis nos couches à trois types de 
traitement thermique sous vide : I) chauffage pendant 
l’évaporation ; II) chauffage après l’évaporation sans 
que la couche soit exposée à l’air ; III) chauffage sous 
vide après exposition à l’air. 

Let IT sont obtenus par contact du support avec un 
ruban de tantale dont la température est indiquée par 
un thermocouple chromel-alumel. Pour III, nous utili- 
sions un four constitué par un bloc d’acier dans lequel 
est incorporé un enroulement de tungstène, et la tem- 
pérature est indiquée par un thermocouple chromel- 
alumel fixé sur le bloc. 

Les couches d'épaisseur non uniforme ont présenté 
souvent une surface craquelée (fig. 2, A) ; ces craque- 
lures, déjà amorcées avant chauffage, s’accentuant 


No4 


sous l'effet, de celui-ci. 1] semble que ce phénomène de 
fracture n intéresse que les couches de forte épaisseur 
et soit lié à l’existence de tensions dans les couches [3]. 


IRTCHNCE 


Il n’affecte pas l’adhérence car il n’y a pas pelliculage ; 
de plus après dissolution de la couche, on constate 
que le quartz ou le verre sous-jacent est lui-même 
craquelé sur une faible profondeur (fig. 2, B). 


IT. Structure des couches. — On sait que les couches 
évaporées sur support non chauffé sont amorphes. 
Nous avons examiné, par diffraction électronique en 


réflexion, l’influence sur la structure des trois types de 


traitement thermique que nous venons de décrire. 
Remarquons tout de suite que ce que l’on observe par 
cette méthode est la structure d’une tranche super- 
ficielle de quelques dizaines d’angstrôoms seulement, 
et que l’on n’obtient pas ainsi directement des rensei- 
gnements sur l’intérieur de la couche. Le traitement 
thermique 1, effectué à plus de 600 °C donne des 
couches cristallisées, avec légère orientation des cris- 
tallites : axe [110] perpendiculaire au plan du support ; 
les couches présentaient alors sous certaines angles un 
aspect laiteux. Le traitement thermique IT effectué 
à 700 °C environ, donne un début très net de cristal- 
lisation et il est très vraisemblable qu'avec une tempé- 
rature un peu plus élevée, on observerait une cristal- 
lisation complète. Par contre le traitement III, effectué 
à différentes températures entre 720 et 780 °C pendant 
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au moins deux heures, ne donne pas de cristallisation. 
Nous avons alors porté la température à 970 °C pendant 
deux heures trente, toujours sans résultat ; on cons- 
tate seulement une diminution de l’épaisseur par rééva- 
poration. Ces dernières expériences semblent en con- 
tradiction avec les résultats de différents auteurs [4], [5] 
mais un nouvel examen par diffraction en réflexion, 
après un décapage chimique qui enlève une partie 
importante de la couche, révèle que l’intérieur de la 
couche est effectivement cristallisé. Remarquons que 
les auteurs cités examinaient par transmission des 
couches relativement minces et sans support. 


ITT. Pouvoir réflecteur dans l’ultra-violet proche et loin- 
tain. — Les mesures du pouvoir réflecteur ont été effec- 
tuées dans un monochromateur à vide décrit par ailleurs 
[6]. Le domaine spectral s'étend de 1 000 à 3 000 À grâce 
à l’emploi d’une lampe à hydrogène sans fenêtre. Une 
étude préliminaire [1] sur des couches minces de sili- 
cium évaporées à partir d’un ruban de Ta dans un 
vide moyen de quelques 10-5 mm de mercure donnait, 
dans une région spectrale limitée à 1 000-1 600 À, des 
résultats en accord avec les mesures présentées ici. 
Il nous a paru nécessaire de vérifier ces résultats en 
éliminant l'influence éventuelle des traces de Ta dans 
les couches. L’emploi du canon à électrons nous l’a 
permis. 

La valeur du pouvoir réflecteur À des couches 
amorphes dépend beaucoup des conditions d’évapo- 
ration, alors que la forme de la courbe À — (à) reste 
sensiblement la même. On constate en particulier 
qu’une évaporation prolongée et continue donne des 
résultats inférieurs à ceux que l’on obtient par une 
évaporation intermittente. Les courbes A et B de la 
figure 3 illustrent la différence entre ceux deux cas. 


A. Évaporation prolongée 
B. » 


GC. Couche cristallisée 


intermittente 


La courbe C représente le pouvoir réflecteur d’une 
couche évaporée sur un support de quartz chauffé 
à 700 oC pendant l’évaporation ; la couche présentait 
un aspect laiteux. L’examen par diffraction électro- 
nique montre une cristallisation avec orientation préfé- 
rentielle des cristallites (axe [110] perpendiculaire au 
plan du support) ; pour pouvoir effectuer la mesure 
optique, il a fallu procéder à un léger polissage méca- 
nique (alumine sur feutre) qui rend amorphe la partie 
superficielle (un décapage chimique révèle que l’inté- 
rieur de la couche reste cristallisé). Le pouvoir réflec- 
teur des couches cristallisées est plus élevé que celui 
des couches amorphes et présente deux maxima : l’un 
faible à 5,3 eV, l’autre plus fort à 4,5 eV. Les discon- 
tinuités du pouvoir réflecteur peuvent être reliées à des 
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transitions électroniques mais il semble difficile dans 
l’état actuel de la théorie des bandes de donner un 
calcul quantitatif satisfaisant. 

Lettre reçue le 9 février 1959. 
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VARIATION 
DE LA RÉSISTIVITÉ DE COUCHES MINCES 
DE Ge ET Te PAR LA PRESSION 


Par Jean RoBin, 
Laboratoire des Hautes Pressions, Bellevue (Seine-et-Oise). 


On sait que, à température constante, la pression a 
une influence sur la résistivité intrinsèque de différents 
corps semi-conducteurs (voir par exemple [1], [2], [31). 
Cette influence se traduit surtout par une variation 
(dE, 1dP)r de la largeur Æ, de la bande interdite avec 
la pression. La plupart des travaux effectués jusqu’à 
maintenant dans ce domaine portent soit sur des 
poudres. soit surtout sur des échantillons massifs mono- 
cristallins. 

Le présent travail rapporte quelques résultats préli- 
minaires que nous avons obtenus en soumettant des 
couches minces de germanium et de tellure à une 


Frétds 


pression de gaz neutre. Pour cela nous avons construit 
un récipient résistant à la pression et possédant une 
embase, représentée figure 1, comportant 12 entrées de 
courant de type conique. L’isolement par du téflon de 


N° 4 
ces entrées, représenté en figure 2, varie de 1011 Q 
jusqu’à plus de 5.1015 Q (limite des mesures possibles 
avec l’Iso-R-mètre Lemouzy utilisé). 


Intérieur 
de Ja 
cuve 


Embase d'acier 


Passage de courant 
Fc. 2: 


Nous avons opéré sur des couches minces évaporées 
assez rapidement sous vide sur des supports de verre 
à extrémités métallisées. Les épaisseurs des couches 
ont été déterminées par méthode interférentielle. Il ne 
semble pas que, jusqu’à maintenant, des travaux de ce 
genre aient fait l’objet de publications importantes. 


a) Germanium. — Nous avons utilisé du germanium 
«spectroscopiquement pur » dont la résistivité p était 
de 30 Q cm. On trouve que la résistivité des couches 
évaporées est de 3 à 5.104Q.cm. Or, pour des couches 
de germanium amorphe on devrait s’attendre [4] à une 
résistivité beaucoup plus grande (de l’ordre de 108 à 
1011 Q.cm). On peut donc penser que la structure 
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de nos couches pourrait être très microcristalline. 
Ê épaisseur des ue varie entre 1 100 et 1 600 À. 

OuS avons mesuré la résistance de ces couches jusqu’à 
2 000 kg/cm? à 18 oC. La figure 3 représente pue a 
fonction de la pression. La courbe en trait plein corres- 
pond aux montées en pression et celle en pointillé aux 
pressions décroissantes. Les deux branches de courbe 


L 


ETS 


N°4 


sont parallèles ; l0g:, p croit linéairement avec la pres- 
sion. À partir de ces valeurs de log, p on calcule que : 


DES IdDP)r = 8.10—$6 eV.cm!?/kg. 


La résistance de la couche après mise sous pression 
est environ 5 % plus élevée qu’elle ne l’était avant, 
probablement par suite d’une légère modification de 
structure. 

Les valeurs de (d0£,/dP)» obtenues sur des échan- 
tillons monocristallins, généralement du type n (voir 
par exemple [3], [5], [6]) sont de l’ordre de 5,5 10-5 ce 
qui est de même signe et du même ordre de grandeur 
que les valeurs que nous obtenons. 


b) Tellure. — Le tellure utilisé était également de 
qualité « spectroscopiquement pure » et était de plus 
purifié au Laboratoire par distillations sous vide. Les 
couches de tellure étaient cristallisées et leurs épais- 
seurs variables de 400 à 1 200 À. La résistivité de ces 
couches variait entre 0,112 et 0,120 Q.em. 

Les mesures ont, comme pour le germanium, été 


faites jusqu’à 2 000 kg/cm? et ramenées à 18 °C par 
correction de l’effet de température. La figure 4 donne, 
pour une seule montée en pression, log: p en fonction 
de P pour deux des échantillons étudiés ; plusieurs 
mesures effectuées à même pression et température 
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montrent que p varie légèrement avec le temps, avec 
tendance à une diminution de l’hystérésis. Les résul- 
tats obtenus sont les mêmes et leur reproductibilité est 
également très bonne pour plusieurs montées et des- 
centes en pression. Chaque montée en pression (courbe 
en trait plein) s'effectue en une journée et chaque 
descente (en pointillé) en une quinzaine de jours. Étant 
donné la légère évolution de p avec le temps, les 
courbes en pointillé sont certainement plus proches de 
l’état d'équilibre d’adsorption ; on peut même sup- 
poser que si les montées et descentes en pression étaient 
faites en des temps beaucoup plus longs le phénomène 
d’hystérésis ne serait plus décelable. 

On constate que la courbe donnant log,, p décroît 
lorsque la pression augmente et présente une légère 
concavité vers le haut. 

De nos courbes expérimentales il résulte que la 
valeur moyenne de (04, /dP)r entre 1 et 2 000 kg/cem? 
est : 

(DE, [dP)}r = —11.10—$ eV .cm?/kg. 


Ce coefficient serait un peu plus faible vers les hautes 
pressions et un peu plus fort vers les faibles pressions. 

Les valeurs trouvées avec des échantillons mono- 
cristallins taillés soit parallèlement, soit perpendi- 
culairement à l’axe C sont respectivement de l’ordre 
de — 19.106 et — 15.106 eV.cm?/kg (voir par 
exemple [7] et [8]). 

Les variations de la largeur de la bande interdite Æ, 
du germanium, et du tellure obtenues par nos travaux 
sur des couches minces polycristallines (et même très 
microcristallincs) se rapprochent beaucoup de celles 
obtenues par les auteurs antérieurs sur des échantillons 
massifs monocristallins. Evidemment l'influence du 
support, tant par sa nature que sa compressibilité et 
l’adhésion de la couche sur lui n’est pas connue et il ne 
peut en être tenu compte dans ce travail préliminaire. 
Il en est de même pour la variation avec la pression 
des barrières de potentiel entre les microcristaux. 

Cependant ces résultats sont suffisamment encou- 
rageants pour qu'on puisse espérer, avec des couches 
plus parfaites, et notamment monocristallines, pouvoir 
tirer des conclusions plus solides à partir de mesures de 
résistivité sous pression, en particulier sur les pro- 
priétés de surface des semi-conducteurs et leur évolu- 
tion, qui sont difficiles à obtenir par des études effec- 
tuées sur des monocristaux massifs. 


Lettre reçue le 9 février 1959. 
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COEFFICIENT DE CONVERSION INTERNE 
DU RAYONNEMENT DE 92 keV DU **“Pa 


Par R. FoucHER, 
Laboratoires de Physique Nucléaire d'Orsay (Seine-et-Oise). 


La mise en évidence d’effets de la structure nucléaire 
sur les coefficients de conversion interne [1], [2] peut 
être favorisée dans les noyaux lourds impair-impair 
déformés possédant des rayonnements M, interdits 


28Th 


2341 


Hire 


1e 2e 


dont l’énergie est inférieure à 150 keV pour lesquels un 
mélange M; + E, ne peut expliquer une conversion 
trop faible. 

Dans cette optique, nous avons repris en spectro- 
métrie y l’étude du ?*4Pa (UX,). Les premières esti- 
mations du coefficient B7 de conversion L apparaissent 
effectivement faibles par rapport aux prévisions de 
Shv et Rose [3]. 

Du schéma de désintégration de l’'UX, (fig. À) est 
tirée l’équation (1) donnant $Z en fonction du nombre 
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de photons X7, 63 keV, 92 keV émis dans l’angle 
solide utilisé et déduits de l’analyse du spectre (fig. 2) 


L | or NET) 
B1(92) — Ur) N02) 
wz(E,) N(63) ; ar, 29 , 
— D Non [269 + EU + x(631] (1) 


az(29), «(29), «r(63), «(63) sont les coefficients de 


conversion L et L + M + ... des rayonnements 

électriques de 29 et 63 keV [4], [5], [6], [7]: © (M), 

w (E,) sont les coefficients de fluorescence globaux 

des couches L excités par conversion interne M, et E:. 

Il faudra évidemment déduire du spectre figure 2, les. 
photons Xr, Xx, Compton, de freinage, de la tran- 

sition UX,-2%4U. 


Conditions expérimentales. — Les sources de 1 à 6uC 
ont été extraites d’un minerai nouvellement traité par 
MM. Muxart et Conte, de l’Institut du Radium. Les 
détecteurs y sont des cristaux de INa(T1) dont l’un 
a 35 mm de diamètre, 7 mm de hauteur et est recou- 
vert d’uné feuille d'aluminium de 20 x rendue diffu- 
sante par 500 ug/em? de magnésie. Les photomulti- 
plicateurs dont des EMI 5914 — les analyseurs et 
circuits de coïncidences seront publiés par ailleurs. 
Nous avons utilisé également un analyseur RIDL à. 
256 canaux et un analyseur 100 canaux du même type ; 
afin d’éviter les additions X,-63 les angles solides sont 
de 2 à 5 %. Pour arrêter les 6 de 2 MeV de l'UX, des 
écrans de 5 ou 7 mm de bérylium sont interposés. 


Étalonnage du spectre. Corrections. — L’étalonnage 
en énergie a été réalisé avec des raies de 18 ; 22 ; 24; 
34; 47; 88,5 ; 140 keV. Les largeurs à mi-hauteur de 
ces raies suivent très bien la relation : 


b 
Lin = c+ Paie L'—15,5 % pour £ = 88,5 keV. 


Les largeurs relatives des raies complexes Xz 
du Ra D et ?*#Pa s’écartent très nettement de cette 
relation. Leur forme a été déterminée avec précision ; 
elle est donnée par l’addition des raies Lx, LB, Ly de 
fluorescence L. En admettant qu’elle est formée pour 
l'essentiel de 50 % de photon d'énergie 13 keV et 50 % 
d'énergie 17 keV, on peut calculer l’absorption dans 


les écrans. Pour le bérylium on trouve a — 7 — 450: 


l'expérience donne 1,54. Le pourcentage des raies La, 
LB, Ly a été pris d’après les résultats de Hoff-Olsen- 
Mann sur les X, émis lors de la capture L du 2#5Np [8], 
En tenant compte des autres facteurs d'absorption, 
d'efficacité, d’addition X,-63 on trouve : 


a(Xr) = 1,72 a(63) — 1,22  a(92) — 1,27. 


L’absorption dans la source qui n’était pas tota- 
lement dépourvue de matière a été étudiée en mesurant 
des sources d'intensité et de dimension variées. 

Les rapports pic d’échappement/pic photoélectrique 
ont été déterminés expérimentalement sur les Y 
de PUX, et celui de 88,5 keV du Te. 


R(63) —12à14 %  R(92) = 4 à 6 4. 


_N°& 


Enfin les proportions de raies X, et Xx (2#{U) con- 
tenues dans le spectre ont été déterminées par coïnci- 
dences £-spectre et X;-spectre. Nous trouvons 
N(Xr) = N(Xi) = 15% du 92 keV. On voit sur 
l'équation 1 qu’une erreur de 30 %, sur ce chiffre fait 
varier f} de quelques pour cent. 


Valeurs de &z ; &w utilisées. — Pour «z(29), les 
valeurs 10 ou 11 ont été proposées [4], [6]. Nous avons 
remesuré ce coefficient en effectuant les coïncidences 
63-29 keV,63-X , nous trouvons «y — 60. Nous revien- 
drons ailleurs sur cette divergence qui fait peu varier 
le rapport 


az (29) : 
—— — () 9; F5 
1 + a(29) D,95 + 0,05. 


Pour «7(63) les mêmes auteurs trouvent 0,45 et 
0,42 ; notre mesure qui peut également être sures- 
timée donne 0,39 (pour 47 — 0,46). Nous prendrons 
donc la valeur théorique de Rose et Sliv pour un £; 
pur : «y — 0,31 — alors «(63) sera voisin de 0,40. 
Enfin pour & , qui est la grandeur la plus importante, 
nous prendrons une limite inférieure. D’après les 
valeurs de Kinsey [9] qui ne sont pas en contradiction 
avec les résultats expérimentaux récents [10], on 
aurait 7, — 0,41 pour une excitation M,, o7 — 0,46 
pour une excitation Æ, ; en prenant 


OL(E) — &@z(M;) — 0,40 


. nous sommes certains de ne pas minorer B£. 


Analyse du spectre. — Les formes des raies Xz, 
63, 92 keV sont bien établies. L’imprécision de l’ana- 
lyse provient de l’évaluation du fond continu. Une 
mesure effectuée en arrêtant les rayonnements directs 
d'énergie inférieure à 100 keV donne une limite infé- 
rieure, La connaissance de la forme et de la largeur à 
mi-hauteur des raies permet d'évaluer un maximum 
de ce fond. Dans ces conditions pour 5 mesures nous 
obtenons pour le premier terme de l’équation (1) les 
valeurs limites 3,85 et 4,20, pour le deuxième 
terme 1,36 et 1,53, donc : 


2,3 < BE < 2,9. 


Étant donné les valeurs prises pour &,, le sens des 
corrections, et comme nous ne trouvons pas de déficit 
de photons X dans les mesures parallèles de «(29) 
et «(63) [11] au contraire — la valeur maximum 
obtenue nous semble être bien une limite. Ce résultat 
n’est pas en contradiction avec les estimations anté- 
rieures qui sont de l’ordre ou inférieures à 2, [4],[6],[7]. 
Une détermination plus directe de $£ est évidemment 
souhaitable. Une mesure du rapport Li/Lrr aix est 
actuellement en cours. 


COMPARAISON AVEC LES RÉSULTATS THÉORIQUES 


Rose SLIV Exr. 


gl &,7 ‘4,2 2,3 à 2,9 


Ainsi, le coefficient de conversion expérimental de 
cette transition M, interdite du noyau impair-impair 
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déformé 234Pa serait nettement inférieur aux valeurs 
théoriques, ce qui indiquerait un important effet de la 
structure nucléaire sur ce coefficient de conversion. 


Lettre reçue le 25 février 1959. 
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ANALYSE PAR SPECTROCHIMIE D'’ÉMISSION 
DANS L’ULTRAVIOLET LOINTAIN 
DE POUDRES CONTENANT 
LES ÉLÉMENTS SUIVANTS : 
GERMANIUM, ARSENIC, SÉLÉNIUM, BROME, 
ÉTAIN, ANTIMOINE, TELLURE. 


Par Mlle Germaine BALLOFFET, 
M. Jacques Romanp et M. Boris Vopar, 


Laboratoire des Hautes Pressions, Bellevue (Seine-et-Oise). 


On a déjà parlé [1] de l’application de l’étincelle 
glissante dans le vide au dosage de poudres contenant 
des éléments difficiles à analyser par spectrographie 
classique. Des dosages du soufre et du phosphore en 
particulier dans des poudres à base d’alumine ont 
donné une sensibilité meilleure que dans le cas de 
l’analyse de ces deux éléments dans les aciers. Il nous a 
paru intéressant de vérifier que les raies «ultimes » [2] 
des éléments ayant un nombre d’électrons optiques 
variant de 4 à 7 sur la 4€ et la 5€ couche permettaient 
bien de tracer des courbes de dosage. Cette étude a été 
faite à partir de poudres de concentrations connues 
fabriquées par synthèse. En même temps il était pos- 
sible de préciser les teneurs limite décelables dans les 
mêmes conditions d’excitation et en présence des 
mêmes matériaux de base pour les divers éléments 
étudiés. 

Rappelons tout d’abord les transitions susceptibles 
de donner les raies les plus sensibles pour ces éléments. 
On a indiqué dans le tableau n° 1 le doublet énus par 
l'ion ne possédant plus qu’un électron optique dans 
l’état de base s lorsque cet électron passe à l’état 


r(5—3) ainsi que la raie émise par l'ion ayant 


2 électrons s dont l’un passe au niveau p. 
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N° 4 


TABLEAU 1 


LONGUEURS D’ONDE DES RAIES SENSIBLES 


ÉLÉMENT GERMANIUM 
Degré d’ionisation IV 
À À 
AS 21/0 FES &p 2Paa 1 188,99 
48 28172 — 4D Pipe 1 229,81 
Degré d’ionisation III 
1 À - 
48? 185 — 4Sp 1P; 1 088,46 
ÉLÉMENT ÉTAIN 
Degré d’ionisation IV 
x À 
5$ 29112 — 9p Paye 1 314,45 
58 2$ ya — 5P Pa 1 437,63 
Degré d’ionisation IIL 
À À 
5s? 189, — 5sp 1P9 1 251,43 


Les études analytiques quantitatives concernant ces 
éléments ont été faites dans les conditions que l’on va 
indiquer. 

On a fabriqué des poudres à base d'oxyde de titane 
choisi à cause du soin apporté à sa préparation (taux 
d’impuretés 107) (*) et du nombre très réduit de ses 
- raies dans notre région spectrale. Ces poudres for- 
maient deux catégories que nous désignerons par À 
etB ;les poudres A contenant les éléments de la 4€ ligne 
de la classification (Ge-As-Se-Br) les poudres B ceux 
de la 5€ ligne (Sn-Sb-Te-I). La teneur de chacun des 
* 4 éléments contenus dans une poudre donnée était la 
même et les différentes poudres ont été faites par divi- 
sions successives et de façon à obtenir des concen- 
trations à intervalles réguliers avec une échelle loga- 
rithmique. Ces poudres ont été tassées dans le trou 
axial (3 mm de diamètre) d’électrodes en aluminium 
pur ne contenant aucun des éléments étudiés. 

Les mesures ont été faites sur des spectres 
de 100 étincelles obtenus avec la source à 3 électrodes 
décrite récemment [3] qui n'utilise qu’une électrode 
contenant la poudre (anode) au lieu de deux dans le 
cas de l’étincelle glissante, et dont le réglage est plus 
simple. L’anode d'aluminium était terminée par un 
cône de sorte que la décharge se produisait bien sur la 
poudre constituant le haut du cône. Cette décharge 
était celle d’un condensateur de 2 nf chargé sous une 
tension de 20 KV. 

On a obtenu une émission des éléments contenus 
dans la poudre, assez reproductible d’un spectre à 
l’autre, chaque dosage étant fait à partir d’une surface 
renouvelée, le bout de l’électrode contenant la poudre 
étant limé. Ceci permet de penser que le mode de fabri- 
cation des poudres était satisfaisant puisqu'il donne 
une homogénéité suffisante. On a étudié 10 poudres 


(*) Nous sommes heureux de remercier ici M. Lapluye 
du Laboratoire du Magnétisme de Bellevue qui nous a 
aimablement fourni l’oxyde de titane pur. 


ARSENIC SÉLÉNIUM BROME 
v VI VII 
967,69 844,15 736,09 

4 029,50 886,82 779,58 
IV V VI 
892,68 PECRE 661,05 

ANTIMOINE TELLURE IopE 

+ VI VII 

1 104,32 951,01 ? 

1 226,00 4 071,40 
IV V VI 

4 042,21 895,20 


dans chaque catégorie A et B, la poudre la plus con- 
centrée contenant 1 % de chaque élément, la moins 
concentrée 0,005 %,. Les courbes de dosage corres- 
pondent aux poudres les plus concentrées et ne com- 
portent que 8 ou 9 points pour rester dans les limites 
de densité optique imposées par le film photogra- 
phique. Sur la figure 1 on trouve les courbes corres- 
pondant aux poudres À, sur la figure 2 celles des 
poudres B ; elles sont obtenues en portant en abscisse 
le logarithme de la concentration et en ordonnées le 
logarithme Ô du rapport des déviations du micro- 
photomètre pour la raie de l’élément et pour la raie du 
du titane prise comme étalon. Nous n’avons pu iden- 
tifier de raies appartenant à l’iode dont les raies sen- 
sibles théoriques sont également absentes dans la 
littérature (voir tableau n° 1). Les courbes sont donc 
tracées pour chacun des 7 autres éléments et corres- 
pondent à la raie de plus courte longueur d’onde (raie 
la plus intense) du doublet des raies ultimes (courbe en 
trait plein) et à la raie la plus sensible émise par l'ion 
de degré d’ionisation inférieur d’une unité (courbe en 


trait pointillé) et indiquée aussi dans le tableau n° 1. 


La raie sensible de Sev n’est pas utilisable car elle est 
située dans un quintuplet très intense de Oy à 760 À. 
On a utilisé comme étalon interne la même raie du 
titane (indiquée sur la figure) pour le dosage de tous 
les éléments de chaque catégorie de poudres, afin de 
comparer la position et la pente respectives des diffé- 
rentes courbes de dosage. : We 


Les courbes de la figure 1 mettent en évidence la 


variation de la sensibilité de la méthode depuis le ger- 
manium jusqu’au brome. La courbe correspondant à la 
raie du doublet de raies ultimes (trait plein) se déplace 
vers les grandes valeurs de à c’est-à-dire les faibles 
intensités de raies de l’élément à doser depuis le ger- 
manium jusqu'au brome ; en même temps la pente 
diminue et la limite de concentration sur chaque courbe 
se déplace vers les fortes concentrations. Remarquons 
que la raie de degré d’ionisation inférieur d’une unité 


Li 


No 4 


est toujours moins intense que la raie ultime sauf pour 
le brome ce qui indique nettement que l’on arrive avec 
cet élément à des conditions défavorables pour l’exci- 


étalon Ti,, 128932 


2 3 4 
100 180. 330 560 1000 1800 3300.5600 10000 Cppm 


RTC 


tation des raies ultimes. On verra plus loin le rôle de 
l’énergie dépensée en regard de l’énergie d’ionisation et 
d’excitation de chaque raie. 

Les courbes de la figure 2 (poudres B) accusent une 
différence beaucoup moins nette entre les 3 éléments 
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étudiés dont les raies sensibles ont des intensités très 
voisines, permettant d’atteindre par suite, des concen- 
trations limites du même ordre. D’autre part, on 
n'obtient pas dans tous les essais le même rapport 
d'intensité entre la raie ultime et la raie de degré d'ioni- 
sation inférieur, dont la densité optique est très 
variable ; on peut avoir pour les 2 raies la même den- 


6 


09 étalon Ti,,, 776,82 


3 4 
180 330 560 1000 1800 3300 5600 10000 Cppm 
Free: 


2 
100 


sité optique, donnant des droites de pentes analogues 
comme pour l’étain et le tellure, ou pour la raie de 
degré d’ionisation inférieur une densité beaucoup plus 
faible et variant peu d’où une courbe non rectiligne 
comme celle de l’antimoine Sbiy sur la figure 2. Cette 


TABLEAU 2 


TENEURS LIMITES DÉCELABLES PAR PHOTOMÉTRIE PHOTOGRAPHIQUE 


ÉLÉMENT GERMANIUM 
Raie ultime x À 1 188,99 
Densité optique 0,19 
Concentration ppm 56 
Énergie e. v. 68 

ÉLÉMENT ÉTAIN 
Raie ultime à À 1 314,45 
Densité optique 0,17 
Concentration ppm 100 


Énergie e. v. 62 


ARSENIC SÉLÉNIUM BRroME 
987,69 844,15 661,05 
0,15 0,2 0,22 
56 100 330 
120 495 280 
ANTIMOINE TELLURE 
1 104,32 951,01 
0,22 0,11 
100 100 
108 470 


- 


>» 


répartition essentiellement variable de l’énergie entre 
les deux degrés d’ionisation indique que l’on favorise 
accidentellement quelquefois lémission de la raie 
ultime ce qui donne une courbe de pente plus forte et 
des points mieux répartis (courbe de dosage avec la 
raie de Sby). Ceci conduit à penser que l’on utilise une 
énergie instantanée par étincelle voisine d’une valeur 
critique pour l'excitation des raies indiquées. Il semble 
qu’en augmentant nettement cette énergie on favo- 
riserait l'émission desraies les plus ionisées, c’est-à-dire 
des raies ultimes au détriment de la raie de degré 
d’ionisation inférieur. En diminuant l'énergie Jes 
dosages pourraient être faits avec cette dernière (cas du 
brome figure 1) mais conduiraient à une sensibilité 
plus faible. 

Pour avoir une idée plus précise de la teneur limite 
décelable nous avons repris les dosages des poudres 
les moins concentrées en utilisant un temps de pose 
plus grand et toujours la raie ultime de plus courte 
longueur d’onde. Le tableau n° 2 indique les concen- 
trations limites pour chaque élément étudié ainsi que 
la densité optique correspondante pour cette teneur. 
Les résultats obtenus dans ces dosages portant sur des 
ions isoélectroniques à 1 électron sur la 4° ou 5° couche 
suggèrent de relier ces variations nettes de sensibilité 
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à l'énergie nécessaire pour obtenir les transitions inté- 
ressantes de tels ions. On à indiqué ces énergies en e. V. 
dans le tableau n° 2 également. Mais si ces différentes 
énergies d’ionisation sont en relation avec les diffé- 
rences de sensibilité d’un élément au suivant dans la 
classification, elles ne peuvent rendre compte du com- 
portement essentiellement différent des éléments de 
la 4e et de la 5° ligne. Rappelons en outre que l’on à 
obtenu avec les éléments de la 3° ligne une très bonne 
sensibilité (quelques ppm) pour le soufre et le phos- 
phore dans les poudres et pour le silicium dont le 


raies ultimes apparaissent dans tout métal commercial: : 


dit pur. La probabilité d’ionisation d’un élément 
donné semble donc être liée non seulement au nombre 
d'électrons de la couche externe, mais aussi au nombre 


total d'électrons. 
Lettre reçue le 2 mars 1959. 
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